Определить суммы всех элементов матрицы, лежащих на диагоналях Pascal(Паскаль)

const
  n = 3;
 
var
  j, i: integer;
  s: real;
  a: array [1..n, 1..n] of real;
 
begin
  for i := 1 to n do
  begin
    for j := 1 to n do
    begin
      a[i, j] := random * n;
      write(a[i, j]:8:2);
    end;
    writeln;
  end;
  s := 0;
  for i := 1 to n do
  begin
    s := s + a[i, i];
    s := s + a[i, n - i + 1];
  end;
  if Odd(n) then 
    s := s - a[n div 2 + 1][n div 2 + 1];
  writeln('Sum: ', s);  
end.