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

/注意最后一个测试点！当区间小于阈值时，输出边端点其一即可,写的时候一直不知道阈值是多少，后来才想明白，因为题中要求保留两位小数，所以阈值定为0.01即可/
代码如下：

#include<stdio.h>
#include<math.h>
double a1, a2, a3, a4;
double f(double x);
int main()
{
	 
	double a, b;
	scanf("%lf %lf %lf %lf %lf %lf", &a1, &a2, &a3, &a4, &a, &b);
	double x, c, d;
    c=(a+b)/2;
    if(f(a)*f(b)<0)
    {
        while(f(c)!=0&&b-a>=0.01)
        {
            if(f(c)*f(a)>0)
            {
                a=c;
                c=(a+b)/2;
            }
            else if(f(c)*f(b)>0)
            {
                b=c;
                c=(a+b)/2;
            }
        }
        if(f(c)==0)
            printf("%.2f",c);
        else 
            printf("%.2f",b);
    }
	else if(f(a)*f(b)==0)
    {
        if(f(a)==0) printf("%.2f",a);
        else if(f(b)==0) printf("%.2f",b);
    }
	return 0;
}
double f(double x)
{
	double aa = a1 * pow(x, 3) + a2 * pow(x, 2) + a3 * pow(x, 1) + a4;
	return aa;
}