fortran实现jacobi迭代

下面是Fortran语言实现Jacobi迭代的代码：

program jacobi
implicit none

integer, parameter :: n = 3 ! 矩阵的维数
integer :: i, j, k, iter ! 循环变量和迭代次数
real :: a(n,n), b(n), x(n), xnew(n), eps ! 矩阵A、向量B、解向量X、新解向量Xnew、误差eps

! 初始化矩阵A和向量B
a = reshape([4., 1., 1., 1., 4., 1., 1., 1., 4.], [n,n])
b = [6., 6., 6.]

! 初始化解向量X和误差eps
x = [0., 0., 0.]
eps = 1.e-6

! Jacobi迭代
iter = 0
do while (maxval(abs(xnew - x)) > eps)
  iter = iter + 1
  x = xnew
  do i = 1, n
    xnew(i) = b(i)
    do j = 1, n
      if (i /= j) then
        xnew(i) = xnew(i) - a(i,j)*x(j)
      end if
    end do
    xnew(i) = xnew(i) / a(i,i)
  end do
end do

! 输出结果
print *, "迭代次数：", iter
print *, "解向量：", xnew

end program jacobi

这个程序实现了一个3x3的线性方程组的Jacobi迭代求解。在主程序中，我们首先定义了矩阵A和向量B，然后初始化解向量X和误差eps。接着进行Jacobi迭代，直到新解向量Xnew和旧解向量X的差的最大值小于误差eps为止。每次迭代中，我们用旧解向量X来计算新解向量Xnew，然后将Xnew赋值给X，继续下一次迭代。最后输出迭代次数和解向量Xnew。