Function numlib::integrate::three_eighths_simpson

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