%La méthode de Newton pour une application f : [u,v]->R
%va fournir un couple (r,n), en partant d'un point a de [u,v],
%où r est proche d'un zéro de f et n est le nombre d'itérations utilisées
%Appel de cette fonction par : [x,n]=newton('f','fprime',2)
function [x,n] = newton(f,fprime,a)

x=a;
tol=0.000001; %tolérance
nmax=10000;   %nombre maximum d'itérations

for n=1:nmax
   y = x-feval(f,x)/feval(fprime,x);
      if (abs(y-x) <= tol),
         break
      end
   x=y;
end