用Fortran语言编写两个4*4矩阵乘的程序

IMPLICT NONE
REAL::A(4,4),B(4,4),C(4,4)=0
INTEGER::I,J,K
READ(*,*) A,B !输入方式也可改为读文件
DO I= 1,4
DO J=1,4
DO K=1,4
C(I,J)=C(I,J)+A(I,K)*B(K,J)
ENDDO
ENDDO
ENDDO
WRITE(*,*) C
END