Function numlib::integrate::simpson_quadrature

pub fn simpson_quadrature(a: f64, b: f64, n: i64, f: fn(_: f64) -> f64) -> f64