Составить процедуру нахождения произведения матрицы А на матрицу B.

program prim_1;
uses crt;
const
n=4;
m=3;
type
matr=array[1..m,1..n] of integer;
matr1=array[1..n,1..m] of integer;
matr2=array[1..n,1..n] of integer;
var
a:matr;
b:matr1;
f:matr2;
i,j,o:integer;
 
procedure proizvedenie(i,j,o:integer;a:matr;b:matr1;f:matr2);
var
r:matr;
z:matr1;
u:matr2;
x,y,c:integer;
begin
randomize;
for x:=1 to m do
for c:=1 to n do
r[x,c]:=random(50)-30;
for x:=1 to m do begin
writeln;
for c:=1 to n do
write(r[x,c]:5);
end;
writeln;
writeln;
 
randomize;
for c:=1 to n do
for y:=1 to m do
z[c,y]:=random(40)-20;
for c:=1 to n do begin
writeln;
for y:=1 to m do
write(z[c,y]:5);
end;
writeln;
writeln;
 
u[x,y]:=0;
for x:=1 to n do
for y:=1 to n do
u[x,y]:=u[x,y]+r[x,c]*z[c,y];
for x:=1 to n do begin
writeln;
for y:=1 to n do
write(u[x,y]:5);
end;
writeln;
end;
 
BEGIN
clrscr;
proizvedenie(i,j,o,a,b,f);
end.

program prim_1;
uses crt;
const nmax=20; {максимальный размер матриц}
type matr=array[1..nmax,1..nmax] of integer;
procedure Vvod(x,y:byte;var m:matr);
var i,j:byte;
begin
for i:=1 to x do
for j:=1 to y do
m[i,j]:=random(5)+1;{взял чтоб не получались очень большие числа и чтоб не было нолей}
end;
 
procedure Vyvod(x,y:byte;var m:matr);
var i,j:byte;
begin
for i:=1 to x do
  begin
    for j:=1 to y do
    write(m[i,j]:4);
    writeln;
  end;
end;
procedure UmnMatr(x,y:byte;m1,m2:matr; var m:matr);
var k,i,j:byte;
begin
for k:=1 to x do
for j:=1 to x do
   begin
     m[k,j]:=0;
     for i:=1 to y do
       begin
         m[k,j]:=m[k,j]+m1[k,i]*m2[i,j];
       end;
   end;
end;
 
var a,b,c:matr;
    n,m:byte;
begin
clrscr;
writeln('Vvedite razmery pryamougolnyh matric m,n:');
repeat
write('m=');readln(m);
write('n=');readln(n);
until(m>0)and(m<=nmax)and(m>0)and(m<=nmax)and(m<>n);
Writeln('Matrica A:');
Vvod(m,n,a);
Vyvod(m,n,a);
Writeln('Matrica B:');
Vvod(n,m,b);
Vyvod(n,m,b);
UmnMatr(m,n,a,b,c);
writeln('Matrica C=A*B:');
Vyvod(m,m,c);
readln
end.