import numpy as np draws = np.random.default_rng(41).binomial(n=8, p=0.75, size=50) print(draws.shape) print(bool(((draws >= 0) & (draws <= 8)).all())) print(np.random.default_rng(41).binomial(n=8, p=1., size=3).tolist())