Вычислить и вывести в квадратной матрице, порядка N произведение элементов, находящихся выше главной диагонали — Pascal(Паскаль)

uses crt;
const nmax=20;
var a:array[1..nmax,1..nmax] of integer;
    n,i,j:byte;
    p:longint;
begin
clrscr;
randomize;
repeat
write('Размер матрицы=');
readln(n);
until n in [1..nmax];
writeln('Исходная матрица:');
for i:=1 to n do
 begin
   for j:=1 to n do
    begin
     a[i,j]:=random(9)+1;
     write(a[i,j]:4);
    end;
   writeln;
 end;
writeln;
p:=1;
for i:=1 to n do
for j:=i+1 to n do
p:=p*a[i,j];
write('P=',p);
readln
end.