pub fn hessian<const N: usize, F>( function: F, point: Vector<N>, ) -> Result<LinearMap<N, N>, AnalysisError>where F: Fn(Vector<N>) -> f64,