function ab=sr_Enufinmod(dig,N,nu,c)
%SR_ENUFINMOD Recurrence coefficients for a subrange 
%exponential integral weight function on [0,c].
%   AB=SR_ENUFINMOD(DIG,N,NU,C) generates for selected values
%   of nu and c the Nx2 array of the first N recurrence 
%   coefficients in DIG-digit arithmetic for the weight 
%   function E_NU(x) on [0,C], using the modified Chebyshev
%   algorithm with modified moments relative to shifted Jacobi
%   polynomials with Jacobi parameters alpha=0, beta=NU-1,
%   transformed from the interval [0,1] to the interval [0,C].

disp('Be sure to load the correct file mmom_Enufin_nu_c.txt on line 15')
digits(dig); nu=vpa(nu); c=vpa(c);
a=0; b=nu-1;
mmom=loadvpa('mmom_Enufin_1_20.txt',9,50,80,1); mmom=vpa(mmom,dig);
abmm0=sr_jacobi01(dig,2*N-1,a,b);
abmm(:,1)=c*abmm0(:,1); abmm(1,2)=(c^nu)*abmm0(1,2);
abmm(2:2*N-1,2)=c^2*abmm0(2:2*N-1,2);
ab=schebyshev(dig,N,mmom',abmm);
