Перемножить три матрицы

unit matrix;
interface
const kolvostrok=2;
type matr=array [1..kolvostrok,1..kolvostrok];
procedure ololo(mA,mB,mC:matr; var mD,mE,mG:matr);
var i,j,k:integer;
implementation
procedure ololo(mA,mB,mC:matr; var mD,mE,mG:matr);
begin
{----------------------------------------}
for i:=1 to kolvostrok do
for j:=1 to kolvostrok do begin
mD[i,j]:=0;
for k:=1 to kolvostrok do
mD[i,j]:=mD[i,j]+mA[i,k]*mB[k,j];
end;
{--------------------------------------}
for i:=1 to kolvostrok do
for j:=1 to kolvostrok do
mE[i,j]:=0;
mE[i,j]:=mC[j,i];
{---------------------------------------}
for i:=1 to kolvostrok do
for j:=1 to kolvostrok do begin
mG[i,j]:=0;
for k:=1 to kolvostrok do
mG[i,j]:=mG[i,j]+mD[i,k]*mE[k,j];
end;
end.

uses matrix;
const m=2;
var A,B,C,D,E,G:matr;
    mas:real;
    i,j,k:integer;
begin
for i:=1 to m do
for j:=1 to m do begin
writeln('Vvedite i,j element massiva A');
readln(mas);
A[i,j]:=mas;
end;
{--------------------------------}
for i:=1 to m do
for j:=1 to m do begin
writeln('Vvedite i,j element massiva B');
readln(mas);
B[i,j]:=mas;
end;
{---------------------------------}
for i:=1 to m do
for j:=1 to m do begin
writeln('Vvedite i,j element massiva C');
readln(mas);
C[i,j]:=mas;
end;
{-----------------------------------}
writeln('Rezultst raboty modulya');
matrica(A,B,C,D,E,G);
for i:=1 to m do
for j:=1 to m do begin
writeln(G[i,j]);
end;
end.

unit matrix;
interface
const max=20;
type Tmatr=array[1..max,1..max] of integer;
procedure Vvod(var mt:Tmatr;var n:byte);
procedure Vyvod(var m:Tmatr;n:byte;c:string);
procedure Trans(m1:Tmatr;n:byte; var m2:Tmatr);
procedure UmnMatr(m1,m2:Tmatr;n:byte;var m:Tmatr);
implementation
procedure Vvod(var mt:Tmatr;var n:byte);
var i,j:byte;
begin
for i:=1 to n do
for j:=1 to n do
mt[i,j]:=random(9)+1;
end;
 
procedure Vyvod(var m:Tmatr;n:byte;c:string);
var i,j:byte;
begin
writeln(c,':');
for i:=1 to n do
  begin
    for j:=1 to n do
    write(m[i,j]:5);
    writeln;
  end;
writeln;
end;
procedure Trans(m1:Tmatr;n:byte; var m2:Tmatr);
var i,j:byte;
begin
for i:=1 to n do
for j:=1 to n do
m2[i,j]:=m1[j,i];
end;
 
procedure UmnMatr(m1,m2:Tmatr;n:byte;var m:Tmatr);
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;
end.

uses crt,matrix;
var a,b,c,ct,m1,m2:Tmatr;
    n:byte;
begin
clrscr;
randomize;
repeat
write('Razmer matric do ',max,'=');
readln(n);
until n in [1..max];
Vvod(a,n);
Vvod(b,n);
Vvod(c,n);
Vyvod(a,n,'Matrica A');
Vyvod(b,n,'Matrica B');
Vyvod(c,n,'Matrica C');
Trans(c,n,ct);
Vyvod(ct,n,'Transponirovannaya matrica Ct');
write('Press Enter...');
readln;
UmnMatr(a,b,n,m1);
Vyvod(m1,n,'Matrica M1=A*B');
UmnMatr(m1,ct,n,m2);
Vyvod(m2,n,'Matrica M2=A*B*Ct');
write('Press Enter...');
readln
end.