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

二分法求函数根的原理为：如果连续函数(在区间[的两个端点取值异号，即(，则它在这个区间内至少存在1个根$r$，即(。

二分法的步骤为：

• 检查区间长度，如果小于给定阈值，则停止，输出区间中点(；否则
• 如果(，则计算中点的值(；
• 如果(正好为0，则(就是要求的根；否则
• 如果(与(同号，则说明根在区间[，令)，重复循环；
• 如果(与(同号，则说明根在区间[，令)，重复循环。

本题目要求编写程序，计算给定3阶多项式(在给定区间[内的根。

输入格式：

输入在第1行中顺序给出多项式的4个系数a​3​​、a​2​​、a​1​​、a​0​​，在第2行中顺序给出区间端点$a$和$b$。题目保证多项式在给定区间内存在唯一单根。

输出格式：

在一行中输出该多项式在该区间内的根，精确到小数点后2位。

输入样例：

3 -1 -3 1
-0.5 0.5

输出样例：

0.33

我楞没看出来我的代码跟网上找的有什么不同

下面是我的代码，结果是用例1和用例3死活不通过。。。。

测试点	提示	结果	耗时	内存
0	sample 非整数根	答案正确	2 ms	128KB
1	二分正好在中点	答案错误	2 ms	128KB
2	区间端点是根	答案正确	1 ms	340KB
3	中点的f值很小但不一定是根	答案错误	2 ms	128KB

#include<stdio.h>
int a3, a2, a1, a0;
float f(float g);

int main()
{

  float a, b, r;
  scanf("%d %d %d %d", &a3, &a2, &a1, &a0);
  scanf("%f %f", &a, &b);
  r = f((a+b)/2);
  if(f(a)==0)
    {
        b = a;
    }
    if(f(b)==0)
    {
        a = b;
    }
    while((b-a)>0.001 && f(a)*f(b)<0)
  {
    //printf("a:%f\tb:%lf\n", a, b);
    //printf("(a+b)/2:%f, r:%f\n",(a+b)/2.0, r);
    
    if(r == 0)
    {
        break;
    }
    if(r*f(a)<0)
    {
        b = (a+b)/2;
    }
    if(r*f(b)<0)
    {
        a = (a+b)/2;
    }
      r =  f((a+b)/2);
   }
  printf("%.2f", (a+b)/2);
  system("pause") ;
}

float f(float g )
{
  float rlt;
  rlt = a3*g*g*g + a2*g*g + a1*g + a0;
  return rlt;
}

下面是网上找的代码，一下子就通过了

#include <stdio.h>
float f(float x);
float a3, a2, a1, a0;


int main()
{
    float a,b,c;
    scanf("%f%f%f%f",&a3,&a2,&a1,&a0);
    scanf("%f%f",&a,&b);
    if(f(a)==0)
        c=a;
    if(f(b)==0)
    c=b;
    while(f(a)*f(b)<0&&b-a>=0.001)
    {
       c=(a+b)/2;
       if(f(c)*f(a)>0)
        a=c;
       if(f(c)*f(b)>0)
        b=c;
       if(f(c)==0)
       break;
    }


    printf("%0.2f",c);
    system("pause");
}
float f(float x)
{
    float result;
    result = a3*x*x*x + a2*x*x + a1*x + a0;
    return result;
}