import numpy as n
import matplotlib.pyplot as pl

a=5
b=25
data=n.random.randint(a,b,100000)
y,binEdges = n.histogram(data,bins=10)
bincenters=0.5*(binEdges[1:]+binEdges[:-1])
err=n.sqrt(y)
# szerokosc slupka w histogramie
w=0.9*(b-a)/10.
# Zakres poziomej osi
pl.xlim(5,26)
# przygotowanie histogramu z bledami
pl.bar(bincenters,y,color='g',width=w,yerr=err)
pl.title('liczby losowe wygenerowane z plaskiego rozkladu')
pl.show()

h1 = n.random.poisson(5, 100000)
# przygotowanie histogramu z rozkładu Poissona
count1, bins1, ignored1 = pl.hist(h1, 14, density=True, color='r', label='Rozklad Poissona')
# losowanie z rozkładu trójkątnego
h2 = n.random.triangular(15, 20, 25, 100000)
# przygotowanie histogramu z rozkładu trójkątnego
count2, bins2, ignored2 = pl.hist(h2, 50, density=True, color='b', label='Rozklad trojkatny')
pl.ylabel('Gestosc prawdopodobienstwa')
pl.xlabel('argument')
pl.legend()
pl.show()