pub fn gamma_pdf(x: Matrix<f64>, k: f64, theta: f64) -> Matrix<f64>
PDF of the Gamma distribution for matrices