function  [c,nit]= dicho(f,a,b)
% dicho calcule un zéro de la fonction f dans l'intervalle [a,b]
% au moyen de la méthode de dichotomie
% juqu'à la précision machine
c=(a+b)/2;
nit=0;
while c>a & c<b,
  if feval(f,a)*feval(f,c) > 0,%alors la racine est dans ]c,b[
    a=c;
  else
    b=c;
  end;
c=(a+b)/2;
nit=nit+1;
end
