Скорость и траектория при двухмерном движении

uses graphABC;
//константы по условию
const tk=6;
      t1=2;
      x1=2;
      y1=6;
var a:array[1..4,0..tk]of real; //массив скоростей и координат
    t,x0,y0:integer;
    xmax,xmin,ymax,ymin,mx,my:real;
    s1,s2:string;
begin
//вычислим и заполним скорости в точках
for t:=0 to tk do
 begin
  a[1,t]:=2+20/(t*t+4);
  a[2,t]:=6+60/(t*t+36);
 end;
//координаты в точке t1(x1,y1)
a[3,2]:=x1;
a[4,2]:=y1;
//вычислим координаты слева от t1
for t:=1 downto 0 do
 begin
  a[3,t]:=a[3,t+1]-a[1,t];
  a[4,t]:=a[4,t+1]-a[2,t];
 end;
//вычислим координаты справа от t1
for t:=3 to tk do
 begin
  a[3,t]:=a[3,t-1]+a[1,t];
  a[4,t]:=a[4,t-1]+a[2,t];
 end;
//найдем минимумы, максимумы и выведем данные на экран
xmax:=a[3,0];
xmin:=a[3,0];
ymax:=a[4,0];
ymin:=a[4,0];
for t:=0 to tk do
 begin
  writeln('t=',t,' vx=',a[1,t]:7:3,' x=',a[3,t]:7:3,' vy=',a[2,t]:7:3,' y=',a[4,t]:7:3);
  if a[3,t]>xmax then xmax:=a[3,t];
  if a[3,t]<xmin then xmin:=a[3,t];
  if a[4,t]>ymax then ymax:=a[4,t];
  if a[4,t]<ymin then ymin:=a[4,t];
 end;
setwindowsize(500,500);
centerwindow;
mx:=(windowwidth-20)/(xmax-xmin);
x0:=round(windowwidth*(-xmin)/(xmax-xmin));
my:=(windowheight-140)/(ymax-ymin);
y0:=round(windowheight*(ymax)/(ymax-ymin))-100;
line(0,y0,windowwidth,y0);
line(x0,0,x0,windowwidth-100);
for t:=1 to trunc(xmax)+1 do
 begin
  line(x0+round(t*mx),y0-3,x0+round(t*mx),y0+3);
  line(x0-round(t*mx),y0-3,x0-round(t*mx),y0+3);
  textout(x0+round(t*mx),y0+10,inttostr(t));
  textout(x0-round(t*mx),y0+10,inttostr(-t));
 end;
for t:=1 to trunc(ymax)+1 do
if t mod 2=0 then
 begin
  line(x0-3,y0-round(t*my),x0+3,y0-round(t*my));
  textout(x0-20,y0-round(t*my),inttostr(t));
  if t<11 then
   begin
    line(x0-3,y0+round(t*my),x0+3,y0+round(t*my));
    textout(x0-30,y0+round(t*my),inttostr(-t));
   end;
 end;
setpencolor(clBlue);
for t:=0 to tk do
if t=0 then moveto(x0+round(a[3,t]*mx),y0-round(a[4,t]*my))
else lineto(x0+round(a[3,t]*mx),y0-round(a[4,t]*my));
setbrushcolor(clRed);
for t:=0 to tk do
circle(x0+round(a[3,t]*mx),y0-round(a[4,t]*my),3);
setbrushstyle(bsClear);
str(xmin:0:3,s1);
str(xmax:0:3,s2);
textout(50,windowheight-80,'Минимальная координата х='+s1+' максимальная='+s2);
str(ymin:0:3,s1);
str(ymax:0:3,s2);
textout(50,windowheight-60,'Минимальная координата y='+s1+' максимальная='+s2);
 
end.