二分法求函数的零点

#include<stdio.h>

#include<math.h>

double target = 1e-7;

double f(double x){

    return pow(x,5)-15*pow(x,4)+85*pow(x,3)-225*pow(x,2)+274*x-121;

}

double fun(double a,double b){

    double mid;

    while(a<b){

        mid=(a+b)/2.0;

        if(fabs(f(mid))<target)

         return mid;

        if(f(mid)*f(a)>0){

            a=mid;

        }else

          b=mid;

    }

    return (a+b)/2.0;

}

int main()

{

      double a=1.5;

      double b=2.4;

      printf("%.6f", fun(a,b));

return 0;

}