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

digits(dig); 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,s1);
end
