pub fn pascal_P(k: u32, p: f64, n: u32) -> f64
This function computes the cumulative distribution functions P(k), Q(k) for the Pascal distribution with parameters p and n.