二分法求方程根有感而发

/*本题目要求编写程序，计算给定3阶多项式在给定区间[a,b]内的根。
输入格式：
输入在第1行中顺序给出多项式的4个系数 
在第2行中顺序给出区间端点a和b。题目保证多项式在给定区间内存在唯一单根。
输出格式：
在一行中输出该多项式在该区间内的根，精确到小数点后2位。*/ 
#include <stdio.h>
double a3,a2,a1,a0,a,b; //定义全局变量方便f函数使用 
double f(double i)
{
    return a3*i*i*i+a2*i*i+a1*i+a0;
}
int main()
{
    scanf("%lf%lf%lf%lf",&a3,&a2,&a1,&a0);
    scanf("%lf%lf",&a,&b);
    double mid=a+(b-a)/2;
    //先判断两点是否为根 
    if(f(a)==0)
    {
        printf("%.2lf",a);
        return 0;
    }
    else if(f(b)==0)
    {
        printf("%.2lf",b);
        return 0;
    }
    //题目的小暗示：要求的精度为小数点后两位，所以阈值为0.001或0.0001
	//即使中值的f很小始终不等于0，从while出循环后。 
	//因为阈值为0.001，f==0处小数点前两位已经确定，无需再二分，符合答案要求 
    while(a<b-0.001) //重点！！ 
    {
        mid=a+(b-a)/2;
        if(f(mid)==0)
            break;
        if(f(mid)*f(a)<0)
            b=mid;
        else if(f(mid)*f(b)<0)
            a=mid;
    }
    printf("%.2lf",mid);
    return 0;
}