pub fn laplace_pdf(x: f64, a: f64) -> f64
This function computes the probability density p(x) at x for a Laplace distribution with width a, using the formula given above.