import numpy as np;
import matplotlib. pyplot as plt;
#生成模拟数据
np.random.seed(2910); #可以设定种子，多次执行有不同的结果。
mu, sigma = 100, 15; #按正态分布的均值100，标准差15
x = mu + sigma * np.random.randn(10000);
#直方图
n, bins, patches = plt.hist(x, 50, density=True, facecolor='r', alpha=0.75); #r为红色，b为蓝色，g为绿色。
plt. xlabel('XXX'); #可以修改标签。
plt. ylabel('Probability');
plt. title('Histogram of ORSCI');
plt. text(60,.025,r'$\mu=100,\ \sigma=15$');
plt. axis([40,160,0,0.03]);
plt. grid(True);
plt. show();

略