function mom=smom_hrgbinet(dig,N,a)
%SMOM_HRGBINET Moments of the half-range generalized Binet 
%weight function. 
%   MOM=SMOM_HRGBINET(DIG,N,A) generates in DIG-digit arithmetic 
%   the 1x2N array MOM of the first 2N moments of the weight
%   function w(x)=-log(1-A*exp(-x)) on [0,Inf], 0 < A <= 1. 

digits(dig); mom=vpa(zeros(2*N,1)); a=vpa(a); s1=vpa(1); 
for k=1:2*N
  skm1=vpa(k-1); skp1=vpa(k+1); f=prod(s1:skm1); 
  mom(k)=f*polylog(skp1,a);
end
