%SMOMENU
%
function mom=smomEnu(dig,N,nu)
digits(dig); 
mom=vpa(zeros(2*N,1)); snu=vpa(nu);
mom(1)=1/snu;
for k=1:2*N-1
  mom(k+1)=(k*(snu+k-1)/(snu+k))*mom(k);
end
