二分逼近/牛顿迭代——一元高次非线性方程求解

#include"iostream"
#include"cstdio"
#include"cstring"
#include"cstdlib"
#include"cmath"
#define precision 0.00000001

using namespace std;

int count=0;

//一元多次方程求解
double function(double x)
{
	return 2*x*x+3.2*x-1.8;
}

double binary(double x,double y)
{
	double mid=(x+y)/2.0;
	while(fabs(x-y)>precision)
	{
		if(function(mid)==0) break;
		if(function(x)*function(mid)<0) y=mid;
		else x=mid;
		mid=(x+y)/2;
		count++;
	} 
	return mid;
}

double function2(double x)   //返回导数,高精度代倒数 
{
	return (function(x+0.000005)-function(x-0.000005))/0.00001; 
}

double newtown(double x)
{
	double x1=x-function(x)/function2(x);
	while(fabs(x1-x)>precision)
	{
		if(function(x1)==0) break;
		x=x1;
		x1=x-function(x)/function2(x);
		count++;
	}
	return x1;
}

int main()
{
	cout<<binary(-0.8,8.0)<<endl;
	cout<<count<<endl;
	count=0;
	cout<<newtown(8)<<endl;
	cout<<count<<endl;
	count=0;
	cout<<newtown(-8)<<endl;
	cout<<count<<endl;
	return 0;
} 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值