2405: 牛顿迭代法求根

Description

用牛顿迭代法求根。方程为 ax3+bx2+cx+d=0 。系数a,b,c,d的值一次为1,2,3,4,由主函数输入。求x在1附近的一个实根。求出根后由主函数输出。结果保留两位小数。

Input

系数a,b,c,d的值

Output

x在1附近的一个实根

Sample Input

1 2 3 4

Sample Output

-1.65

HINT

主函数已给定如下,提交时不需要包含下述主函数

/* C代码 */
int main()
{
double solut(double ,double ,double ,double );
double a,b,c,d;
scanf("%lf%lf%lf%lf",&a,&b,&c,&d);
printf("%.2f",solut(a,b,c,d));;
return 0;
}

参考解答

#include<stdio.h>
#include<math.h>

int main()
{
double solut(double,double,double,double );
double a,b,c,d;
scanf("%lf%lf%lf%lf",&a,&b,&c,&d);
printf("%.2f",solut(a,b,c,d));;
return 0;
}

double solut(double a,double b,double c,double d)
{
double x=1,x0;
do
{
x0=x;
x=x0-(((a*x+b)*x+c)*x+d)/((3*a*x+2*b)*x+c);
}
while(fabs(x-x0)>=1e-5);
return x;
}