Функции: вычисление длины отрезка по координатам и площади треугольника по трем сторонам

var
  x1, y1, x2, y2, x3, y3: real;
 
function len(x1, y1, x2, y2: real): real;
begin
  len := sqrt(sqr(x2 - x1) + sqr(y2 - y1));
end;
 
function S(a, b, c: real): real;
var
  p: real;
begin
  p := (a + b + c) / 2;
  S := sqrt(p * (p - a) * (p - b) * (p - c));
end;
 
begin
  writeln('Введите x1, y1, x2, y2, x3, y3');
  readln(x1, y1, x2, y2, x3, y3);
  writeln(S(len(x1, y1, x2, y2), len(x2, y2, x3, y3), len(x1, y1, x3, y3)));
end.

function Len(x1, y1, x2, y2 : Real) := Sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2));
function Area0(a,b,c,d,e,f : Real) := Sqrt((Len(a,b,c,d)+Len(c,d,e,f)+Len(e,f,a,b))*(-Len(a,b,c,d)+Len(c,d,e,f)+Len(e,f,a,b))*(Len(a,b,c,d)-Len(c,d,e,f)+Len(e,f,a,b))*(Len(a,b,c,d)+Len(c,d,e,f)-Len(e,f,a,b))) / 4; // Герон
function Area1(a,b,c,d,e,f : Real) := Abs(a*d+c*f+b*e-b*c-d*e-a*f) / 2; // универсальный способ нахождения площади многоугольника
function Area2(a,b,c,d,e,f : Real) := Abs(a*(d-f)+c*(f-b)+e*(b-d)) / 2; // через определитель матрицы способ 1
function Area3(a,b,c,d,e,f : Real) := Abs((c-a)*(f-b)-(e-a)*(b-d)) / 2; // через определитель матрицы способ 2
begin
  var a,b,c,d,e,f : Real;
  (a,b,c,d,e,f) := (ReadReal('Введите через пробел координаты вершин треугольника :'),ReadReal,ReadReal,ReadReal,ReadReal,ReadLnReal);
  WriteLn('Площадь треугольника Героном   = ', Area0(a,b,c,d,e,f));
  WriteLn('Площадь треугольника 1 формула = ', Area1(a,b,c,d,e,f));
  WriteLn('Площадь треугольника 2 формула = ', Area2(a,b,c,d,e,f));
  WriteLn('Площадь треугольника 3 формула = ', Area3(a,b,c,d,e,f));
end.