function [mom,digm]=dig_mEnufin(N,nu,c,dig0,dd,nofdig)
%   DIG_MENUFIN determines the number digm of digits  
%   needed to obtain to an accuracy of nofdig digits the
%   2Nx1 array mom of moments of the subrange exponential
%   integral weight function on [0,c].

i=dig0-dd; 
maxerrm=1;
while maxerrm>.5*10^(-nofdig) 
  i=i+dd; digm=i; digm 
  if digm==dig0
    m0=smom_Enufin(digm,N,nu,c);
  else
    m1=smom_Enufin(digm+dd,N,nu,c);
    serrm=abs((m1-m0)./m1);
    errm=double(serrm);
    maxerrm=max(max(errm)) 
    m0=m1;
  end
end
mom=vpa(m1,nofdig); 
