Вычислить произведение чисел, принадлежащих промежутку A, B и лежащих на главной диагонали и выше ее

...
mmul:=1;
For i:=1 to N do
  For j:=i to N do
    If (M[i,j]>=A) and (M[i,j]<=B) then mmul:=mmul*M[i,j];
...

program work8;
var
  a: array[1..20,1..20] of real;
  ar,p,n,m,i,j,s: real;
begin
  randomize;
  repeat
    write('Введите размеры матрицы n, m через пробел: ');
    readln(n,m);
  until (n in [1..20]) and (m in [1..20]);
  writeln('Матрица:');
  begin
 p:=1;
for i:=1 to n do
  for j:=i to n do
    if (ar[i,j]>=a) and (ar[i,j]<=b) then p:=p*ar[i,j];
    end;
    end.

program work8;
var
  ar: array[1..20,1..20] of real;
  p,a,b: real;
  i,j,n,m:integer;
begin
 
  randomize;
  repeat
    write('Введите размеры матрицы n, m через пробел: ');
    readln(n,m);
  until (n in [1..20]) and (m in [1..20]);
  write('Введите A: ');
  readln(A);
  write('Введите B: ');
  readln(B);
 
  for i:=1 to n do
    begin
      for j:=1 to n do
        begin
          ar[i,j]:=random(n);
          write(ar[i,j]:8:2);
        end;
      writeln;
    end;
 
  writeln('Результат:');
  p:=1;
  for i:=1 to n do
    for j:=i to n do
      if (ar[i,j]>=a) and (ar[i,j]<=b) then
        p:=p*ar[i,j];
  writeln(p:8:2);
  readln;
end.