使用二分法求解f(x)=x^3-x-1=0在区…

//使用二分法求解f(x)=x^3-x-1=0在区间(1，2)上的解
#include<iostream.h>
#include<math.h>
const double eps=0.000001; //对于精度的控制
//为方便的代码的书写，故编写此函数，以此求得此函数值
double f(double x)
{
 return pow(x,3)-x-1;  //当函数发生改变时，修改此处即可
}
//解的具体实现的函数
void solequ(double min,double max,double eps)
{
 double temp;
 if(f(min)*f(max)<0)//判断方程在此区间是否有解
 {     //有解的情况的求解过程
  while(f(max)!=0)
  {
   temp=(min+max)/2;
   if(f(temp)*f(min)<0)
   {
    max=temp;
    if(fabs(f(max)-f(min))<eps)
     break;
   }
   else
   {
    min=temp;
    if(fabs(f(max)-f(min))<eps)
     break;
   }
  }
  cout<<"方程的解为: "<<temp<<endl;
 
 }
 else   //无解的处理
  cout<<"此方程在此区间上无解"<<endl;

}
//主函数入口
int main()
{
 double min,max,temp;
 cout<<"请输入积分下限: ";
 cin>>min;
 cout<<"请输入积分上限: ";
 cin>>max;
 if(min>max)  //对于用户输入次序相反的处理
 {
  temp=max;
  max=min;
  min=temp;
 }
 solequ(min,max,eps);
 return 0;
 
}