FORTRAN两实数矩阵相乘程序.docxVIP

!两实数矩阵相乘程序 PROGRAM Matrix_multiply !定义动态数组A,B,S REAL,ALLOCATABLE :: A(:,:) REAL,ALLOCATABLE :: B(:,:) REAL,ALLOCATABLE :: S(:,:) !矩阵维数定义 INTEGER M,N,L,I,J,K PRINT *,请输入矩阵A的大小,格式M N READ *,M,N PRINT *,请输入矩阵B的大小,格式N L READ *,N,L ALLOCATE(A(M,N)) ALLOCATE(B(N,L)) ALLOCATE(S(M,L)) !输入A,B矩阵 PRINT *,请输入矩阵A,M,行,N,列 DO I=1,M PRINT *,请输入第,I,行 READ *,A(I,:) ENDDO PRINT *,请输入矩阵B,N,行,L,列 DO I=1,N PRINT *,请输入第,I,行 READ *,B(I,:) ENDDO !计算A,B相乘 DO I=1,M DO J=1,L S(I,J)=0 DO K=1,N S(I,J)=S(I,J)+A(I,K)*B(K,J) ENDDO ENDDO ENDDO !打印输出 DO I=1,M PRINT *,S(I,:) ENDDO END