Умножение квадратных матриц делается так:
const n=5;
type ar=array[1..n,1..n] of integer;
procedure p3(var m1,m2,m3:ar);
var i,j,k,s:integer;
begin
for i:=1 to n do
for j:=1 to n do
begin
s:=0;
For k := 1 to n do s:=s+m1[i,k]*m2[k,j];
m3[i,j]:=s;
end;
end;