function ab=sr_radtrans_dis(nofdig,N)
%SR_RADTRANS_DIS Recurrence coefficients for a radiative transfer
%weight function.
%   AB=SR_RADTRANS_DIS(NOFDIG,N) uses a 1-component discretization
%   procedure to compute in dig-digit arithmetic the first N
%   recurrence coefficients for the radiative transfer weight
%   function on [0,c] and places them to NOFDIG digits into the
%   Nx2 array AB. The values of c and dig have to be entered
%   when prompted. 

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