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

r = 5

binom_01 = np.empty(r)
binom_02 = np.empty(r)
binom_03 = np.empty(r)
binom_04 = np.empty(r)
binom_05 = np.empty(r)
binom_06 = np.empty(r)
binom_07 = np.empty(r)
binom_08 = np.empty(r)
binom_09 = np.empty(r)
binom_095 = np.empty(r)

# repeat r times:
for j in  range(r):
  binom_01[j] = binom.cdf(k=j, n=5, p=0.1)
  j=j+1
  
for j in  range(r):
  binom_02[j] = binom.cdf(k=j, n=5, p=0.2)
  j=j+1

for j in  range(r):
  binom_03[j] = binom.cdf(k=j, n=5, p=0.3)
  j=j+1

for j in  range(r):
  binom_04[j] = binom.cdf(k=j, n=5, p=0.4)
  j=j+1

for j in  range(r):
  binom_05[j] = binom.cdf(k=j, n=5, p=0.5)
  j=j+1

for j in  range(r):
  binom_06[j] = binom.cdf(k=j, n=5, p=0.6)
  j=j+1

for j in  range(r):
  binom_07[j] = binom.cdf(k=j, n=5, p=0.7)
  j=j+1

for j in  range(r):
  binom_08[j] = binom.cdf(k=j, n=5, p=0.8)
  j=j+1

for j in  range(r):
  binom_09[j] = binom.cdf(k=j, n=5, p=0.9)
  j=j+1

for j in  range(r):
  binom_095[j] = binom.cdf(k=j, n=5, p=0.95)
  j=j+1

binom_n5 = pd.DataFrame({'0.10':binom_01,'0.20':binom_02,'0.30':binom_03,'0.40':binom_04,'0.50':binom_05,'0.60':binom_06,'0.70':binom_07,'0.80':binom_08,'0.90':binom_09,'0.95':binom_095 })
round(binom_n5,3)
print("Probability Distribution Table of Binomial with 'n=5 and p=0.10, ..., 0.95' : ", f'\n{round(binom_n5,3)}\n')