function ab=r_Enumod(dig,N,nu,a)
%R_ENUMOD Recurrence coefficients for the exponential integral
%weight function 
%   R_ENUMOD(N,NU,A) generates in double-precision arithmetic
%   the first N recurrence coefficients for the weight function
%   E_NU(X) on [0,INF] from its modified moments relative to
%   generalized Laguerre polynomials with parameter A and places
%   them into the Nx2 array AB.

dig=32;
smmom=smmom_Enu(dig,N,nu,a); 
mmom=double(smmom),
abm=r_laguerre(2*N-1,a);
ab=chebyshev(N,mmom',abm);
