
function demo_derivata

% funzione
f=@(x) sin(x)-cos(x);

% derivata: D(sin(x))=cos(x), D(cos(x))=-sin(x)
Df=@(x) cos(x)+sin(x);

% punti
x=0:0.01:1;

% intervallo di definizione (comprende [0,1])
a=-1; b=2;

% passo
h=10^(-3);

% lunghezza "x"
L=length(x);

% calcolo derivate seconde (numericamente)
for k=1:L
    xk=x(k);
    [y1(k), err_est(k)] = rapp_incr_simm (f,xk,h,a,b);
end

% valutazione delle derivate seconde "esatte"
Dfx=feval(Df,x);

% massimo errore
ae=abs(Dfx-y1);
max_ae=max(ae);

% stima del massimo errore
max_est=max(err_est);

% stampa
fprintf('\n \t max_ae: %1.3e max_est: %1.3e \n',max_ae,max_est);

% salva su file
fid=fopen('errori.txt','w'); 
fprintf(fid,'\n \t max_ae : %1.15e',max_ae); 
fprintf(fid,'\n \t max_est: %1.15e',max_est); 
fclose(fid); 


% RISULTATI:
% >> demo_derivata
%    max_ae: 2.357e-07 max_est: 2.357e-07 
% >> 
%
% SU FILE:
%  	 max_ae : 2.356998440866676e-07
%  	 max_est: 2.356700388173039e-07
