Function integrate 

pub fn integrate<F>(f: F, a: f64, b: f64) -> f64
where
    F: Fn(f64) -> f64,

Integrates a function from a to b. Uses the Tanh-Sinh quadrature over [-1, +1] and then transforms to an integral over [a, b].