function mom=smom_hrh(dig,N)
%SMOM_HRH Moments of the half-range Hermite weight function.
%   MOM=SMOM_HRH(DIG,N) generates the first 2N moments in DIG-digit
%   arithmetic of the weight function W(X)=EXP(-X^2) on [0,INF].
   
digits(dig); mom=vpa(zeros(2*N,1)); 
for k=1:2*N
  mom(k)=gamma(vpa(k/2))/2;
end
