![](https://csdnimg.cn/release/download_crawler_static/88819213/bg6.jpg)
stop
end
(
2
)对于
3 2
94 389 294 0x x x+ - + =
,
按照上述数学原理
,
编制的程序如下
program newton
implicit none
real :: x(0:50),fx(0:50),f1x(0:50)!分别为自变量 x,函数 f(x)和一阶导数 f1(x)
integer :: k
write(*,*) ”x(0)=”
read(*,*) x(0) !输入变量:初始值 x(0)
open(10,file='1。txt’)
do k=1,50,1
fx(k)=x(k—1)**3+94*x(k—1)**2—389*x(k—1)+294
f1x(k)=3*x(k—1)**2+188*x(k—1)-389
x(k)=x(k-1)—fx(k)/f1x(k) !牛顿法
write(*,'(I3,1x,f11。6)') k,x(k) !输出变量:迭代次数 k 及 x 的值
write(10,'(I3,1x,f11。6)') k,x(k)
if(abs(x(k)—x(k-1))<5e-6) exit !终止迭代条件
end do