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);