Составить программу, демонстротирующию затухающие движения горизонтально брошенного мячика - Pascal ABC

uses graphABC;
const t=0.001; //точность сравнения
      n=20;//количество точек
type point=record
           x,y:integer;
           end;
function Plosh(a,b,c:point):real;//площадь треугольника по координатам
begin
Plosh:=abs(a.x*(b.y-c.y)+b.x*(c.y-a.y)+c.x*(a.y-b.y))/2;
end;
{проверка внутри или нет}
function Prin(a,b,c,d:point):boolean;
var s,s1,s2,s3:real;
begin
s:=Plosh(a,b,c);{площадь данного треугольника}
s1:=Plosh(a,b,d);{1 маленький}
s2:=Plosh(a,c,d);{2}
s3:=Plosh(b,c,d);{3}
Prin:=abs((s1+s2+s3)-s)<t;{если площади равны с заданной точностью}
end;
 
var m:array[1..n] of point;
    xc,yc,ms,i,j,k,p,imn,jmn,kmn,kv,kn,mn:integer;
    s1,s2:string;
begin
randomize;
{центр экрана}
xc:=windowwidth div 2;
yc:=windowheight div 2;
{создаем массив точек, чтобы не считать большие числа, возьмем
координаты по х=[-32..32] по y=[-24..24], потом умножим на 10 при рисовании}
for i:=1 to n do
 begin
  m[i].x:=random(xc div 5)-(xc div 10);
  m[i].y:=random(yc div 5)-(yc div 10);
 end;
{найдем вершины треугольника с минимальной разницей внутри и снаружи}
mn:=n;
for i:=1 to n-2 do
for j:=i+1 to n-1 do
for k:=j+1 to n do
if Plosh(m[i],m[j],m[k])>0 then //если треугольник
 begin
  kv:=0;
  kn:=0;
  for p:=1 to n do
  if not(p in [i,j,k]) then //вуршины треугольника не считаем
   begin
    if (Prin(m[i],m[j],m[k],m[p]))then kv:=kv+1
    else kn:=kn+1;
   end;
  if abs(kv-kn)<mn then
   begin
    mn:=abs(kv-kn);
    imn:=i;
    jmn:=j;
    kmn:=k;
   end;
 end;
ms:=10;{масштаб}
{нарисуем треугольник}
setpencolor(clGreen);
circle(xc+m[imn].x*10,yc-m[imn].y*10,4);
floodfill(xc+m[imn].x*10,yc-m[imn].y*10,clGreen);
circle(xc+m[jmn].x*10,yc-m[jmn].y*10,4);
floodfill(xc+m[jmn].x*10,yc-m[jmn].y*10,10);
circle(xc+m[kmn].x*10,yc-m[kmn].y*10,4);
floodfill(xc+m[kmn].x*10,yc-m[kmn].y*10,10);
line(xc+m[imn].x*10,yc-m[imn].y*10,xc+m[jmn].x*10,yc-m[jmn].y*10);
line(xc+m[jmn].x*10,yc-m[jmn].y*10,xc+m[kmn].x*10,yc-m[kmn].y*10);
line(xc+m[imn].x*10,yc-m[imn].y*10,xc+m[kmn].x*10,yc-m[kmn].y*10);
{нарисуем точки}
kv:=0;
kn:=0;
for i:=1 to n do
 begin
  if Prin(m[imn],m[jmn],m[kmn],m[i]) then {внутри}
   begin
    setpencolor(clBlue);
    kv:=kv+1;
   end
  else
   begin
    setpencolor(clRed);{снаружи}
    kn:=kn+1;
   end;
  circle(xc+m[i].x*10,yc-m[i].y*10,2);
 end;
str(kv-3,s1);//вершины треугольника не считаем
str(kn,s2);
textout(10,10,'Внутри='+s1+' снаружи='+s2);
end.