function ab=sr_symm_subrange_jacobi(nofdig,N)
%SR_SYMM_SUBRANGE_JACOBI Recurrence coefficients for a 
%symmetgric subrange Jacobi weight function.
%   AB=SR_SYMM_SUBRANGE_JACOBI(NOFDIG,N) uses a 
%   4-component discretization procedure to compute in DIG-digit 
%   arithmetic he first N recurrence coefficients for the weight 
%   function w(x)=(1-x)^A*(1+x)^B on [-C,C], A>-1, B>-1, 0<C<1, 
%   and places them to NOFDIG digits into the Nx2 array AB. The 
%   variables DIG, A, B, C are assumed to have been assigned,
%   and declared global, in the calling program. 

global mc mp iq idelta irout dig a b c AB
eps0=.5*10^(-nofdig);
mc=4; mp=0; iq=0; idelta=1; irout=1; AB=[[-c -c/2];[-c/2 0];[0 c/2];[c/2 c]];
digits(dig); eps0=.5*10^(-nofdig); Mmax=3000;
ab=smcdis(dig,N,eps0,@squadgp_symm_subrange_jacobi,Mmax);
ab=vpa(ab,nofdig);
