import numpy as np
import scipy.stats as stPoisson = st.poissonl = 1.5 / 100 * 48X = Poisson(l)
X.pmf(0)0.4867522559599717np.sum(X.pmf(np.arange(6)))0.9998951831177709Y = Poisson(1.5 / 100 * 24)
1 - np.sum(Y.pmf(np.arange(3)))0.00595077061399496