function ab=sr_laglogmm(dig,nofdig,N,a)
%SR_LAGLOGMM  Recurrence coefficients for a generalized
%Laguerre weight function multiplied by a logarithmically
%singular weight function.
%   AB=SR_LAGLOGMM(DIG,NOFDIG,N,A) generates in DIG-digit 
%   arithmetic, using modified moments, the first N recurrence 
%   coefficients for polynomials orthogonal with respect to the 
%   weight function w(x)=x^A*exp(-x)*(x-1-log(x)) on [0,Inf], 
%   A>-1, and places them to NOFDIG digits into the Nx2 array AB. 
  
digits(dig);
mom=mmomlaglog(dig,N,a);
abm=sr_laguerre(dig,2*N-1,a);
ab=schebyshev(dig,N,mom,abm);
ab=vpa(ab,nofdig);
