Function hessian 

pub fn hessian<F>(f: F, x: &[f64], h: f64) -> Vec<Vec<f64>>
where
    F: Fn(&[f64]) -> f64,

Hessian via central differences (n×n matrix for n-dim function).