Function binomial_sample 

pub fn binomial_sample(
    n_trials: usize,
    p: f64,
    n: usize,
    rng: &mut Rng,
) -> Vec<i64>

Sample n values from Binomial(n_trials, p).

Uses direct simulation for small n_trials (< 25) and normal approximation for larger values.