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

global mc mp iq idelta irout a ablag
a=input('enter the value of a');
ablag=r_laguerre(2000);
mc=1; mp=0; iq=0; idelta=2; irout=0; Mmax=2000; 
eps0=.5e-12;
[ab,Mcap,kount]=mcdis(N,eps0,@quad_algexp,Mmax); [Mcap kount]
