Моделирование упругого столкновения шаров

program Stolknovenie;
 
Uses 
  CRT, GraphABC;
var
  w1x, w2x, w1y, w2y, v1x, v2x, v1y, v2y, k, m1, m2: real;
  fi1, fi2, n1, n2: integer;
  gd, gr: integer;
  i, x1, x2, y1, y2, x, y, x10, x20, y10, y20, t: integer;
 
begin
  x10 := 30;
  y10 := 300;
  x20 := 600;
  y20 := 12;
  m1 := 4;
  m2 := 12;
  v1x := 1;
  v2x := 1;
  v1y := 0;
  v2y := 1;
  fi1 := 0;
  k := 0;
  fi2 := trunc(45 / 57);
  n1 := 1;
  n2 := 1;
  y := 0;
  x := 0;
  delay(10);
  y1 := 0;
  t := 0;
  while (abs((x10 + x1) - (x20 - x2)) > 8) or (abs((y10 + y1) - (y20 + y2)) > 8) do
  begin
    x1 := trunc(v1x * t);
    x2 := trunc(v2x * t);
    y1 := trunc(v1y * t);
    y2 := trunc(v2y * t);
    t := t + 1;
    SetPenWidth(1);
    SetPenColor(clBlue);
    Circle(x10 + x1, y10 + y1, 8);
    SetPenColor(clGreen);
    Circle(x20 - x2, y20 + y2, 8);
    delay(10);
    SetPenWidth(6);
    SetPenColor(clWhite);
    Circle(x10 + x1, y10 + y1, 8);
    Circle(x20 - x2, y20 + y2, 8);
  end;
  x10 := x10 + x1;
  y10 := y10 + y;
  x20 := x20 - x2;
  y20 := y20 + y2;
  t := 0;
  for i := 1 to 300 do
  begin
    t := t + 1;
    W1x := V1x + M2 * (1 + k) / (M1 + M2) * (abs(V1x) * cos(fi1) + abs(V2x) * cos(fi2)) * n1;
    W2x := V2x + M1 * (1 + k) / (M1 + M2) * (abs(V1x) * cos(fi1) + abs(V2x) * cos(fi2)) * n2;
    W1y := V1y + M2 * (1 + k) / (M1 + M2) * (abs(V1y) * cos(fi1) + abs(V2y) * cos(fi2)) * n1;
    W2y := V2y + M1 * (1 + k) / (M1 + M2) * (abs(V1y) * cos(fi1) + abs(V2y) * cos(fi2)) * n2;
    x1 := -trunc(w1x * t);
    x2 := -trunc(w2x * t);
    y1 := -trunc(w1y * t);
    y2 := -trunc(w2y * t);
    SetPenColor(clBlue);
    SetPenWidth(1);
    Circle(x10 + x1, y10 + y1, 8);
    SetPenColor(clGreen);
    Circle(x20 - x2, y20 + y2, 8);
    delay(10);
    SetPenWidth(6);
    SetPenColor(clWhite);
    Circle(x10 + x1, y10 + y1, 8);
    Circle(x20 - x2, y20 + y2, 8);
  end;
  Readkey;
end.