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

program PR;
uses crt;
type
  tt = array [1..10,1..10] of real;
 
procedure SP (A: tt; N:byte);
var i:byte;
    R : real;
begin
  i := 1;
  R := 0;
  repeat
   R := R+A[i,i];
   i := i+1;
  until i>N;
writeln('След=',R:0:2);
end;
var  B,C: tt;
     n,i,j:byte;
begin
randomize;
clrscr;
repeat
Write('Размер матрицы В до 10 n=');
readln(n);
until n in [1..10];
writeln('Матрица:');
for i := 1 to n do
 begin
   for j := 1 to n do
    begin
     b[i,j]:=10*random;
     write(b[i,j]:5:2);
    end;
   writeln;
  end;
writeln;
SP(b,n);
write('Нажмите Enter');
readln;
repeat
Write('Размер матрицы C до 10 n=');
readln(n);
until n in [1..10];
writeln('Матрица:');
for i := 1 to n do
 begin
   for j := 1 to n do
    begin
     c[i,j]:=10*random;
     write(c[i,j]:5:2);
    end;
   writeln;
  end;
writeln;
SP(c,n);
readln
END.