新手求教 C++算牛顿二乘法的根出现的问题

[qq_41580767]

#include<iostream.h>
void main()
{
float x,x0,f,f1;
x=1;
do
{
x0=x;
f=x*x*x*x*x-3*x*x+1;
f1=5*x*x*x*x-6*x;
x=x0-f/f1;
}while(x-x0<0.000001);
cout<<"x="<<x<<endl;
}
我是这么写的，但是不知道为什么算不出结果，求大佬指教一下

[幻夢之葉]

如果没有输出任何东西，就是算法有问题，造成了死循环

[日立奔腾浪潮微软松下联想]

第二次迭代的时候f1就成0了，所以 x=x0-f/f1;这里就出现除以0异常。

[qq_41580767]

请问一下应该怎么改呢？ 我之前用int定义的话 会出现我的cpp运行的时候会停止运行 但是用float double的话又什么结果都没有输出出来……

[paschen]

自己要把算法逻辑理清楚，必要时单步跟踪程序运行，观察程序流程及变量变化情况

[日立奔腾浪潮微软松下联想]

引用 3 楼 qq_41580767 的回复:

请问一下应该怎么改呢？ 我之前用int定义的话 会出现我的cpp运行的时候会停止运行 但是用float double的话又什么结果都没有输出出来……

这应该是你用的公式导致的，可以把 x=x0-f/f1; 改成 x = x0 - f / (f1 + 0.0000001); 试试

[qq_41580767]

x = x0 - f / (f1 + 0.0000001); 改成这样后X=1.14286e+006 这是赋值上的问题吗

[日立奔腾浪潮微软松下联想]

稍微改进一下： if (f1 == 0) x = x0 - (f + 0.0000001) / (f1 + 0.0000001); else x = x0 - f / f1; 说到底是你用的这公式造成的，也就这样了。

[cg111]

[code=c] [#include<iostream.h> void main() { using namespace std; float x,x0,f,f1f2; x=1; do { x0=x; f=x*x*x*x*x-3*x*x+1; f1=5*x*x*x*x-6*x; f2=5*4*x*x*x-6 x=x0-f1/f2; }while(x-x0<0.000001); cout<<"x="<<x<<endl; }/code]

[cg111]

牛顿迭代法公式“
X(n+1)=X(n)-F'(Xn)/F"(Xn)
你的代码迭代公式不对