% # of sample pointsn_samples = 1000;p = ones(3,1)/3;% Low Entropyalpha = 0.5*p;% High Entropy%alpha = 10*p;% draw n_samples random points from the 3-d dirichlet(alpha),% and plot the resultspoints = zeros(3,n_samples);for i = 1:n_samples    points(:,i) = dirichletrnd(alpha);endscatter3(points(1,:)', points(2,:)', points(3,:)', 'r', '.', 'filled');