pub fn sph_reg(l: i32, lambda: f64, x: f64) -> f64
This routine computes the Regular Spherical Conical Function P^{-1/2-l}_{-1/2 + i \lambda}(x) for x > -1, l >= -1.