% SMOMVAB  The first 2N moments to dig decimal places of
%    the weight function 
%    v(x;a,b)=x^{a-1}[ln(1/x)]^{a-b+1},  0 < x \leq 1/e
%
function smom=smomvab(dig,N,a,b)
digits(dig);
if subs(a)<=0
  error('a not in range')
  return
end
for k=1:2*N
  g=feval(symengine,'igamma',a-b+2,k+a-1);
  igam=vpa(g);
  smom(k)=igam/(vpa(k+a-1))^vpa(a-b+2);
end