function [ab,Mcap,kount]=sr_symm_subrange_jacobi(nofdig,N)
%SR_SYMM_SUBRANGE_JACOBI Recurrence coefficients for a 
%symmetgric subrange Jacobi weight function.
%   [AB,MCAP,KOUNT]=SR_SYMM_SUBRANGE_JACOBI(NOFDIG,N) uses a 
%   1-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. The output
%   variables MCAP, KOUNT are inherited from the routine smcdis.m.

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