
function demo_sequenza

% Demo di "calcola_sequenza" che la valuta relativamente ai parametri 
% "x1=1" e "N=71".

x1=1; N=71;

% calcola "sequenza", ovvero "x_1,...,x_100"
y=calcola_sequenza(x1,N);

fprintf('\n \t lunghezza(y): %4.0f   N: %4.0f',length(y),N);

% ultimo valore di "y"

% L=length(y);
% S=y(L);

S=y(end);

fprintf('\n \t S: %1.15e',S);

% valore soluzione
sol=log(1+x1);

% valutazione errore 
err=abs(y-sol);

% stampa "err"
fprintf('\n \t err(end): %1.5e',err(end));
fprintf('\n \n');

% grafico
clf;
semilogy(1:N,err,'ko-','LineWidth',2,...
                       'MarkerEdgeColor','k',...
                       'MarkerFaceColor','g',...
                       'MarkerSize',4)

% COMMENTO: 
% il fenomeno di cancellazione fa si che da un certo punto in poi la
% parentesi (sqrt(...)-1) abbia valore 0.

% RISULTATO.
% >> demo_sequenza
% 
%  	 lunghezza(y):   71   N:   71
%  	 S: 0.000000000000000e+00
%  	 err(end): 6.93147e-01 
% >>