基于visual Studio2013解决C语言竞赛题之0422牛顿迭代法

原创于 2013-12-03 15:28:25 发布 · 2.1k 阅读

1 ·

CC 4.0 BY-SA版权

VC++编程技术同时被 2 个专栏收录

780 篇文章

订阅专栏

VisualStduio2013C语言系列

180 篇文章

订阅专栏

本文介绍如何使用牛顿迭代法在数值范围内求解立方多项式方程的根，并通过实例展示了代码实现过程。重点强调了浮点数处理中的精度问题，确保结果的准确性。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

题目

解决代码及点评

该题考察对浮点数的了解

#include<iostream>
#include<cmath>
/*
用牛顿迭代法求方程在 1.5附近的根(精度在10^-5)
2*x^3-4*x^2+3*x-6=0
*/
using namespace std;
void main()
{
	double diedai(double a,double b,double c,double d,double x);
	double a,b,c,d;
	double x=10000.0;
	cout<<"请依次输入方程四个系数：";
	cin>>a>>b>>c>>d;
	x=diedai(a,b,c,d,x);      // 用户输入四个系数后，调用迭代法就根
	cout<<x<<endl;
	//printf("%f",x);
	system("pause");
}
double diedai(double a,double b,double c,double d,double x)
{
	while(abs(a*x*x*x+b*x*x+c*x+d)>0.000001) // 关键在于，对方程的求值，最后是>0.000001而不是==0是因为浮点数不能==0，所以用精度代替
	{
		x=x-(a*x*x*x+b*x*x+c*x+d)/(3*a*x*x+2*b*x+c);
	}
	return x;
}