function ab=sr_radtrans_dis(dig,nofdig,N,c)
%SR_RADTRANS_DIS Recurrence coefficients for a radiative transfer
%weight function.
%   AB=SR_RADTRANS_DIS(DIG,NOFDIG,N,C) 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. In the script calling this routine, the variable
%   DIG must be assigned and declared global.

global mc mp iq idelta irout sxw_radtrans
ab=loadvpa('coeffradtrans.txt',9,32,100,2);
sxw_radtrans=sgauss(dig,N,ab);
mc=1; mp=0; iq=1; idelta=2; irout=1; Mmax=1000; eps0=1;
[ab,Mcap,kount]=smcdis(dig,N,eps0,@quad_sradtrans,Mmax);
