function demo_newton_for

toll=10^(-6); nmax=1000; esempio=1;

% esempi.
switch esempio
    case 1 
        f=@(x) exp(x)-2+x; 
        f1=@(x) exp(x)+1;
        x0=1; sol=0.4428544010023885;
    case 2 
        f=@(x) sin(x)-x; 
        f1=@(x) cos(x)-1;
        x0=1; sol=0;
end

% newton
[xv,step,flag] = newtonfun_for(f,f1,x0,toll,nmax);

% statistiche
fprintf('\n \t soluzione : %1.15e',xv(end));
fprintf('\n \t tolleranza: %1.15e',toll);
fprintf('\n \t numero it.: %4d',length(xv));
if flag == 0
    fprintf('\n \t La procedura e'' terminata correttamente');
else
    fprintf('\n \t La procedura non e'' terminata correttamente');
end
fprintf('\n \n');
% ---------------------- plot ----------------------
indici=(1:length(step))';
semilogy(indici,step); 