Транспонирование и умножение прямоугольных матриц

uses crt;
const max=20;
type matr=array[1..max,1..max] of integer;
{создание матриц}
procedure Vvod(var mt:matr;x,y:byte);
var i,j:byte;
begin
for i:=1 to x do
for j:=1 to y do
mt[i,j]:=random(9)+1;
end;
{вывод матриц на экран}
procedure Vyvod(var mt:matr;x,y:byte;c:string);
var i,j:byte;
begin
writeln('Матрица ',c);
for i:=1 to x do
 begin
  for j:=1 to y do
  write(mt[i,j]:4);
  writeln;
 end;
end;
{траеспонирование матриц}
procedure Trans(m1:matr;x,y:byte;var m2:matr);
var i,j:byte;
begin
for i:=1 to y do
for j:=1 to x do
m2[i,j]:=m1[j,i];
end;
{умножение исходной матрицы на транспонированную}
procedure Umn(m1,m2:matr;x,y:byte;var m3:matr);
var i,j,k:byte;
begin
for k:=1 to x do
for j:=1 to x do
 begin
  m3[k,j]:=0;
  for i:=1 to y do
  m3[k,j]:=m3[k,j]+m1[k,i]*m2[i,j];
 end;
end;
 
var a,b,at,bt,c:matr;
    n,m:byte;
begin
clrscr;
repeat
write('Количество строк от 1 до ',max,' n=');readln(n);
until n in [1..max];
repeat
write('Количество столбцов до ',max,' m=');readln(m);
until m in [1..max];
Vvod(a,n,m);
Vvod(b,n,m);
Vyvod(a,n,m,'A');
Trans(a,n,m,at);
Vyvod(at,m,n,'At');
Umn(a,at,n,m,c);
Vyvod(c,n,n,'A*At');
write('Press Enter...');
readln;
Vyvod(b,n,m,'B');
Trans(b,n,m,bt);
Vyvod(bt,m,n,'Bt');
Umn(b,bt,n,m,c);
Vyvod(c,n,n,'B*Bt');
write('Press Enter...');
readln
end.

for i:=1 to n do
for j:=1 to m do
c[i,j]:=a[i,j]+b[i,j];