Python-сообщество

 import scipy.stats as st
import pandas as pd
import matplotlib.pylab as plt
import numpy as np
M = 0.9
sigma = 0.3
norm = st.norm(loc=M, scale=sigma) # making a random variable
sample = norm.rvs(size=10000) # making a sample
real_M = np.mean(sample)
real_sigma = np.std(sample)
hist_bins = [real_M-3*real_sigma, real_M-2*real_sigma, real_M-real_sigma,
             real_M, real_M+real_sigma, real_M+2*real_sigma, real_M+3*real_sigma]
hist_data = np.histogram(sample, bins=hist_bins)[0]
hist_data = [int(x*100)/len(sample) for x in hist_data]
%matplotlib inline
x = np.arange(M-4*sigma, M+4*sigma, 0.01)
y = norm.pdf(x)
plt.plot(x,y)
plt.hist(sample, bins=20, normed=True)
plt.gcf().set_size_inches(6.5, 6.5, True)
plt.show()

 plt.hist(sample, bins=20, normed=True, edgecolor='yellow')