function ab=sr_explogalgm1half0inf(nofdig)
%SR_EXPLOGALGM1HALF0INF Recurrence coefficients for a weight 
%function on [0,Inf] with an algebraic/extended-logarithmic 
%singularity with exponent -1/2.
%   AB=SR_EXPLOGALGM1HALF0INF(NOFDIG) generates in DIG-digit
%   arithmetic the first N=65 recurrence coefficients for the 
%   weight function w(x)=x^(-1/2)*exp(-x)*log(1+1/x) on [0,Inf],  
%   using a 2-component discretization procedure, and places 
%   them to NOFDIG digits into the Nx2 array AB. The variable 
%   DIG must be entered when prompted by the routine. 

global mc mp iq idelta irout dig ab_algexplog_fin ab_algexplog_inf
dig=input('enter the value of dig');
N=65; Mmax=65; 
ab_algexplog_fin=loadvpa('coeff_explogalgm1half01.txt',9,32,65,2);
ab_algexplog_inf=loadvpa('coeff_explogalgm1half1inf.txt',9,32,65,2);
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_algexplog,Mmax); 
ab=vpa(ab,nofdig);
[Mcap kount]
