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

for i:=1 to n do
  for j:=1 to n do
    if (j<>i) then
    begin
//вот тут нам надо найти расстояние между точками зная их координаты
//Теперь у нас есть диаметр радиуса, который представлен прямой от x1,y1 до x2,y2 
//находим центр этой прямой, находим точки x0 и y0 нашей окружности пробегаемся ещё раз по
//всему циклу и просматриваем все точки не входят ли они в нашу окружность (нам известны 
//координаты центра, радиус (Диаметр пополам) и координаты проверяемой точки (не забываем
//исключить точки с индексом i и j. Если нашли хотя бы одну точку (входящую в радиус),
//останавливаем цикл, сравниваем радиус с тем, что у нас уже есть 
//(если это первый радиус, то просто его запоминаем ну и позицию центра не забываем 
//сохранять). И продолжаем цикл.
    end;
//После выполнения всего цикла выводим радиус и его центр.

uses crt;
const nmax=15;{макс. кол. точек}
      e=0.0001;{точность определения}
type
Point=record{тип точка}
      x,y:real;
      end;
MassP=array[1..nmax] of Point;{массив точек}
procedure Tochki(var t:MassP;var n:integer);{создание массива точек}
var i:integer;
begin
repeat
write('Количество точек от 3 до ',nmax,' n=');
readln(n);
until n in [3..nmax];
for i:=1 to n do
 begin
  t[i].x:=20*random-10;
  t[i].y:=20*random-0;
 end;
writeln('Координаты точек:');
write('X:');
for i:=1 to n do
write(t[i].x:5:1);
writeln;
write('Y:');
for i:=1 to n do
write(t[i].y:5:1);
writeln;
end;
{вычисление координат центра и радиуса окружности}
procedure Okr(t1,t2,t3:Point;var o:Point;var r:real);
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;
 
 
var t:MassP;
    rc:Point;
    n,i,j,k,imn,jmn,kmn:integer;
    r,mn:real;
begin
clrscr;
randomize;
Tochki(t,n);
Okr(t[1],t[2],t[3],rc,r);
mn:=r;{пока минимальный в точках 1,2,3}
{ищем самый маленький радиус}
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],rc,r);
  if r<mn then
   begin
    mn:=r;
    imn:=i;
    jmn:=j;
    kmn:=k;
   end;
 end;
writeln('Окружность c минимальным радиусом:');
writeln(' координаты центра x0=',rc.x:0:2,' Y0=',rc.y:0:2);
writeln(' радиус R=',r:0:2);
writeln(' номера точек ',imn,' ',jmn,' ',kmn);
readln
end.