Умножение двух матриц

unit mod_1; //íàçâàíèå ìîäóëÿ
 
interface
const nmax=20; //ìàêñ ðàçìåð ìàòðèö
var n:integer; //ðåàëüíûé ðàçìåð ìàòðèö, îí îäèíàêîâûé, ò.ê. óìíîæàåì
type mtr=array[1..nmax,1..nmax] of integer;//òèï ìàòðèö
procedure vvod(var a:mtr);
procedure vyvod(a:mtr;k:integer;s:string);
procedure umn(a,b:mtr;var c:mtr);
function sled(a:mtr):integer;
procedure obmen(var a:mtr);
 
implementation
procedure vvod(var a:mtr);//ââîä ýëåìåíòîâ ìàòðèö
var i,j:integer;
begin
for i:=1 to n do
for j:=1 to n do
a[i,j]:=random(20);
end;
 
procedure vyvod(a:mtr;k:integer;s:string);//âûâîä ìàòðèö íà ýêðàí,
                                          //k-êîëè÷åñòâî ïîçèöèé ïîä ÷èñëà
var i,j:integer;
begin
writeln(s);
for i:=1 to n do
 begin
  for j:=1 to n do
  write(a[i,j]:k);
  writeln
 end;
end;
 
procedure umn(a,b:mtr;var c:mtr); //óìíîæåíèå 2õ êâàäðàòíûõ ìàòðèö
var i,j,k:integer;
begin
for i:=1 to n do
for j:=1 to n do
 begin
  c[i,j]:=0;
  for k:=1 to n do
  c[i,j]:=c[i,j]+a[i,k]*b[k,j];
 end;
end;
 
function sled(a:mtr):integer; //âû÷èñëåíòå ñëåäà(ñóììà ãëàâíîé äèàãîíàëè)
var i,s:integer;
begin
s:=0;
for i:=1 to n do
s:=s+a[i,i];
sled:=s
end;
 
procedure obmen(var a:mtr); //îáìåí ñòðîê
var i1,i2,j,x:integer;
begin
repeat
writeln('Ââåäèòå íîìåðà äâóõ ñòðîê äëÿ îáìåíà');
read(i1,i2);
until(i1 in [1..n])and(i2 in [1..n])and(i1<>i2);
for j:=1 to n do
 begin
  x:=a[i1,j];
  a[i1,j]:=a[i2,j];
  a[i2,j]:=x
 end
end;
end.

uses crt,mod_1;
var a,b,c:mtr;
begin
randomize;
repeat
write('Ðàçìåð ìàòðèö îò 2 äî ',nmax,' n=');
read(n);
until n in [2..nmax];
vvod(a);
vvod(b);
vyvod(a,3,'Ìàòðèöà A');
writeln('Ñëåä ìàòðèöû=',sled(a));
vyvod(b,3,'Ìàòðèöà B');
writeln('Ñëåä ìàòðèöû=',sled(b));
umn(a,b,c);
vyvod(c,5,'Ìàòðèöà C=A*B');
obmen(c);
vyvod(c,5,'Îáìåí ñòðîê ìàòðèöû Ñ');
end.