function pdf=mfberni(p) 
% mfberni(p)  generates a Bernoulli pmf with parameter p 
if(p<0 | p>1) 
   error('invalid value for p') 
end 
pdf=[1 1;0 1-p;1 p];