juzhenchengfa.f90

来自「计算卫星位置的fortran程序」· F90 代码 · 共 21 行

  	SUBROUTINE MATMUL(A,B,BR,BC,C,CR,CC)  !矩阵乘积
	IMPLICIT NONE
	INTEGER BR,BC,CR,CC
	REAL*8 B(BR,BC),C(CR,CC),A(BR,CC)
	INTEGER I,J,K
	!BC若不等于CR，这两个矩阵无法相乘
	IF(BC.NE.CR)THEN
	WRITE(*,*)"Matrix size errir"
	stop
	end if

	DO I=1,BR
	  DO J=1,CC
	     A(I,J)=0
	     DO K=1,CR
	     A(I,J)=A(I,J)+B(I,K)*C(K,J)
	     END DO 
	  END DO
	END DO
	RETURN
	END