Написать программу движения 3 тел Солнца, Земли и Марса

program kepler;
 
uses
  GraphABC;
 
const
  GM = 4 * sqr(pi);
 
var
  x, y, t, tk, dt, vx, vy, ax, ay, r, Hu, Hv, Sx, Sy: real;
 
procedure Init;
begin
  x := 1;
  y := 0;
  dt := 0.00001;
  tk := 2;
  vx := 0;
  vy := 2 * pi * 0.9;
  Sx := Window.Width / 3;
  Sy := Sx;
  Hu := Window.Center.X;
  Hv := Window.Center.Y;
  Line(0, round(Hv), Window.Width, round(Hv));
  Line(round(Hu), Window.Height, round(Hu), 0);
  Pen.Color := clgreen;
end;
 
procedure OutRes(x, y: real; c: color);
var
  v, u: word;
begin
  u := round(x * Sx + Hu);
  v := round(-y * Sy + Hv);
  PutPixel(u, v, c);
end;
 
procedure Euler;
begin
  while(t < tk) do
  begin
    r := sqrt(sqr(x) + sqr(y));
    ax := (-GM) / (sqr(r) * r) * x;
    ay := (-GM) / (sqr(r) * r) * y;
    vx := vx + ax * dt;
    vy := vy + ay * dt;
    x := x + vx * dt;
    y := y + vy * dt;
    t := t + dt;
    OutRes(x, y, clred);
  end;
end;
 
begin
  Init;
  OutRes(t, y, clred);
  Euler;
end.