Найти площадь пирамиды, требуется простенькое решение Опции

[Шушпанчик]

Вообщем задача сформулирована следующим образом:

"Написать программу нахождения площади поверхности пирамиды по координатам ее вершин. Подготовить и использовать процедуру определения площади треугольника по известным вершинам".

Задача простенькая с математической и логической т.з., но вот как её описывать в паскале - напрочь забыто.
Желательно использовать базовый инструментарий, такой как: линейные, разветвляющиеся и циклические алгоритмы, массивы, процедуры и функции.

Заранее большое человеческое спасибо! :p2:

[Guest]

Вообщем вот полный код.
Проблема в одном : после процедуры в значениях выходят нули.

var 
  x: Array [1..5] of Real;
  y: Array [1..5] of Real;
  z: Array [1..5] of Real;

{ Площадь треугольника вычисляем по формуле Герона,
 но сам треугольник задаем не отдельными координатами,
 а вершинами }

Function Plosh(x1, y1, z1, x2, y2, z2, x3, y3, z3: TPoint): Real;
 Var
   a, b, c, p: Real;
 Begin
   a:=sqrt(sqr(y2-y3)+sqr(x2-x3)+sqr(z2-z3));
   b:=sqrt(sqr(y1-y3)+sqr(x1-x3)+sqr(z1-z3));
   c:=sqrt(sqr(y1-y2)+sqr(x1-x2)+sqr(z1-z2));
   p:=(a + b + c)/2;
   Plosh := p * Sqrt((p-a)*(p-B)*(p-c))
 End;

Var
 s: Real;
 i, j: Integer;
 pyramide: TPyramide;

Begin
 { заносим в массив координаты вершин основания пирамиды }
 For i := 1 To 4 Do
     Begin
      if i:=1 then  
       WriteLn( 'Вершина основания ', i );
       Write( 'X = ' ); ReadLn(X[i]);
       Write( 'Y = ' ); ReadLn(Y[i]);
       Write( 'Z = ' ); ReadLn(Z[i]);
     else
       WriteLn( 'Вершина основания ', i );
       Write( 'X = ' ); ReadLn(X[i]);
       Write( 'Y = ' ); ReadLn(Y[i]);
       z[i]:=z[1];
     End;



 WriteLn( 'Верхняя точка пирамиды:' );
   Begin
     Write( 'X = ' ); ReadLn(X[5]);
     Write( 'Y = ' ); ReadLn(Y[5]);
     Write( 'Z = ' ); ReadLn(Z[5]);
   End;

 s := 0;
 For i := 1 To 4 Do
     Begin
       { считаем сумму площадей боковых граней }
       If i <> 4 Then j := Succ(i) Else j := 1;
       s := s + Plosh(x[i], y[i], z[i], x[j], y[j], z[j], x[5], y[5], z[5]);
     End;

 s := s + Plosh(x[1], y[1], z[1], x[2], y[2], z[2], x[3], y[3], z[3]);
 s := s + Plosh(x[2], y[2], z[2], x[3], y[3], z[3], x[4], y[4], z[4]);

 WriteLn( 'Площадь пирамиды s = ', s )

End.