06/06/29 21:02:47
>>557の続きだお
SUBROUTINE MTINP(A, M, N)
DIMENSION A(M,N)
READ(*,*) ((A(I,J), J=1,N), I=1,M)
END
C
SUBROUTINE MTPRT(A, M, N)
DIMENSION A(M,N)
CHARACTER F*20
WRITE(F,'(''(1X,'',I3,''F7.1)'')') N
WRITE(*,F) ((A(I,J), J=1,N), I=1,M)
END
C
SUBROUTINE MTMUL(A, B, C, L, N, M)
DIMENSION A(L,N), B(N,M), C(L,M)
DO 10 J=1, M
DO 20 I=1, L
C(I,J) = 0.
DO 30 K=1, N
C(I,J) = C(I,J) + A(I,K) * B(K,J)
30 CONTINUE
20 CONTINUE
10 CONTINUE
END