Function numlib::integrate::adaptive_simpson

pub fn adaptive_simpson(a: f64, b: f64, f: fn(_: f64) -> f64, eps: f64) -> f64