function ab=sr_schroedinger(dig,N,nofdig,mu)
%SR_SCHROEDINGER Schroedinger weight function.
%   AB=SR_SCHROEDINGER(DIG,N,NOFDIG,MU) generates in DIG-digit   
%   arithmetic the first N recurrence coefficients for the
%   weight function w(x)=x^mu*exp(-x^4/16) on [0,Inf], and
%   places them to NOFDIG digits into the Nx2 array AB.

digits(dig); 
ab=vpa(zeros(N,2)); s2=vpa(2); s4=vpa(4);
ab0=sr_hrfreud(156,nofdig,100,mu,4);
ab(:,1)=s2*ab0(:,1);
ab(1,2)=vpa(2^(mu+1))*ab0(1,2); ab(2:N,2)=s4*ab0(2:N,2);
ab=vpa(ab,nofdig);
