function ab=sr_explogalgm1half01(nofdig,N)
%SR_EXPLOGALGM1HALF01 Recurrence coefficients for a weight 
%function on [0,1] with an algebraic/extended-logarithmic 
%singularity with exponent -1/2.
%   AB=SR_EXPLOGALGM1HALF01(NOFDIG,N) generates in DIG-digit
%   arithmetic the first N recurrence coefficients for the 
%   weight function w(x)=x^(-1/2)*exp(-x)*log(1+1/x) on [0,1],  
%   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 abjac ablog
dig=input('enter the value of dig');
abjac=sr_jacobi01(32,200,0,-1/2);
ablog=loadvpa('coeff_lalgm1halflog1.txt',9,32,200,2);
Mmax=200; 
mc=2; mp=0; iq=1; idelta=2; irout=1; eps0=.5*10^(-nofdig); 
[ab,Mcap,kount]=smcdis(dig,N,eps0,@squad_explogalg01,Mmax); 
digits(dig); ab=vpa(ab,nofdig);
[Mcap kount]
