среднее арифметическое диагоналей матрицы? Опции

[bagbir]

Помогите написать процедуру нахождения среднего арифметического для каждой диагонали квадратной матрицы размерности N*N выше и параллельно главной (гл. диагональ невключая)

[bagbir]

program lab2;

{$APPTYPE CONSOLE}

uses
  SysUtils;

const n=4;
type tmatr= Array [1..n, 1..n] of real;

Procedure vvod(var a:tmatr);
  var i,j: 1..n;
Begin
  For i:=1 to n do
    For j:=1 to n do 
    begin
      Write ('A[' , i , ', ' , j , ']=  ');
      Readln (A[i, j]);
    end;
end;
procedure print(const a:tmatr);
  var i,j:integer;
begin
  Writeln('matritca A');
    For i:=1 to n do 
      begin
        For j:=1 to n do
          Write (A[i, j] : 5 : 1);
          Writeln ;
      end;
end;

Procedure SRARIF(const A:tmatr);
 var i,j: 1..n;
      s:real;
begin
  for i:=2 to n do
    begin
      s:=0;
      for j:=i to n do
        s:=s+A[j-i+1,j];
        WriteLn(i,'  ',s/(n-i+1):5:1);
    end;
    Readln
end;
var
A:tmatr;
begin
  vvod(A);
  print(A);
  SRARIF(A);
end.

Появился еще один вопрос(( надо чтобы ответ выводился через массив в отдельной процедуре!