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