Méthode de dichotomie

#include <stdio.h>
#include <stdlib.h>

double f(double x){
      double y;
      y=(x*x*x)-(x*x)+x+1;// la fonction x^3-x^2+x+1
      return y;}
void dichotomie(){
    double a,b,x;
    double E;// la précision epsilon

      printf("\ndonner l'interval [ a  b ] : ");
      scanf("%lf %lf",&a,&b);
      printf("donner E : ");
      scanf("%lf",&E);
do
{    
  x=(b-a)/2;
  if(f(a)*f(x)<0){
   b=x;
   
  }
   if(f(b)*f(x)<0)
  {
   a=x;
 }
} while(b-a > E);
printf(" %lf ",x);
 }



int main(int argc, char *argv[])
{
  dichotomie();

  system("PAUSE");
  return 0;
}