function sigma=var(pf) 
% var(pf)  variance of a pdf or pmf 
if ispf(pf)==0 
   error('invalid probability function') 
end 
sigma=central(pf, 2);