首页 > 其他好文 > 详细

OJ刷题之《牛顿迭代法求根》

时间：2014-12-30 20:44:20 阅读：231 评论：0 收藏：0 [点我收藏+]

标签：c++ iostream namespace 编程迭代

题目描述

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

输入

系数a,b,c,d的值

输出

x在1附近的一个实根

样例输入

1 2 3 4

样例输出

-1.65

提示

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

/* 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;

}

/* C++代码 */

int main()

{

double solut(double ,double ,double ,double );

double a,b,c,d;

cin>>a>>b>>c>>d;

cout<<setiosflags(ios::fixed);

cout<<setprecision(2);

cout<<solut(a,b,c,d)<<endl;

return 0;

}

代码如下：

#include <iostream>
#include <iomanip>
#include <cmath>
using namespace std;
double solut(double a,double b,double c,double d)
{
    double x=1,x0,f,f1;
    do
    {
        x0=x;
        f=((a*x0+b)*x0+c)*x0+d;
        f1=(3*a*x0+2*b)*x0+c;
        x=x0-f/f1;
    }while (fabs(x-x0)>=1e-2);
    return x;
}
int main()

{

 double solut(double ,double ,double ,double );

 double a,b,c,d;

 cin>>a>>b>>c>>d;

 cout<<setiosflags(ios::fixed);

 cout<<setprecision(2);

 cout<<solut(a,b,c,d)<<endl;

 return 0;

}

运行结果：

技术分享