Арифметические действия с тремя матрицами - Pascal

uses crt;
const nmax=10;
var n:byte;
type matr=array[1..nmax,1..nmax] of integer;
procedure Vvod(var mt:matr);
var i,j:byte;
begin
for i:=1 to n do
for j:=1 to n do
mt[i,j]:=random(10);
end;
 
procedure Vyvod(mt:matr;c:string);
var i,j:byte;
begin
writeln(c,':');
for i:=1 to n do
  begin
    for j:=1 to n do
    write(mt[i,j]:4);
    writeln;
  end;
write('Press Enter');
readln;
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;
 
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;
 
var w,f,e,f3,ef,ct,r:matr;
begin
clrscr;
randomize;
repeat
write('Размер матриц до ',nmax,' n=');
readln(n);
until n in [1..nmax];
Vvod(w);
Vvod(f);
Vvod(e);
Vyvod(w,'Matrica W');{исходная W}
Vyvod(f,'Matrica F');{исходная F}
Vyvod(e,'Matrica E');{исходная C}
UmnSc(f,3,f3);{3*F}
UmnMatr(e,f3,ef);{E*3F}
SumMatr(w,ef,r);
Vyvod(r,'Matrica R=W+E*3F');
end.