множество точек А и арифметические ребусы

uses crt;
const n=20;
type point=record  {точка}
           x,y:real
           end;
     okruz=record  {окружность}
           x,y,r:real
           end;
function peres(p1,p2:point;o:okruz):boolean;{определение пересечения}
var s,ab,h:real;
begin
{площадь треугольника по точкам на прямой и центром окружности}
s:=abs(p1.x*(p2.y-o.y)+p2.x*(o.y-p1.y)+o.x*(p1.y-p2.y));
ab:=sqrt(sqr(p1.x-p2.x)+sqr(p1.y-p2.y));{длина основания=расстояние между точками}
h:=s/ab; {высота треугольника}
peres:=h<=o.r{примем что если прямая касается окружности=пересекает
            иначе знак = убрать}
end;
var p:array[1..2*n] of point;
    o:array[1..n] of okruz;
    i,j,k,per,mx,im,jm:byte;
begin
clrscr;
randomize;
for i:=1 to 2*n do
 begin
  p[i].x:=-9+18*random;
  p[i].y:=-9+18*random;
 end;
writeln('Координаты точек');
for i:=1 to 2*n do
write(i:2,'[',p[i].x:4:1,';',p[i].y:4:1,']   ');
writeln;
for i:=1 to n do
 begin
  o[i].x:=-9+random(18);
  o[i].y:=-9+random(18);
  o[i].r:=5+random(10);
 end;
writeln('Координаты центров и радиусы окружностей');
for i:=1 to n do
write(i:2,'[',o[i].x:4:1,';',o[i].y:4:1,';',o[i].r:5:1,'] ');
writeln;
mx:=0;
for i:=1 to 2*n-1 do{перебираем все пары точек}
for j:=i+1 to 2*n do
 begin
  per:=0;
  for k:=1 to n do{считаем пересечения их с окружностями}
  if peres(p[i],p[j],o[k]) then inc(per);
  if per>mx then
   begin
    mx:=per;{находим максимум и точки прямой}
    im:=i;
    jm:=j;
   end;
 end;
if mx=0 then write('Нет пересечений прямых и окружностей')
else
writeln('С максимальным числом окружностей пересекается прямая, проходящая через точки');
writeln(im:2,'[',p[im].x:4:1,';',p[im].y:4:1,']');
writeln(jm:2,'[',p[jm].x:4:1,';',p[jm].y:4:1,']');
readln
end.