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

uses crt;

var
  a: array [1 .. 10, 1 .. 10] of integer;
  i, j: integer;
  p: longint;

begin
  clrscr;
  randomize;
  p := 1;
  for i := 1 to 5 do
  begin
    for j := 1 to 5 do
    begin
      a[i, j] := random(7) + 1;
      write(a[i, j]:3);
    end;
    writeln;
  end;
  writeln;
  for i := 1 to 5 do
  begin
    for j := 1 to 5 do
    begin
      if i <= j then
        p := p * a[i, j];
    end;
  end;
  writeln(p);
  readln;

end.