import pandas as pd
import numpy as np
#import scipy.stats as stats 
from scipy.stats import poisson

r = 5

poisson_02 = np.empty(r)
poisson_04 = np.empty(r)
poisson_06 = np.empty(r)
poisson_08 = np.empty(r)
poisson_1 = np.empty(r)
poisson_2 = np.empty(r)
poisson_3 = np.empty(r)
poisson_4 = np.empty(r)

# repeat r times:
for j in  range(r):
  poisson_02[j] = poisson.cdf(k=j, mu=0.02)
  j=j+1
  
for j in  range(r):
  poisson_04[j] = poisson.cdf(k=j, mu=0.04)
  j=j+1

for j in  range(r):
  poisson_06[j] = poisson.cdf(k=j, mu=0.06)
  j=j+1

for j in  range(r):
  poisson_08[j] = poisson.cdf(k=j, mu=0.08)
  j=j+1

for j in  range(r):
  poisson_1[j] = poisson.cdf(k=j, mu=0.1)
  j=j+1

for j in  range(r):
  poisson_2[j] = poisson.cdf(k=j, mu=0.2)
  j=j+1

for j in  range(r):
  poisson_3[j] = poisson.cdf(k=j, mu=0.3)
  j=j+1

for j in  range(r):
  poisson_4[j] = poisson.cdf(k=j, mu=0.4)
  j=j+1

poisson_mu = pd.DataFrame({'0.02':poisson_02,'0.04':poisson_04, '0.06':poisson_06,'0.08':poisson_08, '0.1':poisson_1,'0.2':poisson_2,'0.3':poisson_3,'0.4':poisson_4})
round(poisson_mu, 3)
print("Probability Distribution Table of Poisson with 'mu=0.02, ..., 0.4' : ", f'\n{round(poisson_mu, 3)}\n')