% MOMHALFRANGEHERMITE Generates the first 2N moments 
%    to dig decimal places of the half-range Hermite  
%    weight function w(x)=exp(-x^2) on R_{+}.
%
function mf=halfrangehermite(dig,N)
digits(dig); 
for k=1:2*N
  mf(k)=vpa(gamma(vpa(k/2))/2);
end