01_preparation.py

%matplotlib inline 
import matplotlib.pyplot as plt
import numpy as np
import scipy.stats as st
from matplotlib import animation as ani
plt.style.use('ggplot')
plt.rc('text', usetex=True)

02_gen_random_poisson.py

x = st.poisson.rvs(5, size=30)   # 平均=5のポアソン分布に従う乱数を30個生成
print x

03_draw_histgram.py

lam = 3     # 平均値
n = 10000   # 生成する乱数の数
x = st.poisson.rvs(5, size=n)

plt.figure(figsize=(9,7))
plt.title("Poisson distribution.  $\lambda=5$, n={}".format(n))
plt.hist(x, bins=np.max(x), color="lightblue")
plt.show()