Function numerical_gradient 

pub fn numerical_gradient(
    f: impl Fn(&[f64]) -> f64,
    x: &[f64],
    h: f64,
) -> Vec<f64>

Central-difference numerical gradient.