Function regular 

fn regular(v: f64, order: i32, y: f64) -> Result<(Scaled, Scaled)>

(P_ν^{-μ}(y), (1 - y²) dP_ν^{-μ}/dy) for 0 <= y < 1.