數值微分 三點

function numerical_diff = ND3(f,x,h)

f_1=f(x);
f_2=f(x+h);
f_3=f(x+2*h);


numerical_diff = ((-3/2)*f_1+2*f_2-(1/2)*f_3)./h;

end