% Méthode de dichotomie 
clc
clear all
x=1:0.1:10;
%f=inline('cos(x)-x.^3');          % exp(d)
f=inline('x-exp(sin(x))');         % exp(a)
%f=inline('x^3-12*x+1');           % exp(c)
%f=inline('x.^3-12*x.^2-60*x+46'); % exp(b)
plot(x,f(x)), grid 
a=1; fa=f(a);
b=10; fb=f(b);
i=0;
eps=1.0e-3;
if ((fa*fb)<0)
while (b-a) > eps
  x=(a+b)/2; 
  i=i+1;
  if (sign(f(x)) == sign(fa))
     a=x; fa=f(x);
  else
     b=x; fb=f(x);
  end
  fprintf('dans l''ititation i=%d \t la solution est x0=%f \t f(x0)= %f \n',i,x,f(x))
end;
    fprintf('La solution finale est x0 = %f \n',x)
    
else
    disp('"On ne peut pas faire de dichotomie dans cet intervalle !!')
end