Определить радиус и центр окружности....

uses crt;
const nmax=28;
      e=0.01;
type
    Point=record {точка}
              x,y:real;{координаты}
            end;
    Okruz=record {окружность}
              o:Point;{координаты центра}
              r:real;{радиус}
            end;
    MassP=array[1..nmax] of Point; {множество-массив точек}
procedure Tochki(var t:MassP;var n:integer);{создание множества точек}
var i:integer;
begin
  repeat
  write('Количество точек 3 do ',nmax,' n=');
  readln(n);
  until n in [3..nmax];
  for i:=1 to n do
   begin
    t[i].x:=10*random;
    t[i].y:=10*random;
   end;
  write('X:');
  for i:=1 to n do
   write(t[i].x:5:2);
  writeln;
  write('Y:');
  for i:=1 to n do
   write(t[i].y:5:2);
  writeln;
end;
{вычисление параметров окружности, проходящей через три точки}
procedure Okr(t1,t2,t3:Point;var ok:Okruz);
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;
  ok.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;
        ok.o.x:=x;
  ok.r:=sqrt(sqr(t1.x-x)+sqr(t1.y-y));
end;
{определение принадлежности точки окружности}
function Prin(a:Point;ok:Okruz):boolean;
begin
if (sqr(a.x-ok.o.x)+sqr(a.y-ok.o.y))-sqr(ok.r)<-e then Prin:=true
  else Prin:=false
end;
var t:MassP;
    ok:Okruz;
    n,i,j,k,l,kol,mx:integer;
    mxr,mxx,mxy:real;
begin
  clrscr;
  randomize;
  Tochki(t,n);
  mx:=0;
  mxr:=0;
  mxx:=0;
  mxy:=0;
  for i:=1 to n-2 do
   for j:=i+1 to n-1 do
    for k:=j+1 to n do
     begin
      Okr(t[i],t[j],t[k],ok);
      kol:=0;
      for l:=1 to n do
       if Prin(t[l],ok)then kol:=kol+1;
       if kol>mx then
         begin
          mx:=kol;
          mxr:=ok.r;
          mxx:=ok.o.x;
          mxy:=ok.o.y;
         end;
     end;
  if mx=0 then write('Окружности, содержащей точки внутри нет!')
    else
      begin
       writeln('Окружность: координаты центра X=',mxx:0:2,' Y=',mxy:0:2,' R=',mxr:0:2);
       writeln('Внутри точек=',mx);
      end;
  readln;
end.