我是一个FORTRAN新手,现有一个迭代公式,想拿一次迭代结果与其前一次迭代结果进行对比如果对比结果达到要求,则迭代结束,否则继续迭代,式例(a-b)/a

Question

我是一个FORTRAN新手,现有一个迭代公式,想拿一次迭代结果与其前一次迭代结果进行对比
如果对比结果达到要求,则迭代结束,否则继续迭代,式例(a-b)/a

Answer 1

举个最简单的例子,Jacobi迭代解x+y=2,2x+3y=5
program Jacobi
integer::N,NMAX=10000
real::X,Y,X0,Y0,ERROR=0.0001
X0=0.0
Y0=0.0
do N=1,NMAX
X=2-Y0
Y=(5.0-2*X0)/3.0
if (abs((X-X0)/X0)