
% NOME:    COGNOME:    MATRICOLA:    POSTAZIONE:   

function f1x = differenza_centrale (f,x,h)

% Oggetto:
% Differenza centrale
% Input:
% f: funzione di cui si vuole approssimare la derivata nei punti "x";
% x: punti in cui si vuole approssimare la derivata (x puo' essere un vettore);
% h: passo (e' un numero);
% Output:
% f1x: approssimazione delle derivate nel punti "x".

% test su "h"
if h == 0, f1x=NaN*ones(size(x)); end

% approssimazione della derivata
f1x=(-feval(f,x+2*h)+8*feval(f,x+h)-8*feval(f,x-h)+feval(f,x-2*h))/(12*h);

