
function demo_derivata_terza

% funzione (vettoriale)
f=@(x) exp(x)-x.^3;

% derivata terza (vettoriale)
D3f=@(x) exp(x)-6;

% intervallo di definizione (esempio [0,2])
a=0; b=2;

% punto in cui valutare la derivata terza
x1=1;

% passo
h=10^(-3);

% calcolo derivata terza (numericamente)
y1h = derivata_terza (f,x1,h,a,b);

H=h/10;
y1H = derivata_terza (f,x1,H,a,b);

err_est=abs(y1h-y1H);

% valutazione della derivata terza (analiticamente)
y1=feval(D3f,x1);

% errore assoluto
ae=abs(y1-y1h);
re=abs(y1-y1h)/abs(y1);

% stampa
fprintf('\n \t y1h     : %1.5f',y1h);
fprintf('\n \t y1      : %1.5f',y1);
fprintf('\n \t ae      : %1.3e',ae);
fprintf('\n \t re      : %1.3e',re);
fprintf('\n \t err_est : %1.3e',err_est);
fprintf('\n');

% stampa su file
fid=fopen('esperimento.txt','w');
fprintf(fid,'\n \t y1h     : %1.5f',y1h);
fprintf(fid,'\n \t y1      : %1.5f',y1);
fclose(fid);


% RISULTATI:
%
% >> demo_derivata_terza
% 
%   y1h     : -3.28036
%   y1      : -3.28172
%   ae      : 1.360e-03
%   re      : 4.145e-04
%   err_est : 1.461e-03
% 
% >>
