rand_props2.m

function [prop_vec] = rand_props(num_prey)

rand_nums=sort(rand(num_prey-1,1)); %sets up vec for sorting
prop_vec=zeros(num_prey,1);
prop_vec(1,1)=rand_nums(1,1);

for j=2:num_prey-1
   prop_vec(j,1)=rand_nums(j,1)-rand_nums(j-1,1);
end
prop_vec(num_prey,1)=1-sum(prop_vec);