
function demo

f=@(x) sin(x); f1=@(x) cos(x);
f1x_esatta=feval(f1,pi/4);
for k=1:6
    h(k)=10^(-k);
    hL=h(k);
    f1x_da=differenza_in_avanti(f,pi/4,hL);
    err(k)=abs(f1x_esatta-f1x_da);
    fprintf('\n \t %1.5e',err(k));
end

semilogy(h,err,'bo');

FID=fopen('risultati.txt','w');
fprintf(FID,'\n %1.5e',err);
ST=fclose(FID);

fprintf('\n \n')

% >>demo
% 
%  	 3.65038e-02
%  	 3.54729e-03
%  	 3.53671e-04
%  	 3.53565e-05
%  	 3.53554e-06
%  	 3.53442e-07
%
% >>




function f1x=differenza_in_avanti(f,x,h)

f1x=(f(x+h)-f(x))/h; 

% NB: si puo' mettere il "." pensando "x" e "h" come vettori di uguale 
% dimensione.