function Msk = maskv(N,p)

inds = ones(1,p);
weight = 1;

nP = factorial(N+p-1)/factorial(N-1)/factorial(p);
Msk = zeros(int64(nP),1);


for n=1:nP
    
    Msk(n) = weight;
    
    [weight,inds] = nxtseqdf(inds,N);
    
end    
