5-18 二分法求多项式单根 (20分)

#include <cstdio>
#include <iostream>
#include <cstring>
#include <cstdlib>
#include <algorithm>
#include <map>
using namespace std;
#define INF 0x3f3f3f3f

float a0,a1,a2,a3;
float fun(float t)
{
    float s;
    s=a3*pow(t,3.0)+a2*pow(t,2.0)+a1*t+a0;
    return s;
}
int main()
{
    float a,b,x;
    cin>>a3>>a2>>a1>>a0;
    cin>>a>>b;
    while(b-a>=0.001)
    {
        x=(a+b)/2;
        if(fun(x)==0)
        {
            printf("%.2f\n",x);
            return 0;
        }
        else if(fun(x)*fun(a)>0)
            a=x;
        else
            b=x;

    }
    printf("%.2f\n",x);
    return 0;
}