// Methode de Newton pour la resolution de f(x)=0
// Entree : f, df la fonction et sa derivee.
// xn+1 = xn −f(xn)/f'(xn)
// x0 approximation initiale
// epsi > 0 le test d'arret est |f(x)|< epsi

eps = 0.000001;
K   = 100;

x=x0 ;
k=1;
fx=f(x);
dfx=df(x); 
X=x0 ;
E=abs(fx);

while abs(fx) > eps & k<K do
 x=x-fx/dfx;
 fx=f(x);
 dfx=df(x);
 X=[X,x];
 E=[E,abs(fx)]
 k=k+1;
end
            
disp([k,X,E]);
          


