function ab=sr_hrfreud(dig,nofdig,N,mu,nu)
%SR_HRFREUD Recurrence coefficients for the half-range Freud 
%weight function.
%   AB=SR_HRFREUD(DIG,NOFDIG,N,MU,NU) generates in DIG-digit 
%   arithmetic the first N recurrence coefficients for the 
%   weight function w(x)=x^MU*exp(-x^NU) on [0,Inf], MU>-1,
%   NU>0, and places them to NOFDIG digits into the Nx2 
%   array AB.
 
digits(dig);
mom=smom_hrfreud(dig,N,mu,nu);
ab=schebyshev(dig,N,mom');
ab=vpa(ab,nofdig);
