Найти сумму элементов главной диагонали матрицы А:- Pascal(Паскаль)

1 2 3 0,9
6 4 1 1,2
5 -1,1 5 -3,6
5 -1,1 5 -3,6

Uses
 CRT;
Const
 a:array[1..4,1..4] of real=((1, 2, 3, 0.9),
                             (6, 4, 1, 1.2),
                             (5, -1.1, 5, -3.6),
                             (5, -1.1, 5, -3.6));
Var
 sum:real;
 i,j:integer;
begin
sum:=0;
for i:=1 to 4 do
 begin
  for j:=1 to 4 do
    write(a[i,j]:6:1);
   sum:=sum+a[i,i];
  writeln;
 end;
write('SUM:=',sum);
end.