Function simpson 

pub fn simpson<F: Fn(f64) -> f64>(f: &F, a: f64, b: f64, panels: usize) -> f64

Composite Simpson’s rule for a scalar integrand over [a, b] with panels subintervals (forced even, ≥ 2). The general-N version of the fixed 100-chunk solver in super.