二分法求多项式根

#include <stdio.h>
#include <math.h>
int a[4];
  double c[2];
double fun(double x)
{
  return a[0]*x*x*x+a[1]*x*x+a[2]*x+a[3];
}




int main(int argc, char* argv[])
{
          int i;
  double x=10,x1=1,f1=1;


for( i=0;i<4;i++)
{
scanf("%d",&a[i]);
}
  for( i=0;i<2;i++)
{
scanf("%lf",&c[i]);
}


  while(fabs(f1)>(1.0e-5))
  {
   x=x1;


   if(fun(c[0])*fun(c[1])<0)
   {
                   x1=(c[0]+c[1])/2;
                         f1=fun(x1);






   }
   if(f1*fun(c[0])>0)
           c[0]=x1;
     else if(f1*fun(c[1])>0)
                 c[1]=x1;


  }
    printf("%0.3f",x1);  //精确到小数点3位


        return 0;
}