pub fn legendre(n: u32, x: f64) -> f64
Legendre polynomial P_n(x). (n+1)P_{n+1} = (2n+1)x P_n − n P_{n-1}.
P_n(x)
(n+1)P_{n+1} = (2n+1)x P_n − n P_{n-1}