function mom=smom_gprudnikov(dig,N,nu,alpha)
%SMOM_GPRUDNIKOV Moments of the generalized Prudnikov weight function.
%   MOM=SMOM_GPRUDNIKOV(DIG,N,NU,ALPHA) generates in DIG-digit 
%   arithmetic the first 2N moments (the first moment if N=0) of the
%   generalized Prudnikov weight function W(X), where 
%                W(X)=X^ALPHA*RHO_NU(X)
%   and
%
%          RHO_NU(X)=2*X^(NU/2)*K_NU(2*sqrt(X))  on [0,INF],
%
%   where K_NU is the second-kind modified Bessel function of real
%   order NU.

digits(dig);
if N>0
  mom=vpa(zeros(2*N,1));
end
vnu=vpa(nu); valpha=vpa(alpha); v1=vpa(1);
if N==0
  mom=gamma(valpha+vnu+v1)*gamma(alpha+v1);
else
  for k=1:2*N
    km1=k-1; vkm1=vpa(km1);
    mom(k)=gamma(vkm1+valpha+vnu+v1)*gamma(vkm1+valpha+v1);
  end
end
