pub fn legendre_twofloat(n: usize) -> Rule<Compensated<f64, f64>>
Create a Gauss-Legendre quadrature rule with n points on [-1, 1] (Df64 version).