Нахождение касательных

uses crt;
const nmax=88;
      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..1117*nmax] of Okruz;//массив окружностей
 
procedure Tochki(var t:MassP;var n:integer);//создание множества точек
var i:integer;
    f:text;
begin
  repeat
    write('Kolichestvo tochek ot 3 do ',nmax,' n=');
    readln(n);
  until n in [3..nmax];
for i:=1 to n do
  assign(f,'1.txt');
  reset(f);
for i:=1 to n do
  readln(f,t[i].x,t[i].y);
  close(f);
begin
  {t[i].x:=random(100);
  t[i].y:=random(100);}
 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;
   for i:=1 to k do begin
    writeln('Окружность ',[i],': ','X=',ok[i].o.x:5:2,'   ','Y=',ok[i].o.y:5:2,'   ','R=',ok[i].r:5:2);
 end;
end;
var t:MassP;
    ok:MassO;
    n,k,i,j:integer;
begin
clrscr;
randomize;
Tochki(t,n);
{можно для наглядности вывести координаты точек}
for i:=1 to n do begin
write('X= ',t[i].x:5:1,'  ');
end;
writeln;
for i:=1 to n do
write('Y= ',t[i].y:5:1,'  ');
writeln;
for i:=1 to n do
write;
MassOkr(t,n,ok,k);
 
for i:=1 to n-1 do 
for j:= i+1 to n do
begin
//проверка точки лежащей на окружности и нахождение касательной 
if (sqr(t[i].x-ok[i].o.x)+ sqr(t[i].y-ok[i].o.y)=sqr(ok[i].r) )
and  
(cos(pi/2)=(((-t[i].x+ok[i].o.x)*(-t[j].x+ok[i].o.x)+
             (-t[i].y+ok[i].o.y)*(-t[j].y+ok[i].o.y))/
             ((sqrt(sqr(-t[i].x+ok[i].o.x)+sqr(-t[i].y+ok[i].o.y))*
             (sqrt(sqr(-t[j].x+ok[i].o.x)+sqr(-t[j].y+ok[i].o.y)))))))
             then 
writeln('Координаты прямой =',[i],' ',t[i].x:5:1,':',t[i].y:5:1,' и',t[j].x:5:1,':',t[j].y:5:1)
else ;//writeln('Координаты прямой =',[i],' ',t[i].x:5:1,':',t[i].y:5:1,' и ',t[j].x:5:1,':',t[j].y:5:1);
end;
end.