简单二分法，为什么会死循环呢？？

[firerun]

#include <iostream>
#include <math.h>
using namespace std;

class Bifid{
	private:
		double x1, x2;

	public:
		Bifid(double a, double b) { x1=a; x2=b;}
		
		double fun(double x){
			return 1.0*pow(x,2) + 0.0*x - 5.0;
			}			

		double erfenfa() {
			if (fun(x1)==0) return x1;
			if (fun(x2)==0)	return x2;
			if (fun(x1)*fun(x2) > 0) {
				cout << "wu jie" << endl;
				return 0;	}
		
			static double mid=(x1+x2)/2.0;
					
			while(fabs(fun(mid))>1e-10) {
				cout << x1 << ' ' << x2 << endl;  //检测程序另外加的
				mid=(x1+x2)/2.0;
			  if(fun(x1)*fun(mid) < 0 )	x1=mid;			
			  else x2=mid; }	
			return mid;	
			}	

	};

int main()
{
	Bifid a(0.0, 200.0);
	cout << a.erfenfa() << endl;
	
	return 0;
}
	

[Jncryang]

#include <iostream>
#include <math.h>
using namespace std;

class Bifid{
    private:
        double x1, x2;

    public:
        Bifid(double a, double b) { x1=a; x2=b;}
        
        double fun(double x){
            return 1.0*pow(x,2) + 0.0*x - 5.0;
            }            

        double erfenfa() 
		{
            if (fun(x1)==0) return x1;
            if (fun(x2)==0)    return x2;
            if (fun(x1)*fun(x2) > 0) 
			{
                cout << "wu jie" << endl;
                return 0;    
			}
        
            static double mid=(x1+x2)/2.0;
                    
            while(fabs(fun(mid)) < 1e-10) 
			{
                cout << x1 << ' ' << x2 << endl;  //检测程序另外加的
                mid=(x1+x2)/2.0;
				if(fun(x1)*fun(mid) < 0 )    
					x1=mid;            
				else x2=mid; 
			}    
            return mid;    
		}    

    };

int main()
{
    Bifid a(1.0, 200.0);
    cout << a.erfenfa() << endl;
    return 0;
}

循环始终为真导致死循环！

[firerun]

可是while这块我也检查不出什么问题啊？？

            while(fabs(fun(mid))>1e-10) {
                cout << x1 << ' ' << x2 << endl;  //检测程序另外加的
                mid=(x1+x2)/2.0;
              if(fun(x1)*fun(mid) < 0 )    x1=mid;            
              else x2=mid; }    
            return mid;    
            }    

[birchlin]

有没有比二分发更快的查询方法?

[zxwl]

在看代码

不过首先：实数比较是否为0，不能：
==0

应该：
//精度有你来定
const int BOND = 0.000001
//......
if（num > -BOND && num < BOND）//这样判断

[zxwl]

抱歉
const float BOND = 0.000001

[firerun]

[Quote=引用 2 楼 firerun 的回复:]
可是while这块我也检查不出什么问题啊？？

C/C++ code while(fabs(fun(mid))>1e-10) {
cout << x1 << ' ' << x2 << endl; //检测程序另外加的
mid=(x1+x2)/2.0;
if(fun(x1)*fun(mid) > 0 ) x1=mid; //##小于0改为大于0就ok
else x2=mid; }
return mid;
}
[/Quote]

[zwspider]

float fabs(float x); // float 的精度是7

while(fabs(fun(mid))>1e-10) // 条件可能为FALSE吗？死循环

建议你自己写个double的fabs()