function pmf=mfbino(n, p) 
% mfbino(n, p)  generates a binomial pmf with parameters n, p 
if p<0 | p>1 
   error('invalid value for p') 
end 
if n<=0 
   error('invalid value for n') 
end 
newx=[0:n]; 
n=ones(1,length(newx))*n; 
newp=cmb(n, newx).*(p.^newx).*((1-p).^(n-newx)); 
pmf=[1 1; newx' newp'];