function ab=r_Enufinmod(N,nu,c)
%R_ENUFINMOD Recurrence coefficients for a subrange 
%exponential integral weight function on [0,c].
%   AB=R_ENUFINMOD(N,NU,C) for selected values of NU and C 
%   generates in double-precision arithmetic the Nx2 array 
%   of the first N recurrence coefficients for the weight  
%   function w(x)=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 in')
disp('line 14 of this file')
mmom=loadvpa('mmom_Enufin_1_20.txt',9,50,80,1);
%mom=vpa(mmom,15);
mom=double(mmom);
abmm0=r_jacobi01(2*N-1,0,nu-1);
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=chebyshev(N,mom',abmm); ab=double(ab);
