function mom=smom_binet(dig,N)
%SMOM_BINET Moments of the Binet weight function.
%   MOM=SMOM_BINET(DIG,N) generates in DIG-digit arithmetic
%   the (2N)x1 array MOM of the first 2N moments of the weight
%   function w(x)=-log(1-exp(-|x|) on [-Inf,Inf].

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