return (f(x0i + h) - f(x0i - h)) / 2.0