Function rgsl::randist::dirichlet::dirichlet_pdf

pub fn dirichlet_pdf(alpha: &[f64], theta: &[f64]) -> f64

This function computes the probability density p(\theta_1, … , \theta_K) at theta[K] for a Dirichlet distribution with parameters alpha[K], using the formula given above.