function [mmom,digmm]=dig_mmEnufin(N,nu,c,dig0,dd,nofdig)
%   DIG_MMENUFIN determines the number digmm of digits  
%   needed to obtain to an accuracy of nofdig digits the
%   2*Nx1 array mmom of modified moments relative to monic
%   Jacobi01 polynomials on [0,c] of the subrange exponential
%   integral weight function on [0,c].

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