function x=dichot(init1,init2,N);

% x=dichot(init1,init2,N);
%
% approche par la méthode de dichotomie le zéro réel de 
% la fonction polynomiale 
% x -- > x^5 + x^3/2 + 1 
% (compris entre -1 et -1/2).  

x1=init1;
x2=init2;
for i=1:N 
   y=(x1+x2)/2;
   u=(x1^5 + x1^3/2 +1)*(y^5 + y^3/2 +1); 
   if u < 0
   x1=x1;
   x2=y;
   else 
       x1=y;
       x2=x2; 
   end
end
x=x2;
