pub fn binomial_pmf(n: u64, k: Matrix<u64>, p: f64) -> Matrix<f64>
PMF of the Binomial(n, p) distribution for matrices