function [ab,Mcap,kount]=sr_hypexp(nofdig,N)
%SR_HYPEXP Recurrence coefficients for a hyperexponential weight
%function
%   [AB,MCAP,KOUNT]=SR_HYPEXP(NOFDIG,N) generates in DIG-digit 
%   arithmetic the first N recurrence coefficients for the weight 
%   function w(x)=exp(-exp(x)) on [0,Inf], using a 4-component
%   discretization procedure, and places them to NOFDIG digits into 
%   the Nx2 array AB. The variable DIG must be entered when prompted. 
%   The output variables MCAP, KOUNT are inherited from the routine 
%   smcdis.m.

global mc mp iq idelta irout uv AB dig eps0 
dig=input('enter the value of dig')
digits(dig);
eps0=.5*10^(-nofdig);
mc=4; mp=0; iq=0; idelta=2; irout=1; Mmax=2000;
AB=[[vpa(0) vpa(3/4)];[vpa(3/4) vpa(3/2)];[vpa(3/2) vpa(3)]; ...
   [vpa(3) Inf]];
[ab,Mcap,kount]=smcdis(dig,N,eps0,@squadgp,Mmax);
ab=vpa(ab,nofdig);
