Matrix multiplication by subroutine subprogram

Find the multiplication of any of the matrix with (3 by 3) dimension by applying subroutine subprogram

Solution:

       SUBROUTINE AL(A,B,C,M,N,L)

       DIMENSION A(M,N),B(M,N),C(M,N)

       DO I=1,M

       DO J=1,N

       C(I,J)=0.0

       Do k = 1,L

       C(I,J)=C(I,J)+A(I,L)*B(L,J)

       ENDDO

       ENDDO

       ENDDO

       RETURN

       END

       DIMENSION A(3,3),B(3,3),C(3,3)

       PRINT*, 'ENTER THE VALUE OF MATRIX A'

       DO I=1,3

       DO J=1,3

       READ*, A(I,J)

       ENDDO

       ENDDO

       PRINT*, 'ENTER THE VALUE OF MATRIX B'

       DO I=1,3

       DO J=1,3

       READ*, B(I,J)

       ENDDO

       ENDDO

       CALL AL(W,Q,E,3,3,3)

       PRINT*, C(1,1),'   ',C(1,2),'   ',C(1,3)

       PRINT*, C(2,1),'   ',C(2,2),'   ',C(2,3)

       PRINT*, C(3,1),'   ',C(3,2),'   ',C(3,3)

       END