Function adaptive_simpson 

pub fn adaptive_simpson<F: Fn(f64) -> f64>(
    f: &F,
    a: f64,
    b: f64,
    tol: f64,
    depth: usize,
) -> f64

Adaptive Simpson’s rule.