pub fn geometric_P(k: u32, p: f64) -> f64
This function computes the cumulative distribution functions P(k), Q(k) for the geometric distribution with parameter p.