В заданном множестве точек найти координаты вершин всех треугольников

{$R *.dfm}
//определение площади треугольника по координатам вершин
function Plosh(a,b,c:Tpoint):real;
begin
Plosh:=abs((a.x-c.x)*(b.y-c.y)-(b.x-c.x)*(a.y-c.y))/2;
end;
//проверка точки на принадлежность треугольнику
function Prin(a,b,c,d:Tpoint):boolean;
var s,s1,s2,s3:real;
begin
s:=Plosh(a,b,c);{площадь данного треугольника}
if s<t then result:=false
else
 begin
  s1:=Plosh(a,b,d);{1 маленький}
  s2:=Plosh(a,c,d);{2}
  s3:=Plosh(b,c,d);{3}
  result:=abs(s1+s2+s3-s)<t;//если площади равны с заданной точностью
 end;
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.ScrollBars:=ssVertical;
end;
procedure TForm1.Button1Click(Sender: TObject);
var a:array of TPoint;
    n,i,j,k,p,q,r:integer;
    f:boolean;
begin
randomize;
val(Edit1.Text,n,k);
if(k<>0)or(n<6)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+random(21);
  a[i].Y:=-10+random(21);
  with StringGrid1 do
   begin
    cells[0,i+1]:=inttostr(i+1);
    cells[1,i+1]:=inttostr(a[i].X);
    cells[2,i+1]:=inttostr(a[i].Y);
    rowcount:=rowcount+1;
   end;
  end;
StringGrid1.RowCount:=StringGrid1.RowCount-1;
Memo1.Clear;
Memo1.Lines.Add('Пары треугольников, входящих друг в друга');
f:=false;
for i:=0 to n-3 do//перебираем точи для 1 треугольника
for j:=i+1 to n-2 do
for k:=j+1 to n-1 do
for p:=0 to n-3 do//то же ддя второго
for q:=p+1 to n-2 do
for r:=q+1 to n-1 do
if not(p in [i,j,k])and not(q in [i,j,k])and not(r in [i,j,k])//если не совпадают
//и точки 2 принадлежат первому
and Prin(a[i],a[j],a[k],a[p])and Prin(a[i],a[j],a[k],a[q])
and Prin(a[i],a[j],a[k],a[r])then
 begin
  f:=true;
  Memo1.Lines.Add('[('+inttostr(a[i].X)+';'+inttostr(a[i].Y)+');'+
                   '('+inttostr(a[j].X)+';'+inttostr(a[j].Y)+');'+
                   '('+inttostr(a[k].X)+';'+inttostr(a[k].Y)+')]'+
                   '   [('+inttostr(a[p].X)+';'+inttostr(a[p].Y)+');'+
                   '('+inttostr(a[q].X)+';'+inttostr(a[q].Y)+');'+
                   '('+inttostr(a[r].X)+';'+inttostr(a[r].Y)+')]');
  end;
if not f then Memo1.Lines.Add('Таких пар нет');
end;

const t=0.001;//точность сравнения вещественных чисел

if s<0.001 then result:=false

{$R *.dfm}
//проверка принадлежности точки треугольнику
function prin(t1,t2,t3,t:TPoint):boolean;
//положение точки относительно прямой, проходящей через 2 точки
function odin(tt1,tt2,tt:TPoint):integer;
begin
result:=tt.X*(tt2.Y-tt1.Y)+tt.Y*(tt1.X-tt2.X)+tt1.Y*tt2.X-tt1.X*tt2.Y;
end;
var q1,q2,q3:integer;
begin
q1:=odin(t1,t2,t);
q2:=odin(t2,t3,t);
q3:=odin(t3,t1,t);
result:=((q1>0)and(q2>0)and(q3>0))
or((q1<0)and(q2<0)and(q3<0));
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.ScrollBars:=ssVertical;
end;
procedure TForm1.Button1Click(Sender: TObject);
var a:array of TPoint;
    n,i,j,k,p,q,r:integer;
    f:boolean;
begin
randomize;
val(Edit1.Text,n,k);
if(k<>0)or(n<4)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+random(21);
  a[i].Y:=-10+random(21);
  with StringGrid1 do
   begin
    cells[0,i+1]:=inttostr(i+1);
    cells[1,i+1]:=inttostr(a[i].X);
    cells[2,i+1]:=inttostr(a[i].Y);
    rowcount:=rowcount+1;
   end;
  end;
StringGrid1.RowCount:=StringGrid1.RowCount-1;
Memo1.Clear;
Memo1.Lines.Add('Пары треугольников, входящих друг в друга');
f:=false;
for i:=0 to n-3 do//перебираем точи для 1 треугольника
for j:=i+1 to n-2 do
for k:=j+1 to n-1 do
for p:=0 to n-3 do//то же ддя второго
for q:=p+1 to n-2 do
for r:=q+1 to n-1 do
if not(p in [i,j,k])and not(q in [i,j,k])and not(r in [i,j,k])//если не совпадают
and prin(a[i],a[j],a[k],a[p])//и принадлежат точки 2-го 1-му
and prin(a[i],a[j],a[k],a[q])
and prin(a[i],a[j],a[k],a[r]) then
 begin
  f:=true;
  Memo1.Lines.Add('[('+inttostr(a[i].X)+';'+inttostr(a[i].Y)+');'+
                   '('+inttostr(a[j].X)+';'+inttostr(a[j].Y)+');'+
                   '('+inttostr(a[k].X)+';'+inttostr(a[k].Y)+')]'+
                   '   [('+inttostr(a[p].X)+';'+inttostr(a[p].Y)+');'+
                   '('+inttostr(a[q].X)+';'+inttostr(a[q].Y)+');'+
                   '('+inttostr(a[r].X)+';'+inttostr(a[r].Y)+')]');
  end;
if not f then Memo1.Lines.Add('Таких пар нет');
end;