function pmf=mfhyper(m, n, k) 
% mfhyper(m, n, k)  generates a hypergeometric pmf with parameters m, n, k 
% 
%    Like lotto: 
%       m  number of "good" numbers 
%       n  number of "bad" numbers 
%       k  how many numbers you pick 
if m<0 
   error('invalid value for m') 
end 
if n<0 
   error('invalid value for n') 
end 
if k<0 
   error('invalid value for k') 
end 
x=[0:k]; 
m=ones(size(x))*m; 
n=ones(size(x))*n; 
k=ones(size(x))*k; 
p=hypergeo(m+n, m, k, x); 
pmf=[1 1; x' p'];