import numpy as np
from scipy import stats as st

mu = 4.5*2

pn = 1

for n in range(0,12):
	pn = pn - st.poisson.pmf(n, mu=mu)


print('La probabilidad es:', pn)