二分法.cpp

数值分析的实验报告

#include<iostream>
#include <iomanip>
#include <cmath>
#define f(x) (x*x*x-x-1)
using namespace std;
int main()
{
	int exact;
	cout<<"请输入您要的精度值(结果的有效小数位):"<<endl;
	cin>>exact;
	double e=pow(10,-exact);
	double x,a=1,b=2;
	for(int i=0;;i++,a=b)
	{
		x=(a+b)/2;
		if(f(x)*f(b)<0)
			a=x;
		else
			b=x;
		if(fabs(a-b)<e)
			break;
	cout<<"经过第"<<i+1<<"次运算:"<<"x"<<i<<"="
			<<setiosflags(ios::fixed)<<setprecision(exact)<<a
			<<",x"<<i+1<<"="
			<<setiosflags(ios::fixed)<<setprecision(exact)<<b
			<<endl;
	}
	cout<<"共运算了"<<i<<"次,结果为:";
	
	cout<<setiosflags(ios::fixed)<<setprecision(exact)<<x<<endl;
	return 0;
}