function mom=smom_sqbinet(dig,N)
%SMOM_SQBINET Moments of the square Binet weight function.
%   MOM=SMOM_SQBINET(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 [-Inf,Inf]. 

digits(dig); mom=vpa(zeros(1,2*N));
s1=vpa(1); s2=vpa(2); s3=vpa(3); s4=vpa(4); 
mom(1)=s4*zeta(s3);
for k=3:2:2*N-1
  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)=vpa(2)*f*(skp1*zeta(skp2)-s);
end
