import numpy as np
import matplotlib.mlab as mlab
import matplotlib.pyplot as plt
from scipy.stats import norm
# 产生数据
mu =100# mean of distribution
sigma =15# standard deviation of distribution
x = mu + sigma * np.random.randn(10000)
num_bins =50# 区间个数
n, bins