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

lsmCC

1. // 这题其实得用二分来做，为实数的二分
   
2. #include <bits/stdc++.h>
   
3. using namespace std;
   
4. double a,b,c,d;
   
5. double solve(double x)
   
6. {
   
7.     return a*x*x*x+b*x*x+c*x+d;
   
8. }
   
9. int main()
   
10. {
    
11.     cin>>a>>b>>c>>d;
    
12.     for(double i=-100;i<=100;i++)  //即在-100到100中一个一个区间去判断
    
13.     {
    
14.         double x1=i,x2=i+2.0;
    
15.         if(solve(x1)==0.0)printf("%.2lf",i); //如果带进去x1==0，则说明就是此方程的根
    
16.         if(solve(x1)*solve(x2)<0)//标明此处必有一个根，然后去判断
    
17.         {
    
18.             while(x2-x1>=0.001)
    
19.             {
    
20.                 double mid=(x1+x2)/2.0;
    
21.                 if(solve(x1)*solve(mid)<=0.0)x2=mid;
    
22.                 else x1=mid;
    
23. 
    
24. 
    
25.             }
    
26.             printf("%.2lf",x1);
    
27.         }
    
28. 
    
29.     }
    
30.     return 0;
    
31. 
    
32. 
    
33. }

复制代码

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

isdkz

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

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