Найти касательные к окружностям

uses crt;
const nmax=59;//больше точек при переборе не потянет
      e=0.01;//точность вычислений
type
Point=record//тип точка с координатами
      x,y:real;
      end;
Okruz=record//тип окружность, координаты центра и радиус
      o:Point;
      r:real;
    end;
MassP=array[1..nmax] of Point;//массив точек
MassO=array[1..117*nmax] of Okruz;//массив окружностей
procedure Tochki(var t:MassP;var n:integer);//создание множества точек
var i:integer;
begin
repeat
write('Kolichestvo tochek ot 3 do ',nmax,' n=');
readln(n);
until n in [3..nmax];
for i:=1 to n do
 begin
  t[i].x:=20*random-0;
  t[i].y:=20*random-0;
 end;
end;
 
procedure Okr(t1,t2,t3:Point;var o:Point;var r:real);//вычисление параметров окружности, проходящей через 3 точки
var a,b,x,y:real;
    k0,k1,k2,m0,m1,m2:real;
begin
k0:=sqr(t1.x)-sqr(t2.x)+sqr(t1.y)-sqr(t2.y);
k1:=2*(t1.y-t2.y);
k2:=2*(t1.x-t2.x);
m0:=sqr(t1.x)-sqr(t3.x)+sqr(t1.y)-sqr(t3.y);
m1:=2*(t1.y-t3.y);
m2:=2*(t1.x-t3.x);
a:=k2*m0-k0*m2;
b:=k2*m1-k1*m2;
if b=0 then exit;
y:=a/b;
o.y:=y;
if abs(m2) > e then x:=(m0-y*m1)/m2
else
  if abs(k2) > e then x:=(k0-y*k1)/k2
  else exit;
o.x:=x;
r:=sqrt(sqr(t1.x-x)+sqr(t1.y-y));
end;
procedure MassOkr(t:MassP;n:integer;var ok:MassO;var k:integer);//создание массива окружностей, 
//проходящих через любые 3 точки множества
var i,j,l:integer;
begin
k:=0;
for i:=1 to n-2 do
for j:=i+1 to n-1 do
for l:=j+1 to n do
 begin
  k:=k+1;
  Okr(t[i],t[j],t[l],ok[k].o,ok[k].r);
 end;
end;
function Prin(a,b:Point;c:real):boolean;//определение принадлежности точки кругу
begin
if (sqr(a.x-b.x)+sqr(a.y-b.y))-c*c<-e then Prin:=true
else Prin:=false
end;
//нахождение параметров окружности, включающей максимальное число точек
procedure Max(t:MassP;n:integer;ok:MassO;k:integer;var ro:Okruz;var kv,kn:integer);
var i,j,mx,k1:integer;
begin
mx:=0;ro.o.x:=0;ro.o.y:=0;ro.r:=0;
for i:=1 to k do//перебираем все окружности
 begin
  k1:=0;//количество точек внутри
  for j:=1 to n do//перебираем все точки
  if Prin(ok[i].o,t[j],ok[i].r) then k1:=k1+1;//если внутри окружности, считаем
  if k1>mx then
   begin
    mx:=k1;//находим максимальное
    ro.o:=ok[i].o;//координаты центра
    ro.r:=ok[i].r;//радиус
   end;
 end;
kv:=mx;//количество внутри
kn:=n-kv-3;//количество снаружи(3 точки на окружности)
end;
 
var t:MassP;
    ok:MassO;
    ro:Okruz;
    n,k,i,vn,sn:integer;
    gd,gm,x0,y0:integer;
    s1,s2:string;
begin
clrscr;
randomize;
Tochki(t,n);
{можно для наглядности вывести координаты точек}
for i:=1 to n do
write(t[i].x:5:2);
writeln;
for i:=1 to n do
write(t[i].y:5:2);
writeln;
MassOkr(t,n,ok,k);
Max(t,n,ok,k,ro,vn,sn);
for i:= 1 to n do 
if abs(t[i].x)=ro.r then writeln('Касательная в точке',t[i].x)
else writeln('Касательных нет'); 
writeln('Окружность: координаты центра x0=',ro.o.x:0:2,' Y0=',ro.o.y:0:2,' R=',ro.r:0:2);
writeln('Внутри точек=',vn,'  снаружи=',sn);
readln;
end.