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