pub fn legendre_Plm_e(l: i32, m: i32, x: f64) -> (Value, Result)
This routine computes the associated Legendre polynomial P_l^m(x) for m >= 0, l >= m, |x| <= 1.