function ab=sr_algexp(nofdig,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 dig ablag
tic
dig=input('enter the value of dig');
a=input('enter the value of a');
Mmax=2000; ablag=sr_laguerre(dig,Mmax);
mc=1; mp=0; iq=0; idelta=2; irout=0; 
eps0=.5*10^(-nofdig);
[ab,Mcap,kount]=smcdis(dig,N,eps0,@squad_algexp,Mmax); [Mcap kount]
digits(dig); ab=vpa(ab,nofdig);
toc
