peroxide::special::legendre
pub fn unit_gauss_legendre_quadrature<F>(f: F, n: usize) -> f64 where F: Fn(f64) -> f64,