転置行列と行列の積の求め方

Question

学校でフォートラン77を使用しています。 Aの転置行列とBの行列の積のD行列を作りたいのですが、以下のプログラムをどのように変えれば良いでしょうか。具体的に教えてもらえると嬉しいです

C2345678
dimension A(100,100), B(100,100), C(100,100)
*
read(,) N, L, M
*
write(,) '**** [A] '
do i = 1, N
read(,) (A(i,j),j=1,L)
end do
*
write(,) '** [B] '
do i = 1, L
read(,) (B(i,j),j=1,M)
end do
*
call AMB( A,N,L,B,M,C )
*
*
write(,) '** [C] =[A][B] **'
do i = 1, N
write(,) (c(i,j),j=1,M)
end do
*
stop
end

---

• SUBROUTINE AMB

---

SUBROUTINE AMB( A,N,L,B,M,C )
DIMENSION A(100,100), B(100,100), C(100,100)
*
DO I = 1, N
DO J = 1, M
C(I,J) = 0.0
DO K = 1, L
C(I,J) = C(I,J) + A(I,K)*B(K,J)
END DO
END DO
END DO
*
RETURN
END