pub fn complex_step_derivative<F>( function: F, x: f64, ) -> Result<DerivativeEstimate, AnalysisError>where F: Fn(Complex64) -> Complex64,