fortran 对一个方形矩阵进行转置

Fortran中可以使用do循环和临时变量来对一个方形矩阵进行转置。下面是一个简单的示例程序：

program transpose_matrix
  implicit none
  
  integer :: i, j
  integer, dimension(3,3) :: my_matrix
  integer :: temp
  
  ! 初始化矩阵
  my_matrix = reshape([1,2,3,4,5,6,7,8,9], [3,3])
  
  ! 打印原始矩阵
  write(*,*) "Original Matrix:"
  do i = 1, 3
    write(*,*) my_matrix(i,:)
  end do
  
  ! 转置矩阵
  do i = 1, 3
    do j = i+1, 3
      temp = my_matrix(i,j)
      my_matrix(i,j) = my_matrix(j,i)
      my_matrix(j,i) = temp
    end do
  end do
  
  ! 打印转置后的矩阵
  write(*,*) "Transposed Matrix:"
  do i = 1, 3
    write(*,*) my_matrix(i,:)
  end do
  
end program transpose_matrix

在上面的示例程序中，我们使用reshape函数初始化了一个3x3的矩阵`my_matrix`。然后，我们使用两个嵌套的do循环和一个临时变量`temp`来对矩阵进行转置。最后，我们使用一个循环来打印转置后的矩阵。