Вычислить площади треугольников и выбрать наибольшую из площадей

Var p,S,Smax : Real;
     i : Byte;
...
Smax:=0;
For i:=1 to N do
Begin
  p:=(a[i]+b[i]+c[i])/2;
  S:=sqrt(p*(p-a[i])*(p-b[i])*(p-c[i]));
  Writeln(i,' S= ',S)
  If S>Smax then Smax:=S;
End;
Writeln('максимальная площадь: ',Smax);
...

uses crt;
function trg(a,b,c:real):boolean;
begin
trg:=(a<b+c)and(b<a+c)and(c<a+b)
end;
function pls(a,b,c:real):real;
var p:real;
begin
p:=(a+b+c)/2;
pls:=sqrt(p*(p-a)*(p-b)*(p-c))
end;
const nmax=10;
var a,b,c:array[1..nmax] of real;
    n,i,imx:integer;
    p,mx:real;
begin
clrscr;
randomize;
repeat
write('Количество треугольников от 1 до ',nmax,' n=');
readln(n);
until n in [1..nmax];
for i:=1 to n do
 begin
  a[i]:=1+20*random;
  b[i]:=1+20*random;
  c[i]:=1+20*random;
 end;
write('  ');
for i:=1 to n do
write(i:6);
writeln;
write('a:');
for i:=1 to n do
write(a[i]:6:1);
writeln;
write('b:');
for i:=1 to n do
write(b[i]:6:1);
writeln;
write('c:');
for i:=1 to n do
write(c[i]:6:1);
writeln;
write('s:');
mx:=0;
imx:=0;
for i:=1 to n do
 begin
  if trg(a[i],b[i],c[i]) then p:=pls(a[i],b[i],c[i])
  else p:=0;
  write(p:6:1);
  if p>mx then
   begin
    mx:=p;
    imx:=i
   end;
 end;
writeln;
write('Максимальная площадь у треугольника ',imx,' = ',mx:0:1);
readln
end.