function ab=r_algexpabslog(N)
%R_ALGEXPABSLOG Recurrence coefficients for a weight function
%with an algebraic/exponential/absolutelogarithmic singularity.
%   AB=R_ALGEXPABSLOG(N,A) uses a 2-component discretization
%   procedure to compute the Nx2 array AB of the first N
%   recurrence coefficients for the weight function 
%   w(x)=x^A*exp(-x)*|log(x)| on [0,Inf]. The value of A has to 
%   be entered when prompted.

global mc mp iq idelta irout a abjl ablag
a=input('enter the value of a');
abjl=double(loadvpa('coeff_lalgm1halflog1.txt',9,32,100,2));
ablag=r_laguerre(100);
mc=2; mp=0; iq=1; idelta=2; irout=1; Mmax=100; 
eps0=.5e-8;
[ab,Mcap,kount]=mcdis(N,eps0,@quad_algexpabslog,Mmax); [Mcap kount]
