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

来源:互联网 发布:cm4d软件 编辑:程序博客网 时间:2024/05/16 01:21

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