Function polymers::math::integrate_1d

pub fn integrate_1d(
    f: &dyn Fn(&f64) -> f64,
    x_min: &f64,
    x_max: &f64,
    num_points: &u128,
) -> f64