type ma3x=array[1..50,1..50] of integer;
............
procedure pr_ma3x(m1,t,n2:integer;a1,a2:ma3x;var pr:ma3x);
{n1=m2=t}
var i,j,k:integer;
begin
 for i:=1 to m1 do begin
   for j:=1 to n2 do begin
     for k:=1 to t do pr[i,j]:=pr[i,j]+a1[i,k]*a2[k,j]
   end
 end
end;

Const
 sizeX = 5;
 sizeY = 3;

Type
 TOneMatrix =
   Array[1 .. sizeX, 1 .. sizeY] Of Real;
 TTwoMatrix =
   Array[1 .. sizeY, 1 .. sizeX] Of Real;
 TResult =
   Array[1 .. sizeX, 1 .. sizeX] Of Real;

(* matrix multiplication *)
Procedure matrixMult(Var m: TResult;
         a: TOneMatrix; b: TTwoMatrix);
 Var
   i, j, k: Integer;
 Begin
   For i := 1 To sizeX Do
     For j := 1 To sizeX Do
       Begin
         m[i, j] := 0;
         For k := 1 To sizeY Do
           m[i, j] := m[i, j] + a[i, k] * b[k, j]
       End;
 End;

Var
 m1: TOneMatrix;
 m2: TTwoMatrix;
 r: TResult;
Begin
 {заполняешь матрицу m1(5x3)}
 {заполняешь матрицу m2(3x5)}
 matrixMult(r, m1, m2);
 {получаешь матрицу r(5x5)}

End.