clear

%-------------------------------------------------
% Questo M-file implementa il metodo di Newton
% Le funzioni fun1 e fun1d contengo la definizione
% della funzione e della derivata, rispettivamente.
%-------------------------------------------------
x0=input('Dammi x0 = ');

k=0;
tol=1.e-10;
kmax=100;
xx(1)=x0;

x1=x0-fun1(x0)/fun1d(x0);
xx(2)=x1;
while abs((x1-x0)) > tol*abs(x1) & k <= kmax,
    x0=x1;
    x1=x0-fun1(x0)/fun1d(x0);
    k=k+1;
    xx(k+2)=x1;
end
disp('La radice cercata e` : ')
x0
disp('Numero iterazioni effettuate : ')
k
disp('Sequenza delle iterate ');
xx

