В заданной прямоугольной таблице найти сумму диагональных элементов - Free Pascal

var
  A: array of array of longint;
  F: text;
  N, ans, k: longint;
 
procedure InputData;
begin//Ввод матрицы из файла input.txt . 1ая строка файла - размеры матрицы (N).
  Assign(F, 'input.txt');
  Reset(F);
  Readln(F, N);
  SetLength(A, N + 1, N + 1);
  for N := 1 to N do for k := 1 to N do Read(F, A[N, k]);
  close(F);
end;
 
procedure ProcessData;
begin
  for K := 1 to N do ans := ans + A[k, k];
  for k := N downto 1 do ans := ans + A[k, k];
end;
 
procedure OutPutData;
begin
  Write('Sum = ', Ans);
end;
 
begin
  ans := 0;
  InputData;
  ProcessData;
  OutPutData;
end.