function ab=sr_recgamma(dig,N)
%SR_RECGAMMA Recurrence coefficients for the reciprocal gamma function.
%   AB=SR_RECGAMMA(DIG,N) generates in DIG-digit arithmetic 
%   the Nx2 array AB of the first N recurrence coefficients 
%   for the orthogonal polynomials relative to the weight function 
%   W(X)=1/GAMMA(X) on [0,INF] using a multicomponent discretization
%   procedure.

global mc mp iq idelta irout DM uv AB dig
mc=10; mp=0; iq=0; idelta=1; irout=1; Mmax=1000;
eps0=.5*10^(-dig);
AB=[[vpa(0) vpa(5)];[vpa(5) vpa(10)];[vpa(10) vpa(15)];[vpa(15) vpa(20)]; ...
  [vpa(20) vpa(25)];[vpa(25) vpa(30)];[vpa(30) vpa(35)];[vpa(35) vpa(40)]; ...
  [vpa(40) vpa(45)];[vpa(45) Inf]];
ab=smcdis(dig,N,eps0,@squadgp,Mmax);
