我是一个FORTRAN新手,现有一个迭代公式,想拿一次迭代结果与其前一次迭代结果进行对比
发布网友
发布时间:2022-05-18 03:48
我来回答
共1个回答
热心网友
时间:2023-10-08 23:02
举个最简单的例子,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)<ERROR.and.abs((Y-Y0)/Y0)<ERROR) then
write(*,*)"迭代次数:",N
write(*,*)"X=",X," Y=",Y
exit
end if
X0=X
Y0=Y
end do
if (N==NMAX+1) then
write(*,*) "没得到结果!"
end if
pause
end program Jacobi