Function jacobian 

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