function X=moment(pf, n) 
% moment(pf, n)  calculates nth moment of a pdf or pmf 
if ispf(pf)==0 
   error('invalid probability function') 
end 
x=xvals(pf); 
p=pvals(pf); 
g=(x.^n).*p; 
if isdiscrt(pf)==1 
   X=sum(g); 
else 
   X=trapz(x,g); 
end