rand_props2.m

function [prop_mat] = rand_props(num_prey,num_vecs)

rand_nums=sort(rand(num_prey-1,num_vecs),1); %sets up and sorts mat of random nums
prop_mat=zeros(num_prey,num_vecs);
prop_mat(1,:)=rand_nums(1,:);

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