一元三次方程求解

#include<iostream>
#include<cstdio>
using namespace std;
double a,b,c,d;
double f(double x)
{
    return a*x*x*x+b*x*x+c*x+d;
}
int main()
{
        cin>>a>>b>>c>>d;
    double l,r,x1,x2,mid;

    int s=0,i;
    for(i=-100;i<100;i++)
    {
        l=i;
        r=i+1;
        x1=f(l);
        x2=f(r);
        if(!x1)
        {
            printf("%.2lf ",l);
            s++;
        }
        if(x1*x2<0)
        {
            while(r-l>=0.001)
            {
                mid=(r+l)/2;
                if(f(l)*f(mid)<0)r=mid;
                else l=mid;
            }
            printf("%.2lf ",l);
            s++;
        }
        if(s==3)break;
    }
    return 0;
}