Найти произведение элементов двумерного массива находящихся на главной диагонали матрицы размером (10×10) (с программой и алгоритмом)

//Элементы a[i,i] - элементы главной диагонали, в переменной p будем хранить их произведение.

var  B: array[1..100, 1..100] of integer;
i, j, n: integer;  P: longint;
begin
write('Введите размерность: ');
readln(n);
p := 1;
for i := 1 to n do  begin
for j := 1 to n do  begin 
     B[i, j] := Random(100); 
     write(B[i, j]:3);
     if i = j then p := p * B[i, j];
end; 
writeln; 
end;
writeln('Произведение главной диагонали матрицы: ', p);
end.