Из множества координат найти вершины прямоугольника

uses crt;
 
function S(x1,x2,x3,x4,y1,y2,y3,y4:real):real;
var a, b, c, d ,p: real;
    t:text;
begin
 
a:=sqrt(sqr(abs(x2-x1))+sqr(abs(y2-y1)));
b:=sqrt(sqr(abs(x3-x2))+sqr(abs(y3-y2)));
S:=a*b;
 
end;
 
var
  sum, min:real;
  x:array[1..1000] of real;
  y:array[1..1000] of real;
   n,i,j,l,z:integer;
  indexi, indexj, indexl,indexz: integer;
begin
  clrscr;
  write('Kol-vo tochek: ');
  readln(n);
  min:=maxint;
  randomize;
  for i:=1 to n do begin
    X[i]:=random(21)-10;
    Y[i]:=random(21)-10;
  end;
  for i:=1 to n do begin
       for j:=1 to n do begin 
      if (i=j) then continue;
      for l:=1 to n do begin
      if (l=j) or (l=i) then continue;
       for z:=1 to n do begin
       if (z=j) or (z=i) or (z=l) then continue;
        sum := S(X[i], X[j], X[l], X[z], Y[i], Y[j], Y[l], Y[z]);
        if sum = 0 then continue;
       if sum<min then begin
          min:=Sum;
          indexi := i; indexj := j; indexl := l;indexz := z;
        end;
      end;
    end;
  end;
  writeln('min S: ',min:0:2, ' точки: ',
    '[', indexi, '](', x[indexi]:0:0, ',', y[indexi]:0:0, ')',
    '[', indexj, '](', x[indexj]:0:0, ',', y[indexj]:0:0, ')',
    '[', indexl, '](', x[indexl]:0:0, ',', y[indexl]:0:0, ')',
    '[', indexz, '](', x[indexz]:0:0, ',', y[indexz]:0:0, ')');
  readln;
end;
end.

uses crt;
const p=0.0001;{точность сравнения длин}
      m=6;{кол. отрезков=4 стороны+2 диагонали}
      max=16;{максимальное число точек}
type Point=record {тип - точка}
           x,y:real;
           end;
     Mas=array[1..100] of Point;
{вычисление длины отрезка}
function Dlina(a,b:Point):real;
begin
Dlina:=sqrt(sqr(a.x-b.x)+sqr(a.y-b.y));
end;
function IsRect(a,b,c,d:Point):boolean;
var r:array[1..6] of real;
    i,j:integer;
    x:real;
begin
r[1]:=Dlina(a,b);
r[2]:=Dlina(a,c);
r[3]:=Dlina(a,d);
r[4]:=Dlina(b,c);
r[5]:=Dlina(b,d);
r[6]:=Dlina(c,d);
for i:=1 to m-1 do {сортируем длины по возрастанию}
for j:=i+1 to m do
if r[i]>r[j] then
 begin
  x:=r[i];
  r[i]:=r[j];
  r[j]:=x;
 end;
{если первые 4 попарно равны и последние 2 равны и они больше первых}
IsRect:=((abs(r[1]-r[2])<p)and(abs(r[3]-r[4])<p))
and(abs(r[5]-r[6])<p)and(r[5]>r[4])
end;
 
var t:array[1..100] of Point;{массив точек}
    n,i,j,k,l,c:integer;
begin
clrscr;
randomize;
repeat
writeln('Введите количество точек от 4 до ',max,' n=');
readln(n);
until n in [1..16];
for i:=1 to n do
 begin
  writeln('Точка ',i);
  write('x=');readln(t[i].x);{вводить лучше с клавы, чтобы были прямоугольники}
  write('y=');readln(t[i].y)
 end;
clrscr;
writeln('Координаты:');
write('№');
for i:=1 to n do
write(i:4);
writeln;
write('X:');
for i:=1 to n do
write(t[i].x:4:1);
writeln;
write('Y:');
for i:=1 to n do
write(t[i].y:4:1);
writeln;
writeln('Вершины прямоугольников');
c:=0;
for i:=1 to n-3 do
for j:=i+1 to n-2 do
for k:=j+1 to n-1 do
for l:=k+1 to n do
if IsRect(t[i],t[j],t[k],t[l])then
 begin
  c:=1;
  writeln('(',t[i].x:4:1,t[i].y:4:1,') (',t[j].x:4:1,t[j].y:4:1,')'+
         ' (',t[k].x:4:1,t[k].y:4:1,') (',t[l].x:4:1,t[l].y:4:1,')');
 end;
if c=0 then write('Прямоугольников нет');
readln
end.