pub fn bernoulli_pdf(x: u32, p: f64) -> f64
This function computes the probability p(k) of obtaining k from a Bernoulli distribution with probability parameter p, using the formula given above.