%
%SR_ENUMOD Recurrence coefficients for an exponential integral
%weight function 
%   AB=SR_ENUMOD(DIG,N,NU,NOFDIG) generates in DIG-digit
%   arithmetic the first N recurrence coefficients for the
%   weight function w(x)=E_NU(x) on (0,Inf) from its 2N modified
%   moments relative to a monic generalized Laguerre polynomial
%   with Laguerre parameter -1 and places them to NOFDIG digits
%   into the Nx2 array AB. 

function ab=sr_Enumod(dig,N,nu,nofdig)
tic
digits(dig); m1=vpa(-1);
mmom=loadvpa('mmom_Enu_1half.txt',9,50,80,1);
%mmom=loadvpa('mmom_Enu_1.txt',9,50,80,1);
%mmom=loadvpa('mmom_Enu_2.txt',9,50,80,1);
%mmom=loadvpa('mmom_Enu_5.txt',9,50,80,1);
abm=sr_laguerre(dig,2*N-1,m1);
ab=schebyshev(dig,N,mmom',abm);
ab=vpa(ab,nofdig);
toc