function [ab,Ncap,kount]=r_symm_subrange_jacobi(N,eps0)
%R_SYMM_SUBRANGE_JACOBI Recurrence coefficients for symmetric
%subrange Jacobi polynomials.
%   [AB,NCAP,KOUNT]=R_SYMM_SUBRANGE_JACOBI(N,EPS0) generates to
%   to an accuracy EPS0 the Nx2 array AB of recurrence coefficients
%   for the weight function
%               w(x(t;c)),  -1<=t<=1, 0<c<=1,
%   where w is the Jacobi weight function with parameters alpha,
%   beta, and 
%                       x(t;c)=c*t.
%   NCAP and KOUNT are output variables inherited from the routine
%   mcdis. The variables c, alpha,beta must be assigned in the 
%   script calling this routine and declared global. If an error
%   message appears, try to increase Nmax.

global mc mp iq idelta irout AB c alpha beta
mc=1; mp=0; iq=0; idelta=1; irout=1; AB=[-1 1];
Nmax=3000;
[ab,Ncap,kount]=mcdis(N,eps0,@quadgp,Nmax);
