
% nome: cognome: matricola: postazione: 

function demo_II

% dati problema
f=@(x) log(x)-3;
f1=@(x) 1./x;
x1=10;
z1=30;
toll=10^(-12);
nmax=1000;

% metodo iterativo
[xv,zv] = newton_fourier(f,f1,x1,z1,toll,nmax);

% stampa risultati finali
fprintf('\n \t xv(end): %1.15e zv(end): %1.15e err: %1.3e \n \n',xv(end),...
    zv(end),abs(xv(end)-zv(end)));

% grafici
semilogy(1:length(xv),abs(xv-zv));

% salvataggio figura su file
saveas(gcf, 'errore', 'jpg')

%% RISULTATI:
% >>demo_II
% 
% xv(end): 2.008553692318767e+01 zv(end): 2.008553692318767e+01 err: 3.553e-15 
%  
% >>
%
% % NOTA: La soluzione e' ovviamente exp(3)
% format long e; exp(3)
% 
% ans =
% 
%      2.008553692318767e+01

