Написать программу нахождения суммы диагональных элементов A*B.

Sum:=0;
for i:=1 to n do
  Sum:=Sum+C[i,i];

uses crt;
const nmax=8;
var A,B,C:array [1..nmax, 1..nmax] of real;
    n,i,j,k:byte;
    s:real;
begin
clrscr;
repeat
write ('введите размер матриц до ',nmax,' n=');
readln (n);
until n in [1..nmax];
writeln ('введите элементы матрицы А: ');
for i:=1 to n do
for j:=1 to n do
begin
write ('A[',i,',',j,']=');
readln(A[i,j]);
end;
writeln ('введите элементы матрицы B : ');
for i:=1 to n do
for j:=1 to n do
begin
write ('B[',i,',',j,']=');
readln(B[i,j]);
end;
clrscr;
writeln('Matrica A:');
for i:=1 to n do
  begin
    for j:=1 to n do
    write(a[i,j]:6:2);
    writeln;
  end;
writeln;
writeln('Matrica B:');
for i:=1 to n do
  begin
    for j:=1 to n do
    write(b[i,j]:6:2);
    writeln;
  end;
writeln;
{умножение матриц}
for k:=1 to n do
for j:=1 to n do
  begin
    C[k,j]:=0;
    for i:=1 to n do
    C[k,j]:=C[k,j]+A[k,i]*B[i,j];
  end;
writeln ('матрица С=A*B:');
s:=0;
for i:=1 to n do
   begin
     for j:=1 to n do
      begin
       write (C[i,j]:8:2);
       if i=j then s:=s+C[i,j];{сумма элементов главной диагонали}
      end;
     writeln;
   end;
writeln;
write('Сумма диагональных элементов=',s:0:2);
readln
end.