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

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