下面代码实现了调换矩阵任意两列的功能
program main
implicit none
real :: x(4,4)
integer :: i, j, k
call random_seed()
call random_number(x)
do i = 1, size(x,1)
write(*,'(*(f7.3,3x))') x(i,:)
end do
write(*,*) '请输入要交换的列号j,k...'
read(*,*) j, k
x(:,j:k:k-j) = x(:,k:j:j-k)
do i = 1, size(x,1)
write(*,'(*(f7.3,3x))') x(i,:)
end do
end program main