function fp=diff_1_2(f, x, h)   
    x_ph = x + h;
    x_mh = x - h;
    fp=(f(x_ph)-f(x_mh))./(2*h);
end