
function [fp_l, fp_r, fp_c]=derivative_point(x0, deltax, f)

fp_l=(f(x0+deltax)-f(x0))./deltax;
fp_r=(f(x0)-f(x0-deltax))./deltax;
fp_c=(f(x0+deltax)-f(x0-deltax))./(2*deltax);

end