import scipy.stats as ss
import numpy as np
import matplotlib.pyplot as plt
fig = plt.figure(figsize=(8, 8))
ax = fig.add_subplot(111)
alpha, loc, beta = 5, 100, 22
data = ss.gamma.rvs(alpha, loc=loc, scale=beta, size=5000)
myHist = ax.hist(data, 100, normed=True)
rv = ss.gamma(alpha, loc, beta)
x = np.linspace(0, 600)
h = ax.plot(x, rv.pdf(x), lw=2)
plt.show()
python: построение гистограммы с функциональной линией сверху