% MOMFREUD Generates the first 2N moments to dig decimal places
%    of the Freud weight function 
%    w(x)=|x|^alpha*exp(-|x|^beta) on R
%
function mf=momfreud(dig,N,alpha,beta)
digits(dig); 
for k=1:2*N
  if rem(k,2)==0
    mf(k)=0;
  else
    mf(k)=vpa(2*gamma(vpa((k+alpha)/beta))/beta);
  end
end