PROCEDURE MMULT (A,B:MATRIX ; VAR C:MATRIX ; N:INTEGER) ;
    VAR I,J,K:INTEGER ;
    BEGIN
      FOR I:=1 TO N DO
      BEGIN
        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
    END ;