import numpy as np
resultsSum = 0
days = 30
lam = 2 # how many audits in a day
iterations = 100
for i in range(iterations):
s = np.random.poisson(lam=(lam), size=(days))
print("\n")
resultsSum += sum(s)
print(resultsSum / iterations) # how many audits in a month on average (should be close to days*lam)