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