function ab=r_algexpabslog(N,a)
%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 eps0
a=input('enter the value of a')
mc=2; mp=0; iq=1; idelta=2; irout=1; Mmax=500; 
eps0=.5e-8;
[ab,Mcap,koUnt]=mcdis(N,eps0,@quad_algexpabslog,Mmax);
