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

const e=0.001;
type point=record
           x,y:real;
           end;
{вычисление координат центра и радиуса окружности
алгоритм отсюда http://algolist.manual.ru/maths/geom/equation/circle.php
вертикальные линии вычислил напрямую без обмена точек.}
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;
o.y:=a/b;
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;
{Принадлежность точки кругу}
function Prin(o:point;r:real;t:point):boolean;
begin
Prin:=sqr(o.x-t.x)+sqr(o.y-t.y)<=r*r; //квадрат расстояния между точкой и центром
                                      //не больше квадрата радиуса
end;
var a,b:array[1..20] of point;
    n,m,i,j,k,q,v,imn,jmn,kmn:integer;
    o:point;
    r,mn:real;
begin
repeat
write('Введите количество точек 1 множества от 3 до 20 n=');
readln(n);
until n in [3..20];
repeat
write('Введите количество точек 2 множества от 3 до 20 m=');
readln(m);
until m in [3..20];
randomize;
writeln('Первое множество');
for i:=1 to n  do
 begin
  a[i].x:=-9+18*random; //точки в квадрате (-9..9,-9..9))
  a[i].y:=-9+18*random;
  write(i:2,'(',a[i].x:4:1,';',a[i].y:4:1,') ');
  if i mod 5=0 then writeln;
 end;
writeln;
writeln('Второе множество');
for i:=1 to m  do
 begin
  b[i].x:=-7+18*random; //сдвинем на 2 вправо
  b[i].y:=-7+18*random; //и на 2 вверх
  write(i:2,'(',b[i].x:4:1,';',b[i].y:4:1,') ');
  if i mod 5=0 then writeln;
 end;
writeln;
v:=0;//пока нет нужной окружности
imn:=0;
jmn:=0;
kmn:=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(a[i],a[j],a[k],o,r);
  q:=1;
  while(q<=m)and prin(o,r,b[q])do inc(q);
  if q>m then
   begin
    if v=0 then
     begin
      mn:=r;
      imn:=i;
      jmn:=j;
      kmn:=k;
      v:=1;
     end
    else if r<mn then
     begin
      mn:=r;
      imn:=i;
      jmn:=j;
      kmn:=k;
     end;
   end;
 end;
if v=0 then write('Нет точек, образующих заданную окружность')
else writeln('Заданная окружность проходит через точки:');
writeln(imn:2,'(',a[imn].x:4:1,';',a[imn].y:4:1,')');
writeln(jmn:2,'(',a[jmn].x:4:1,';',a[jmn].y:4:1,') ');
writeln(kmn:2,'(',a[kmn].x:4:1,';',a[kmn].y:4:1,') ');
end.