function [x,k]=newton(f,fp,x0,epsilon,Nmax)
// 
//  Methode de Newton selon la formule d appel
//  [x,k]=newton(f,fp,x0,epsilon,Nmax)
//  
xp = x0;
x = [xp];
fxp = f(xp);
k = 0;
while (abs(fxp) > epsilon) & (k <= Nmax)
  fprimxp = fp(xp);
  xp = xp - fxp / fprimxp;
  x = [x, xp];
  fxp = f(xp);
  k=k+1;
end
endfunction