二分.cpp

来自「数值分析程序。。。。。好多程序下载 喜欢的朋友可以下 随时欢迎」· C++ 代码 · 共 31 行

#include <iostream.h>
#include <math.h>
double a,b,e;
double f (double x)
{
        double y;
	    y=pow(x,3)-x-1;
	    return y;
}
void main()
{
	cout<<"请输入有根跟区间a,b及误差限e"<<endl;
	cin>>a>>b>>e;
L:	double y0=f(a); 
	double x=(a+b)/2;
    double y=f(x);
    if(y*y0>0)
		a=x;
	else
		b=x;
	if(b-a>=e)
		goto L;
	else
	{   
		cout<<a<<endl;
	}
}