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