牛顿迭代法、迭代逼近

本文主要是介绍牛顿迭代法、迭代逼近，希望对大家解决编程问题提供一定的参考价值，需要的开发者们随着小编来一起学习吧！

/*zstu  求一元三次方程的根

要求1附近的根，只要把x0初始化为1代入，然后牛顿迭代1000次就出结果了

*/

#include<stdio.h> double c[4],d[3]; double f(double x) { double fac=1,sum=0; int i; for(i=0;i<4;i++) { sum+=fac*c[i]; fac*=x; } return sum; } double df(double x) { double fac=1,sum=0; int i; for(i=0;i<3;i++) { sum+=fac*d[i]; fac*=x; } return sum; } main() { int i; double x0; while(scanf("%lf%lf%lf%lf",&c[3],&c[2],&c[1],&c[0])!=EOF) { d[0]=c[1]; d[1]=c[2]*2; d[2]=c[3]*3; x0=1; i=1; while(i<1000) { i++; x0=x0-f(x0)/df(x0); } printf("%.2f/n",x0); } } 　

这篇关于牛顿迭代法、迭代逼近的文章就介绍到这儿，希望我们推荐的文章对编程师们有所帮助！

---