code 34013;
   real procedure MATMAT(L, U, I, J, A, B); value L,U,I,J;
   integer L,U,I,J; array A,B;
   begin integer K; real S;
       S:= 0;
       for K:=L step 1 until U do S:= A[I,K] * B[K,J] + S;
       MATMAT:= S
   end MATMAT

       eop