%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