
% nome: cognome: matricola: postazione: 

function demo_I

% dati del problema
f=@(x) x.^4-10;
f1=@(x) 4*x.^3;
x1=1;
z1=2;
toll=10^(-6);
nmax=1000;

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

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


%% APPLICAZIONE DI demo_I
% >> demo_I
% 
%  	 xv(end): 1.778279410038923 zv(end): 1.778279402556865 err: 7.48e-09 
%  
% >>