Если точки могут быть вершинами равнобедренного тупоугольного треугольника, вычислите его площадь

uses crt;
const t=0.001;{точнсть сравнения чисел}
{вычисление длины стороны}
function D(x1,y1,x2,y2:real):real;
begin
D:=sqrt(sqr(x1-x2)+sqr(y1-y2));
end;
{вычисление площади}
function S(x1,y1,x2,y2,x3,y3:real):real;
begin
s:=abs((x1-x3)*(y2-y3)-(x2-x3)*(y1-y3))/2;
end;
 
var x1,y1,x2,y2,x3,y3,pl:real;
    a:array[1..4] of real;
    i,j:byte;
    x:real;
    st:string;
begin
clrscr;
writeln('Введите координаты вершин треугольника:');
writeln('1->');
readln(x1,y1);
writeln('2->');
readln(x2,y2);
writeln('3->');
readln(x3,y3);
a[1]:=D(x1,y1,x2,y2);{найдем длины сторон и запишем в массив}
a[2]:=D(x1,y1,x3,y3);
a[3]:=D(x3,y3,x2,y2);
pl:=S(x1,y1,x2,y2,x3,y3);{найдем площадь}
if pl<t then {если она практически 0}
 begin
  write('Это не треугольник, точки лежат на 1 прямой.');
  readln;
  exit;
 end;
for i:=1 to 2 do {рассортируем стороны по возрастаеию}
for j:=i+1 to 3 do
if a[i]>a[j] then
 begin
  x:=a[i];
  a[i]:=a[j];
  a[j]:=x;
 end;
if (abs(a[1]-a[2])>t)and(abs(a[2]-a[3])>t) then{если нет рапвных сторон}
 begin
  write('Это не равнобедренный треугольник.');
  readln;
  exit;
 end;
if sqr(a[1])+sqr(a[2])>=sqr(a[3]) then{если угол против большей стороны
                                       не больше 90 градусов}
 begin
  writeln('Это равнобедренный не тупоугольный треугольник');
  readln;
  exit;
 end;
writeln('Это равнобедренный тупоугольный треугольник');
writeln('Его площадь=',S(x1,y1,x2,y2,x3,y3):0:2);
a[4]:=pl;{добавим в массив площадь и снова отсрортируем}
for i:=1 to 3 do
for j:=i+1 to 4 do
if a[i]>a[j] then
 begin
  x:=a[i];
  a[i]:=a[j];
  a[j]:=x;
 end;
writeln('Его стороны и площадь в порядке возрастания:');
for i:=1 to 4 do
write(a[i]:0:2,' ');
readln
end.