function ab=r_dblexp(N,al,be,ga,a,b)
%    ab=R_DBLEXP(N) generates the first N recurrence coefficients
%    for the monic polynomials orthogonal on [0,Inf] relative to
%    the weight function a(t)=t^ga*exp(-a*t^(-al)-b* ^be), ga>0,
%    a>0, b>0, al>0, be>0.
%
global mc mp iq idelta irout eps0 Mmax Mcap
f0='%5.1f %5.1f %8.4f %5.0f\n';
mc=2; mp=0; iq=1; idelta=2; irout=0;
[ab,Mcap,kount]=mcdis(N,eps0,@quaddblexp,Mmax);
%t0=(a*al/(b*be))^(1/(al+be));
%fprintf(f0,al,be,t0,Mcap)
%if Mcap>Mmax
%  fprintf(f0,al,be,t0,Mcap)
%  error('quit')
%end
