function ab=sr_lsrjacobi(nofdig,N)
%SR_LSRJACOBI Recurrence coefficients for a lower 
%subrange Jacobi weight function.
%   AB=SR_LSRJACOBI(NOFDIG,N) uses a 2-component
%   discretization procedure to compute in DIG-digit arithmetic
%   the first N recurrence coefficients for the weight function 
%   w(x)=(1-x)^A*(1+x)^B on [-1,C], A>-1, B>-1, -1<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 
eps0=.5*10^(-nofdig);
%Mcap, kount
mc=2; mp=0; iq=1; idelta=2; irout=1;
digits(dig); eps0=.5*10^(-nofdig); Mmax=3000; 
ab0=smcdis(dig,N,eps0,@squad_lsrjacobi,Mmax);
ab0=vpa(ab0,nofdig); sc=vpa(c,nofdig); 
s1=vpa(1,nofdig); s2=vpa(2,nofdig); s1h=vpa(1/2,nofdig);
ab(:,1)=((s1+sc)*ab0(:,1)-(s1-sc))/s2;
ab(1,2)=s1h*(s1+sc)*ab0(1,2);
k=(2:N)'; ab(k,2)=(s1h*(s1+sc))^2*ab0(k,2);
