function pmf=mfuni(range) 
% mfuni(range)  generates a uniform pmf over the given range row vector 
R=length(range); 
p=ones(R, 1)./R; 
pmf=[1 1; range' p];