import numpy.random as rd
import numpy as np

# Q1
f=np.zeros(100000)
for k in range (0,10000): 
    n = rd.poisson(40)
    x = rd.binomial (n,1/10)
    f[k] = x
print(f)

# Q2
m = np.min(f)
M = np.max(f)
print(m,M)