inkluzje<-function(n,dx){
#wyznaczanie prawdopodobieńsw inkluzji nie wiekszych od 1
n1<-n
dx1<-dx
N<-nrow(dx)
pi1<-matrix(0,ncol=1,nrow=N)
sp<-sum(dx)
pi1<-dx*(n/sp)
pi1
km<-1
m<-pi1[1]
for (k in 2:N) if (pi1[k]>m) {m<-pi1[k]; km<-k}
if (m>1) {pi1[km]<-1;dx1[km]<-0}
while (m>1)
{n1<-n1-1
sp<-0
for (k in 1:N) sp<-sp+dx1[k]
for (k in 1:N) if (dx1[k]>0) pi1[k]<-dx1[k]*(n1/sp)
m<-pi1[1]
for (k in 2:N) if (pi1[k]>m) {m<-pi1[k]; km<-k}
if (m>1) {pi1[km]<-1;dx1[km]<-0}
}
pi1
}