function h_eff = testder2(h)
% codice per confrontare il 
% grafico della funzione cos(x)
% sull'intervallo [0,2pi]
% con il valore approsimato numericamente 
% della derivata della funzione sin(x)
% mediante la formula del secondo ordine
% f'(x0) = (f(x0+h) - f(x0-h))/2h

%  

n = ceil(2*pi/h);

x = linspace(0, 2*pi, n);

h_eff = x(2) - x(1);

f = sin(x);

%for i = 2:(end-1)
%    fprime(i) = (f(i+1) - f(i-1))/(2*h);
%end

fplus  = f(3:end);
fminus = f(1:end-2);

fprime = (fplus-fminus)./(2*h_eff);

figure

plot(x, cos(x), 'g', x(2:end-1), fprime, 'r');




