function ab=sr_gegexp(nofdig,N)
%SR_GEGEXP Recurrence coefficients for a Gegenbauer weight
%function multiplied by an exponential function.
%   AB=SR_SEGEXP(DIG,NOFDIG,N,LAM,A) uses a 1-component
%   discretization procedure to compute in DIG-digit arithmetic 
%   the first N recurrence coefficients for the weight function 
%   w(x)=(1-x^2)^(LAM-1/2)*exp(-A^x^2) on [-1,1] and places them 
%   to NOFDIG digits into the Nx2 array AB. The values of the 
%   variables DIG, LAM, A have to be entered when prompted. 

global mc mp iq idelta irout dig lam a
dig=input('enter the value of dig')
lam=input('enter the value of lam')
a=input('enter the value of a')
mc=1; mp=0; iq=1; idelta=2; irout=1; Mmax=500;
eps0=.5*10^(-nofdig);
[ab,Mcap,kount]=smcdis(dig,N,eps0,@quad_sgegexp,Mmax); Mcap
ab=vpa(ab,nofdig);
