Compute finite difference gradients

compute_finite_difference_gradients(f, X, h)