Составить программу, моделирующую абсолютно упругий центральный удар двух шаров с массами m1 и m2 движущихся со скоростями v1 и v2

uses crt, graph;
var t1,t2,y:integer;
    gd,gm:integer;
begin
initgraph(t1,t2,'');
t1:=25;
t2:=getmaxX-25;
y:=getmaxY div 2;
circle(t1,y,20);
circle(t2,y,20);
while t1<=getmaxX div 2-20 do
 begin
  delay(50);
  setcolor(black);
  circle(t1,y,20);
  circle(t2,y,20);
  t1:=t1+5;
  t2:=t2-5;
  setcolor(white);
  circle(t1,y,20);
  circle(t2,y,20);
 end;
while t1>=25 do
 begin
  delay(50);
  setcolor(black);
  circle(t1,y,20);
  circle(t2,y,20);
  t1:=t1-5;
  t2:=t2+5;
  setcolor(white);
  circle(t1,y,20);
  circle(t2,y,20);
 end;
outtextxy(250,460,'Press any key...');
readkey
end.

Uses Dos, Crt, Graph;
  Var
    gd, gm,
    Yc,
    ScrWidth, x10, x20, x1, x2, r: Integer;
    s,
    Color                        : Word;
    TStep, Scale,
    v10, v20,
    v1, v2, v11, v21, m1, m2, d  : Single;
    p1, p2, px                   : Pointer;
    KK                           : Char;
    impact                       : Boolean;
Begin
  gd:=detect;
  InitGraph(gd,gm,'d:\bp\bgi');
  ScrWidth:=GetMaxX;
  Yc:=GetMaxY shr 1;
 
  Color:=Yellow;
  SetColor(Color);
  SetFillStyle(1,Color);
  r:=12;
{массы, кг}
  m1:=12; m2:=8;
{начальные скорости, м/с}
  v10:=15; v20:=10;
 
  Line(4,Yc+r,ScrWidth,Yc+r);
  OutTextXY(30,Yc-120,'[a] - догонять; [b] - навстречу.');
  x10:=25; x20:=ScrWidth div 4;
 
 
  kk:=ReadKey;
  If kk='b' then
    begin
      x10:=25; x20:=ScrWidth - 25;
      v10:=15; v20:=-10;
    end;
 
  d:=ScrWidth*0.1;
  Scale:=ScrWidth/d;
 
  s:=ImageSize(0,0,2*r,2*r);
  GetMem(p1,s);
  GetMem(p2,s);
  GetMem(px,s);
 
  GetImage(x10-r,Yc-r,x10+r,Yc+r,px^);
  PieSlice(x10,Yc,0,360,r);
  GetImage(x10-r,Yc-r,x10+r,Yc+r,p1^);
 
  PieSlice(x20,Yc,0,360,r);
  GetImage(x20-r,Yc-r,x20+r,Yc+r,p2^);
  OutTextXY(30,Yc+120,'[1] - старт; [Esc] - выход.');
 
  TStep:=0.05;  {шаг по времени, сек}
  kk:=ReadKey;
  If kk='1' then
    Repeat
      x1:=x10; x2:=x20;
      v1:=v10; v2:=v20;
      impact:=false;
      kk:=#2;
      Repeat
        PutImage(x1-r,Yc-r,px^,CopyPut);
        PutImage(x2-r,Yc-r,px^,CopyPut);
        x1:=x1 + Round(v1*TStep*Scale);
        x2:=x2 + Round(v2*TStep*Scale);
        if not impact then
          if x1+r>=x2-r then  {момент удара. Здесь скорости изменятся в соответствии с законами сохранения}
            begin
              v11:=(2*m2*v2 + (m1-m2)*v1)/(m1+m2);
              v21:=(2*m1*v1 + (m2-m1)*v2)/(m1+m2);
              v1:=v11; v2:=v21;
              impact:=true;
            end;
        PutImage(x1-r,Yc-r,p1^,CopyPut);
        PutImage(x2-r,Yc-r,p2^,CopyPut);
        Delay(Round(TStep*1000));
        If KeyPressed then
          kk:=ReadKey;
      Until (kk=#27) or (x1+r>=ScrWidth) or (x2+r>=ScrWidth);
      if kk=#2 then kk:=ReadKey;
      PutImage(x1-r,Yc-r,px^,CopyPut);
      PutImage(x2-r,Yc-r,px^,CopyPut);
    Until kk<>'1';
  CloseGraph;
End.