function [x,err,I] = Newton1(init,N);

% x=Newton(init,N); 
% 
% Approche par la méthode de Newton initiée à init le zéro réel du polynôme 
% x^5+ x^3/2 +1 compris entre - 1 et -1/2. 

x=init;
err=1;
I=1; 
while (err>eps) && (I<N) 
    y = x - (x.^5 + x.^3/2 +1)./(5*x.^4 + 3*x.^2/2);
    err = max(max(abs(y-x)));
    x=y;
    I=I+1;
end


