C语言OJ项目参考(2405)牛顿迭代法求根

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;
}