function ab=sr_laglog(dig,nofdig,N,a)
%SR_LAGLOG  Recurrence coefficients for a generalized
%Laguerre weight function multiplied by a logarithmically
%singular function.
%   AB=SR_LAGLOG(DIG,NOFDIG,N,A) generates in DIG-digit 
%   arithmetic 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=momlaglog(dig,N,a);
ab=schebyshev(dig,N,mom);
ab=vpa(ab,nofdig);
