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