function y=derivative1(f,x,h)
y=(f(x+h)-f(x))/h;
