Для квадратных матриц A,B и С одинаковой размерности найти матрицу D

uses crt;
const nmax=20;
type matr=array[1..nmax,1..nmax] of integer;
var m,n:byte;
 
procedure Vvod(var mt:matr);
var i,j:byte;
begin
for i:=1 to m do
for j:=1 to n do
mt[i,j]:=random(20);
end;
 
procedure Vyvod(mt:matr;c:string);
var i,j:byte;
begin
writeln(c,':');
for i:=1 to m do
 begin
  for j:=1 to n do
  write(mt[i,j]:6);
  writeln;
 end;
end;
 
procedure SumMatr(m1,m2:matr;var m:matr);
var i,j:byte;
begin
for i:=1 to n do
for j:=1 to n do
m[i,j]:=m1[i,j]+m2[i,j];
end;
 
procedure RazMatr(m1,m2:matr;var m:matr);
var i,j:byte;
begin
for i:=1 to n do
for j:=1 to n do
m[i,j]:=m1[i,j]-m2[i,j];
end;
 
procedure UmnMatr(m1,m2:matr;var m:matr);
var i,j,k:byte;
begin
for k:=1 to n do
for j:=1 to n do
  begin
   m[k,j]:=0;
   for i:=1 to n do
   m[k,j]:=m[k,j]+m1[k,i]*m2[i,j];
  end;
end;
 
procedure UmnSc(m1:matr;s:integer;var m2:matr);
var i,j:byte;
begin
for i:=1 to n do
for j:=1 to n do
m2[i,j]:=m1[i,j]*s;
end;
var a,a1,a2,a3,b,b1,b2,c,d:matr;
begin
clrscr;
randomize;
repeat
write('Количество строк до ',nmax,' m=');
readln(m);
until m in [1..nmax];
repeat
write('Количество столбцов до ',nmax,' n=');
readln(n);
until n in [1..nmax];
Vvod(a);
Vvod(b);
Vvod(c);
Vyvod(a,'Matrica A');{исходная А}
Vyvod(b,'Matrica B');{исходная B}
Vyvod(c,'Matrica C');{исходная C}
UmnSc(a,3,a1);{3*A}
UmnMatr(b,c,b1);{B*C}
UmnSc(b1,4,b2);{4*B*C}
SumMatr(a1,b2,a3);{3A+4BC}
RazMatr(a3,c,d);{3A+4BC+C}
Vyvod(d,'Матрица 3A+4BC-C');
readln
end.