I have generated a distribution (normalized such that the sum is equal to 1) by using the code:
M=500; % Number of samples
z=1;
SUM = 1;
ns = rand(1,M).^z; % random numbers
TOT = sum(ns);
X = (ns/TOT)*SUM; % Re-scaling
hist(X(1,:),100)
For an exponent ##z=1##, the sampling distribution is...