Матрицы. Найти произведение ненулевых элементов, расположенных ниже главной диагонали

uses crt;
const n=7;
type matr = array [1..n,1..n] of integer;
var i,j,k,p: integer;
    a,b,c: matr;
 
procedure prod(m: matr; var p: integer);
begin
 p:=1;
 for i:=1 to n do
  for j:=1 to n do
   if i>j then
    if m[i,j]<>0 then
     p:=p*m[i,j];
end;
 
begin
 clrscr;
 randomize;
 for i:=1 to n do
  begin
   for j:=1 to n do
    begin
     a[i,j]:=random(10)-4;
     write(a[i,j]:3);
    end;
   writeln;
  end;
 writeln;
 prod(a,p);
 writeln('Произведение:',p);
 for i:=1 to n do
  begin
   for j:=1 to n do
    begin
     b[i,j]:=random(10);
     write(b[i,j]:3);
    end;
   writeln;
  end;
 writeln;
 prod(b,p);
 writeln('Произведение:',p);
 for i:=1 to n do
  begin
   for j:=1 to n do
    begin
     c[i,j]:=random(10)-4;
     write(c[i,j]:3);
    end;
   writeln;
  end;
 writeln;
 prod(c,p);
 writeln('Произведение:',p);
 readln;
end.