function ab=sr_hrsqgbinet(dig,nofdig,N,a)
%SR_HRSQGBINET  Recurrence coefficients for the half-range squared 
%generalized Binet weight function.
%   AB=SR_HRSQGBINET(DIG,NOFDIG,N,A) generates in DIG-digit arithmetic
%   the first N recurrence coefficients for orthogonal polynomials   
%   relative to the weight function w(x)=[log(1-A*exp(-x))]^2 on
%   [0,Inf], 0 < A < 1, using a moment-based method. The recurrence
%   coefficients to NOFDIG digits are placed into the Nx2 array AB.

mom=smom_hrsqgbinet(dig,N,a);
ab=schebyshev(dig,N,mom');
ab=vpa(ab,nofdig);
