庞加莱买面包问题

import numpy as np
import scipy.stats as sta

import matplotlib.pyplot as plt

X = sta.norm(loc=950, scale=20)
# 面包的重量列表
wbread = []

for i in range(365):
    x = X.rvs(size=100)
    # 从100 个面包中随便拿一个，我们就拿第 1 个
    wbread.append(x[0])

print(np.mean(wbread))
# 偏度: -0.088739676238455
print('偏度:', sta.skew(wbread))
plt.subplot(1, 2, 1)
plt.hist(wbread)

# 第 2 年，我们每次拿出最大

wbread = []
for i in range(365):
    x = X.rvs(size=100)
    # 从100 个面包中随便拿一个，我们就拿第 1 个
    wbread.append(max(x))

print(np.mean(wbread))
# 偏度: 0.5133333290340636，
# 功能判断出老板的行为，每次都给他最大的面包
print('偏度:', sta.skew(wbread))
plt.subplot(1, 2, 2)

plt.hist(wbread)
plt.show()