Составить подпрограмму–функцию для нахождения суммы элементов матрицы, расположенных выше главной диагонали

uses crt;
const nmax=5;
type matr=array[1..nmax,1..nmax] of real;
var n:byte;{общий размер матриц}
{процедура ввода исходных матриц}
procedure Vvod(var a:matr;n:byte;c:char);
var i,j:byte;
begin
writeln('Матрица ',c);
for i:=1 to n do
 begin
  for j:=1 to n do
   begin
    a[i,j]:=-2+5*random;
    write(a[i,j]:7:2)
   end;
  writeln
 end;
end;
{процедура умножения квадратных матриц}
procedure Umn(a1,a2:matr;var a3:matr;n:byte);
var k,j,i:byte;
begin
for k:=1 to n do
for j:=1 to n do
 begin
  a3[k,j]:=0;
  for i:=1 to n do
  a3[k,j]:=a3[k,j]+a1[k,i]*a2[i,j];
 end;
end;
{функция по условию}
function Summa(a:matr;n:byte):real;
var i,j:byte;
    s:real;
begin
s:=0;
for i:=1 to n-1 do
for j:=i+1 to n do
s:=s+a[i,i];
Summa:=s;
end;
 
var a,b,a2,as,y1,b2,b3,bs,y2,c,x,y:matr;
    i,j,t:byte;
    mx:real;
    s:string;
begin
clrscr;
randomize;
repeat
write('Размер матриц до ',nmax,' n=');
readln(n);
until n in [1..nmax];
Vvod(a,n,'A');
Vvod(b,n,'B');
write('Нажмите Enter для продолжения');
readln;
Umn(a,a,a2,n);{A^2}
Umn(b,b,b2,n);
Umn(b2,b,b3,n);{B^3}
writeln('Матрица X:');
for i:=1 to n do
 begin
  for j:=1 to n do
   begin
    x[i,j]:=a2[i,j]+b3[i,j];
    write(x[i,j]:12:2)
   end;
  writeln
 end;
writeln('Сумма элементов выше главной диагонали=',Summa(x,n):0:2);
repeat
write('Введите натуральное число от 2 до 7 - степень матриц T=');
readln(t);
until t in[2..7];
as:=a;{A^t}
for i:=1 to t-1 do
 begin
  Umn(as,a,c,n);
  as:=c;
 end;
Umn(a,as,y1,n);{A*A^t}
bs:=b;{B^t}
for i:=1 to t-1 do
 begin
  Umn(bs,b,c,n);
  bs:=c;
 end;
Umn(b,bs,y2,n);{B*B^t}
writeln('Матрица Y:');
for i:=1 to n do
 begin
  for j:=1 to n do
   begin
    y[i,j]:=y1[i,j]+y2[i,j];
    write(y[i,j]:12:2)
   end;
  writeln
 end;
writeln('Сумма элементов выше главной диагонали=',Summa(y,n):0:2);
if Summa(x,n)>Summa(y,n) then mx:=Summa(x,n)else mx:=Summa(y,n);
write('Максимальная сумма=',mx:0:2);
readln
end.