function mom=smom_Enu(dig,N,nu)
%SMOM_ENU Moments of the exponential integral weight function.
%   MOM=SMOM_ENU(DIG,N,NU) generates in DIG-digit arithmetic
%   the (2N)x1 array MOM of the first 2N moments of the weight
%   function W(X)=E_NU(X) on (0,INF], nu>0.

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
