#Schemat losowania Hartleya-Rao
losHartleyRao<-function(n,N,pi1)
{skum<-matrix(0,ncol=1,nrow=N+1)
skum[2]<-pi1[1]
for (k in 2:N) skum[k+1]<-skum[k]+pi1[k]
u<-runif(1)
sn<-matrix(0,ncol=1,nrow=n)
z<-1
while (z<=n)
{
k<-1
q<-0
while (q==0)
{
if ((skum[k]<u+z-1)&(u+z-1<=skum[k+1])) {q<-1; sn[z]<-k}
k<-k+1
}
z<-z+1
};sn}