Определить радиус и центр такой окружности, проходящей хотя бы через три разные точки множества точек

implementation
type Point=record
           x,y:double;
           end;
{$R *.dfm}
//определение параметров окружности
procedure Okr(a,b,c:Point;var o:Point;var r:double);
var ab,ac:Point;
    a1,b1,a2,b2,c1,c2:double;
begin
//середины сторон
ab.x:=(a.x+b.x)/2;
ab.y:=(a.y+b.y)/2;
ac.x:=(a.x+c.x)/2;
ac.y:=(a.y+c.y)/2;
//координаты центра
a1:=b.x-a.x;
b1:=b.y-a.y;
c1:=ab.x*(b.x-a.x)+ab.y*(b.y-a.y);
a2:=c.x-a.x;
b2:=c.y-a.y;
c2:=ac.x*(c.x-a.x)+ac.y*(c.y-a.y);
o.x:=(c1*b2-c2*b1)/(a1*b2-a2*b1);
o.y:=(a1*c2-a2*c1)/(a1*b2-a2*b1);
//радиус
r:=sqrt(sqr(a.x-o.x)+sqr(a.y-o.y));
end;
 
procedure TForm1.FormCreate(Sender: TObject);
begin
with StringGrid1 do
 begin
  colcount:=3;
  defaultcolwidth:=40;
  rowcount:=2;
  cells[0,0]:='№пп';
  cells[1,0]:='X';
  cells[2,0]:='Y';
  scrollbars:=ssVertical;
 end;
Memo1.Clear;
end;
procedure TForm1.Button1Click(Sender: TObject);
var a:array of Point;
    n,i,j,k,imn,jmn,kmn,q:integer;
    o,omn:Point;
    r,rmn:double;
begin
randomize;
val(Edit1.Text,n,k);
if(k<>0)or(n<3)then
 begin
  ShowMessage('Неверно введено количество точек');
  Edit1.Clear;
  Edit1.SetFocus;
  exit;
 end;
setlength(a,n);
StringGrid1.RowCount:=2;
for i:=0 to n-1 do
 begin
  a[i].X:=-10+20*random;
  a[i].Y:=-10+20*random;
  with StringGrid1 do
   begin
    cells[0,i+1]:=inttostr(i+1);
    cells[1,i+1]:=format('%.2f',[a[i].X]);
    cells[2,i+1]:=format('%.2f',[a[i].Y]);
    rowcount:=rowcount+1;
   end;
  end;
StringGrid1.RowCount:=StringGrid1.RowCount-1;
q:=0;
rmn:=0;
imn:=-1;
jmn:=-1;
kmn:=-1;
for i:=0 to n-3 do//перебираем точи
for j:=i+1 to n-2 do
for k:=j+1 to n-1 do
 begin
  q:=q+1;
  Okr(a[i],a[j],a[k],o,r);
  if q=1 then
   begin
    rmn:=r;
    omn.x:=o.x;
    omn.y:=o.y;
    imn:=i;
    jmn:=j;
    kmn:=k;
   end
  else if r<rmn then
   begin
    rmn:=r;
    omn.x:=o.x;
    omn.y:=o.y;
    imn:=i;
    jmn:=j;
    kmn:=k;
   end;
 end;
Memo1.Lines.Add('Окружность минимального радуса');
Memo1.Lines.Add('Радиус='+format('%.2f',[rmn]));
Memo1.Lines.Add('Координаты центра: x='+format('%.2f',[omn.x])+'  y='+format('%.2f',[omn.y]));
Memo1.Lines.Add('Образована точками c координатами:');
Memo1.Lines.Add('('+format('%.2f',[a[imn].x])+';'+format('%.2f',[a[imn].y])+')'+
              '  ('+format('%.2f',[a[jmn].x])+';'+format('%.2f',[a[jmn].y])+')'+
              '  ('+format('%.2f',[a[kmn].x])+';'+format('%.2f',[a[kmn].y])+')');
 
end;