brmul.for

来自「Fortran基本算法源程序,线性方程,矩阵运算,非线性方程,极值问题」· FOR 代码 · 共 14 行

	SUBROUTINE BRMUL(A,B,M,N,K,C)
!       实数矩阵乘法
	DIMENSION A(M,N),B(N,K),C(M,K)
	DOUBLE PRECISION A,B,C
	DO 50 I=1,M
	DO 50 J=1,K
	  C(I,J)=0.0
	  DO 10 L=1,N
	    C(I,J)=C(I,J)+A(I,L)*B(L,J)
10	  CONTINUE
50	CONTINUE
	RETURN
	END