1.问题:编写用牛顿迭代法求方程根的函数。方程为a*x*x*x+b*x*x+c*x+d=0,系数a,b,c,d由主函数输入。求x在1附近的一个实根。求出根后,有主函数输出。
牛顿迭代法的公式是:x=x0-f(x0)/f'(x0),设迭代到|x-x0|<=1/100000结束。
2.思路:定义好x与x0两个变量,通过do-while循环进行迭代,每次循环都通过x=x0-f(x0)/f'(x0)公式得到的x再次赋值给x0,并判断|x-x0|<=1/100000的条件是否成立,最终输出x的值。
3.程序流程图:
4.代码实现:
#include <iostream>
#include<math.h>
using namespace std;
int main()
{
float x=1.5,x0;
int a,b,c,d;
cin>>a>>b>>c>>d;
do{
x0=x;
float f=a*x0*x0*x0+b*x0*x0+c*x0+d;
float f1=3*a*x0*x0+2*b*x0+c;
x=x0-f/f1;
}while(fabs(x-x0)>=0.00001);
cout<<x;
return 0;
}