pub fn pascal_pdf(k: u32, p: f64, n: u32) -> f64
This function computes the probability p(k) of obtaining k from a Pascal distribution with parameters p and n, using the formula given above.