pub fn damped_newton<const N: usize, F>( function: F, initial: Vector<N>, tolerance: f64, max_iterations: u32, ) -> Result<NewtonReport<N>, AnalysisError>where F: Fn(Vector<N>) -> Vector<N>,