一元三次方程求解用的二分，为什么结果不对

lsmCC

// 这题其实得用二分来做，为实数的二分
#include <bits/stdc++.h>
using namespace std;
double a,b,c,d;
double solve(double x)
{
    return a*x*x*x+b*x*x+c*x+d;
}
int main()
{
    cin>>a>>b>>c>>d;
    for(double i=-100;i<=100;i++)  //即在-100到100中一个一个区间去判断
    {
        double x1=i,x2=i+2.0;
        if(solve(x1)==0.0)printf("%.2lf",i); //如果带进去x1==0，则说明就是此方程的根
        if(solve(x1)*solve(x2)<0)//标明此处必有一个根，然后去判断
        {
            while(x2-x1>=0.001)
            {
                double mid=(x1+x2)/2.0;
                if(solve(x1)*solve(mid)<=0.0)x2=mid;
                else x1=mid;


            }
            printf("%.2lf",x1);
        }

    }
    return 0;


}

其中有一个地方不懂x2-x1>=0.001

isdkz

这段代码中的 x2-x1>=0.001 是用来控制二分法的迭代次数的条件。在二分法中，我们希會在每次迭代中减小区间的长度，直到满足精度要求。这里的0.001是一个预先设定的精度值，表示当区间长度小于等于0.001时，就认为已经找到了足够接近的解。因此，当 x2-x1 的值小于等于0.001时，就会退出迭代，输出当前的 x1 值作为方程的解。

然而，在实际应用中，这种固定的迭代终止条件并不总是适用的。在一些情况下，可能需要根据具体问题来调整迭代的终止条件，以确保找到足够精确的解。