function ab=sr_exp_abc(dig,nofdig,N,a,b,c)
%SR_EXP_ABC Recurrence coefficients for orthogonal polynomials
%relative to a three-parameter exponential weight function.
%   AB=SR_EXP_ABC(DIG,N,A,B,C) 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)=exp(-B/x^A-C*x^A on [0,Inf],
%   A>0, B>0, C>0.

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