Function numerical_hessian 

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

Finite-difference Hessian matrix (central differences, second order).