# coding: utf-8
import numpy as np
import pylab as pl

data = np.random.normal(0, 1, 100000)
mu = np.mean(data)
sigma = np.std(data)

pl.figure()
pl.subplot(2,1,1)
pl.xlabel('Score', size=14)
pl.ylabel('Frequency', size=14)
pl.ylim(0,10000)
pl.text(-4, 9500, '''
        $\mu=%.1f$
        $\sigma=%.2f$''' % (mu, sigma),
        verticalalignment='top', horizontalalignment='left')
pl.grid(True)
pl.axvline(x=mu, linewidth=1, color='r')
#pl.hist(abc, bins=40, range=(-5, 5), normed=True)
pl.hist(data, bins=40, range=(-4, 4),facecolor='cyan')

pl.subplot(2,1,2)
pl.xlabel('Score', size=14)
pl.ylabel('Cumulative probability (%)', size=14)
pl.ylim(0.0,1.0)
pl.grid(True)
pl.hist(data, bins = 10000, range=(-4, 4), normed=True, cumulative=True, histtype='step')

pl.tight_layout()
#pl.show()
pl.savefig("fig_hist.png", dpi=100)