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