opensrdk_optimization::numerical_diff
pub fn numerical_diff( func: &dyn Fn(&[f64]) -> f64, x: &[f64], h: f64) -> Vec<f64>