function ab=sr_explog_fin(nofdig,N)
%SR_EXPLOG_FIN Recurrence coefficients for a weight function on
%[0,1] containing algebraic, exponential, and extended-logarithmic 
%factors.
%   AB=SR_EXPLOG_FIN(NOFDIG,N) generates in DIG-digit arithmetic
%   the first N recurrence coefficients for the weight function
%   w(x)=x^A*exp(-x)*log(1+1/x) on [0,1], A>=0 an integer, using a 
%   1-component discretization procedure, and places them to NOFDIG 
%   digits into the Nx2 array AB. The variables DIG and A must be 
%   netered when prompted by the routine. 

global mc mp iq idelta irout dig a abexplog01
dig=input('enter the value of dig');
a=input('enter the value of a');
%
% Be sure to insert the numerical value of a into the "alternative
% command" Mcap=ceil((2*N+a)/2-1) of the routine smcdis.m. and check
% that Mcap <= 65 for the chosen values of N and a.
%
abexplog01=loadvpa('coeff_explog0_fin.txt',9,32,65,2);
Mmax=100; eps0=1;
mc=1; mp=0; iq=1; idelta=2; irout=1; 
[ab,Mcap,kount]=smcdis(dig,N,eps0,@squad_explog_fin,Mmax); 
ab=vpa(ab,nofdig);
[Mcap kount]
