面试题：Fortran多维数组的复杂操作

Question

假设有一个二维数组`B(5, 5)`，代表一个5x5的矩阵。编写Fortran程序，计算该矩阵的转置矩阵，并将结果存储在另一个二维数组`C`中，最后输出转置后的矩阵。

Answer 1

program matrix_transpose
    implicit none
    integer :: i, j
    integer, dimension(5, 5) :: B
    integer, dimension(5, 5) :: C

   ! 初始化矩阵B
    B = reshape((/1, 2, 3, 4, 5, &
                  6, 7, 8, 9, 10, &
                  11, 12, 13, 14, 15, &
                  16, 17, 18, 19, 20, &
                  21, 22, 23, 24, 25/), [5, 5])

   ! 计算转置矩阵
    do i = 1, 5
        do j = 1, 5
            C(i, j) = B(j, i)
        end do
    end do

   ! 输出转置后的矩阵
    do i = 1, 5
        do j = 1, 5
            write(*, '(I3)', advance = 'no') C(i, j)
        end do
        write(*, *)
    end do
end program matrix_transpose