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

digits(dig); mom=vpa(zeros(1,2*N));
s1=vpa(1); s2=vpa(2); s3=vpa(3); 
mom(1)=s2*zeta(s3);
for k=2:2*N
  skm1=vpa(k-1); f=prod(s1:skm1);
  skp1=vpa(k+1); skp2=vpa(k+2); 
  nu=vpa(2:k); s=sum(zeta(nu).*zeta(skp2-nu));
  mom(k)=f*(skp1*zeta(skp2)-s);
end
