Function romberg 

pub fn romberg(
    f: impl Fn(f64) -> f64,
    a: f64,
    b: f64,
    max_levels: usize,
    tol: f64,
) -> f64

Romberg integration — Richardson extrapolation on the trapezoidal rule.