plt.hist(samples, range=(0, samples.max()), bins=samples.max() + 1, normed=
    True, histtype='stepfilled')
Is there a bug in binning in matplotlib histograms? Or non-randomness of the rvs method in scipy.stats