# xi**3 empírikus és elmélet eloszlásfv-e

#import scipy
from scipy import linspace,cumsum
from scipy.stats import uniform
import matplotlib.pyplot as plt

# először szimulációval:
N=1000
xi3=uniform.rvs(0,1,N)**3
xi3.sort()
y=linspace(0,1-1/N,N)
plt.step(xi3,y) # lépcsős stílus
#plt.legend('szimulált')


xx=linspace(0,1)
yy=xx**(1.0/3)
plt.plot(xx,yy)

plt.legend(['szimulált','pontos'])
plt.title('xi**3 eloszlása (xi egyenletes [0,1]-en)')
plt.show()