function mom=smom_sqgbinet(dig,N,a)
%SMOM_SQGBINET Moments of the squared generalized Binet weight 
%function.
%   MOM=SMOM_SQGBINET(DIG,N,A,NOFDIG) 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 [-Inf,Inf], 
%   0 < A <= 1.

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