return np.array([(f(x + hi) + f(x - hi)) / 2.0 - f_x for hi in steps])