牛顿迭代求方程的解

源程序：牛顿迭代法解非线性方程的根

！要解的方程为f(x)=x^2-16

program main

implicit none

real,parameter::eps=10**(-6.) ！设置一个精度

real::x0=10,x1,df,g ！df表示微分integer::k=1

integer,parameter::m=1000 ！设置一个循环次数的上限，避免不收敛导致死循环

do

df=2*x0

g=(x0**2+1.5)/(2*x0) ！g为迭代格式g=x - f(x)/f(x)' ，！f(x)'为导数，作为分母不能为零

if(df==0)then ！连续三重if结构嵌套write(*,*)"导数为零,奇异"

else

x1=g

if(abs(x1-x0)

write(*,*)"x=",x1

exit

else

if(k==m)then

write(*,*)"迭代失败" exit

else

k=k+1

x0=x1

end if

end if

end if

end do

end program