%SMOMENUFIN Symbolic counterpart of momEnufin.m
%
function mom=smomEnufin(dig,n,nu,c)
digits(dig);
c=vpa(c,dig);
mom=vpa(zeros(2*n,1));
y=feval(symengine,'igamma',1-nu,c);
y=vpa(y,dig);
for k=1:2*n
  y1=gamma(vpa(k,dig))*c^k*sgammastar(dig,k,c);
  mom(k)=(y1+c^(nu+k-1)*y)/vpa(nu+k-1,dig);
end
