Найти координаты трёх точек, которые образуют треугольник с наименьшей площадью

  uses crt;
  procedure SquareTR(a,b,c,p,x1,x2,x3,y1,y2,y3,S:real);
begin
  a:=sqrt(sqr(abs(x1-x2))+sqr(abs(y1-y2)));
  b:=sqrt(sqr(abs(x2-x3))+sqr(abs(y2-y3)));
  c:=sqrt(sqr(abs(x3-x1))+sqr(abs(y3-y1)));
  p:=(a+b+c)/2;
  S:=sqrt(p*(p-a)*(p-b)*(p-c));
end;
  var X:array[1..100] of real;
      Y:array[1..100] of real;
      a,b,c,p,x1,x2,x3,y1,y2,y3,S,min:real;
      N,i,j,l:integer;
begin
  clrscr;
  write('Kol-vo tochek: ');
  readln(N);
  min:=0;
  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
    x1:=X[i];
    y1:=Y[i];
    for j:=1 to N do begin
      if j<>i then begin
        x2:=X[j];
        y2:=Y[j];
      end;
      for l:=1 to N do begin
        if (l<>i) and (l<>i) then begin
          x3:=X[l];
          y3:=Y[l];
        end;
        if S<min then min:=S;
      end;
    end;
  end;
  writeln('min S: ',min);
  readln
end.

uses crt;
 
function S(x1,x2,x3,y1,y2,y3:real):real;
var a, b, c, p : real;
begin
  a:=sqrt(sqr(abs(x1-x2))+sqr(abs(y1-y2)));
  b:=sqrt(sqr(abs(x2-x3))+sqr(abs(y2-y3)));
  c:=sqrt(sqr(abs(x3-x1))+sqr(abs(y3-y1)));
  p:=(a+b+c)/2;
  S:=sqrt(p*(p-a)*(p-b)*(p-c));
end;
 
var
  X:array[1..100] of real;
  Y:array[1..100] of real;
  sum, min:real;
  N,i,j,l,c:integer;
  indexi, indexj, indexl : integer;
begin
  clrscr;
  write('Kol-vo tochek: ');
  readln(N);
  min:=maxlongint;
  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 = i) or (l = j) then continue; { пропускаем одинаковые индексы }
 
        sum := s(X[i], X[j], X[l], Y[i], Y[j], Y[l]);
        if sum = 0 then continue; { пропускаем нулевую площадь: точки расположены на одной линии }
 
        if sum<min then begin
          min:=Sum;
          indexi := i; indexj := j; indexl := l;
        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, ')');
  readln
end.