一元三次方程求解用的二分,为什么结果不对
// 这题其实得用二分来做,为实数的二分#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 这段代码中的 x2-x1>=0.001 是用来控制二分法的迭代次数的条件。在二分法中,我们希會在每次迭代中减小区间的长度,直到满足精度要求。这里的0.001是一个预先设定的精度值,表示当区间长度小于等于0.001时,就认为已经找到了足够接近的解。因此,当 x2-x1 的值小于等于0.001时,就会退出迭代,输出当前的 x1 值作为方程的解。
然而,在实际应用中,这种固定的迭代终止条件并不总是适用的。在一些情况下,可能需要根据具体问题来调整迭代的终止条件,以确保找到足够精确的解。
页:
[1]