pub fn beta_pdf(x: f64, a: f64, b: f64) -> f64
This function computes the probability density p(x) at x for a beta distribution with parameters a and b, using the formula given above.