pub fn adaptive_integration( f: &dyn Fn(f64) -> f64, a: f64, b: f64, epsabs: f64, epsrel: f64, ) -> (f64, f64)