pub fn simpsons(f: impl Fn(f64) -> f64, a: f64, b: f64, n: usize) -> f64
Simpson’s 1/3 rule. n must be even; if odd, n is incremented by 1.