Определение площади шестиугольника по заданным координатам вершин

type point = record
               x, y:real;
             end;
 
function l(a, b:point):real;
begin
  l:= sqrt(sqr(a.x - b.x) + sqr(a.y - b.y));
end;
 
function s(x, y, z:point):real;
var p, a, b, c:real;
begin
  a:= l(x, y);
  b:= l(x, z);
  c:= l(y, z);
  p:= (a + b + c) / 2;
  s:= sqrt(p * (p - a) * (p - b) * (p - c));
end;
 
function solution(a, b, c, d, e, f:point):real;
begin
  solution:= s(a, b, c) + s(a, c, d) + s(a, d, f) + s(d, e, f);
end;
 
var a, b, c, d, e, f:point;
 
begin
  writeln('Введите координаты точек шестиугольника');
  readln(a.x, a.y);
  readln(b.x, b.y);
  readln(c.x, c.y);
  readln(d.x, d.y);
  readln(e.x, e.y);
  readln(f.x, f.y);
  writeln('Площадь шестиугольника равна: ', solution(a, b, c, d, e, f));
end.

type point = record
               x, y:real;
             end;
 
procedure l(a, b:point; var r:real);
begin
  r:= sqrt(sqr(a.x - b.x) + sqr(a.y - b.y));
end;
 
procedure s(x, y, z:point; var t:real);
var p, a, b, c:real;
begin
  l(x, y, a);
  l(x, z, b);
  l(y, z, c);
  p:= (a + b + c) / 2;
  t:= sqrt(p * (p - a) * (p - b) * (p - c));
end;
 
procedure solution(a, b, c, d, e, f:point; var p:real);
var s1, s2, s3, s4: real;
begin
  s(a, b, c, s1);
  s(a, c, d, s2);
  s(a, d, f, s3);
  s(d, e, f, s4);
  p:= s1 + s2 + s3 + s4;
end;
 
var a, b, c, d, e, f:point;
    pl:real;
 
begin
  writeln('Введите координаты точек шестиугольника');
  readln(a.x, a.y);
  readln(b.x, b.y);
  readln(c.x, c.y);
  readln(d.x, d.y);
  readln(e.x, e.y);
  readln(f.x, f.y);
  solution(a, b, c, d, e, f, pl);
  writeln('Площадь шестиугольника равна: ', pl);
end.