//
// studio dell'errore nel calcolo della derivata prima
// con le differenze in avanti 
//

h = [ 0.1 0.01 0.001 1d-4 1e-7 1e-12]';

// approssimazione della derivata di exp(1) differenze in avanti


df = (exp(1+h)-exp(1))./h;

errass = abs(df-exp(1));

[h errass]