function [mom,dig]=dig_mmEnu(N,nu,a,dig0,dd,nofdig)
%DIG_MMENU Number of digits needed for the first 2*N modified
%moments, based on generalized Laguerre polynomials with 
%parameter A, of the exponential integral weight function 
%E_nu(x) on [0,Inf].

i=dig0-dd; 
maxerr=1;
while maxerr>.5*10^(-nofdig)
  i=i+dd; dig=i
  if i==dig0
    mom0=smmom_Enu(dig,N,nu,a);
  else
    mom1=smmom_Enu(dig,N,nu,a);
    serr=abs((mom1-mom0)./mom1);
    err=double(serr);
    maxerr=max(max(err)) 
    mom0=mom1;
  end
end
mom=vpa(mom1,nofdig);
