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

function length(x1,y1,x2,y2: real): real;
begin
  length:=sqrt(sqr(x2-x1)+sqr(y2-y1));
end;
 
function area(x1,y1,x2,y2,x3,y3: real): real;
var a,b,c,p: real;
begin
  a:=length(x1,y1,x2,y2);
  b:=length(x2,y2,x3,y3);
  c:=length(x1,y1,x3,y3);
  p:=(a+b+c)/2;
  area:=sqrt(p*(p-a)*(p-b)*(p-c));
end;
 
var
  x1,y1,x2,y2,x3,y3: real;
  s1,s2: real;
begin
  writeln('Âû÷èñëåíèå ïëîùàäè 1-ãî òðåóãîëüíèêà');
  write('Ââåäèòå êîîðäèíàòû 1-é âåðøèíû: '); readln(x1,y1);
  write('Ââåäèòå êîîðäèíàòû 2-é âåðøèíû: '); readln(x2,y2);
  write('Ââåäèòå êîîðäèíàòû 3-é âåðøèíû: '); readln(x3,y3);
  s1:=area(x1,y1,x2,y2,x3,y3);
  writeln('Âû÷èñëåíèå ïëîùàäè 2-ãî òðåóãîëüíèêà');
  write('Ââåäèòå êîîðäèíàòû 1-é âåðøèíû: '); readln(x1,y1);
  write('Ââåäèòå êîîðäèíàòû 2-é âåðøèíû: '); readln(x2,y2);
  write('Ââåäèòå êîîðäèíàòû 3-é âåðøèíû: '); readln(x3,y3);
  s2:=area(x1,y1,x2,y2,x3,y3);
  writeln('Ïëîùàäü 1-ãî òðåóãîëüíèêà = ',s1:0:3);
  writeln('Ïëîùàäü 2-ãî òðåóãîëüíèêà = ',s2:0:3);
  if s1>s2 then writeln('1-é òðåóãîëüíèê áîëüøå')
    else if s2>s1 then writeln('2-é òðåóãîëüíèê áîëüøå')
      else writeln('Òðåóãîëüíèêè ðàâíû');
end.