function ab=sr_explog(nofdig,N)
%SR_EXPLOG Recurrence coefficients for a weight function on
%[0,Inf] with an algebraic/extended-logarithmic singularity.
%   AB=SR_EXPLOG(NOFDIG,N) generates in DIG-digit arithmetic
%   the Nx2 array AB of the first N recurrence coefficients for
%   for the weight function w(x)=x^A*exp(-x)*log(1+1/x) on [0,Inf]
%   using a 2-component discretization procedure. The variables 
%   DIG, A (a nonnegative integer) must be entered when prompted. 
%   The two required coefficients files, ab_explog_fin and ab_explog_inf, 
%   for the value of a entered must be uploaded and declared global
%   in the calling script.

global mc mp iq idelta irout eps0 ab_explog_fin ab_explog_inf dig a
dig=input('enter the value of dig');
a=input('enter the value of a');
Mmax=65; 
mc=2; mp=0; iq=1; idelta=2; irout=1; 
%
% The choice eps0=1 below forces the discretization procedures to
% execute only one iteration, which is enough since this single
% iteration already yields exact answers.
%  
eps0=1;
[ab,Mcap,kount]=smcdis(dig,N,eps0,@squad_explog,Mmax); 
ab=vpa(ab,nofdig);
%[Mcap kount]
