%R_SUBJACOBI0
%  This routine generates to a relative 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 and
%            x(t;c)=.5(1+c)t-.5(1-c).
%  If an error message appears, try to increase Nmax.
% 
function [ab,Ncap]=r_subjacobi0(N,eps0,c,alpha,beta)
global mc mp iq idelta irout AB
global a b cg
a=alpha; b=beta; cg=c;
mc=2; mp=0; iq=1; idelta=1; irout=1; AB=[[-1 0];[0 1]];
Nmax=3000;
[ab,Ncap,kount]=mcdis(N,eps0,@quadsubj,Nmax);
