function s=first_derivative(a,h)
s=(f(a+h)-f(a-h))/(2*h);
s