function mmom=smmom_Enu(dig,N,nu)
%   MMOM=SMMOM_ENU(DIG,N,NU) computes in DIG-digit
%   arithmetic the first 2N modified moments, relative
%   to monic generalized Laguerre polynomials with 
%   parameter a=-1, of the exponential integral weight
%   function w(x)=E_nu(x) on [0,Inf].   

tic
digits(dig); sm1=vpa(-1); s0=vpa(0); s1=vpa(1); snu=vpa(nu);
mmom=vpa(zeros(2*N,1));
ablm1=sr_laguerre(dig,2*N-1,sm1); 
abl=sr_laguerre(dig,N); xwl=sgauss(dig,N,abl);
abj=sr_jacobi01(dig,N,s0,snu-s1); xwj=sgauss(dig,N,abj);
for k=1:2*N
  u1=0;
  for i=1:N
    u2=s0;
    for j=1:N
      t=xwj(j,2)*psubk(dig,xwl(i,1)*xwj(j,1),k-1,ablm1);
      u2=u2+t;
    end
    t=xwl(i,2)*u2;
    u1=u1+t;
  end
  mmom(k)=u1;
end
toc
