SUBROUTINE MMULT(A,B,C,I1,I2,I3,N1,N2,N3)
C
C         MATRIX MULTIPLICATION SUBROUTINE
C           A = N1 X N2 MATRIX
C           B = N2 X N3 MATRIX
C           C = A * B
C
      DOUBLE PRECISION A(N1,N2),B(N2,N3),C(N1,N3)
      DO 10 I=1,I1
        DO 10 J=1,I3
          C(I,J)=.0D0
          DO 10 K=1,I2
   10       C(I,J)=C(I,J)+A(I,K)*B(K,J)
      RETURN
      END