function ab=sr_exp_abcd(dig,nofdig,N,a,b,c,d)
%SR_EXP_ABCD Recurrence coefficients for orthogonal polynomials
%relative to a four-parameter exponential weight function.
%   AB=SR_EXP_ABDC(DIG,N,A,B,C,D) generates in DIG-digit
%   arithmetic the Nx2 array AB of the first N recurrence
%   coefficients for polynomials orthogonal with respect to
%   the weight function w(x)=x^D*exp(-B/x^A-C*x^A) on [0,Inf],
%   A>0, B>0, C>0, D real.

digits(dig); a=vpa(a); b=vpa(b); c=vpa(c); d=vpa(d);
mom=smom_exp_abcd(dig,N,a,b,c,d);
ab=schebyshev(dig,N,mom);
ab=vpa(ab,nofdig);
