%R_LSRBINET Recurrence coefficients for the lower sub-range
%Binet weight function
function ab=r_lsrbinet(N)
global c ablog abjac mc mp iq idelta irout
c=1/100; N=65;
ablog=loadvpa('coeff_lalg0log1.txt',9,32,200,2);
ablog=double(ablog);
abjac=r_jacobi01(200);
mc=3; mp=0; iq=1; idelta=2; irout=1; 
eps0=.5e-13; Mmax=200;
[ab,Mcap,kount]=mcdis(N,eps0,@quad_lsrbinet,Mmax);
[Mcap,kount]
for k=1:65
  if k==1
    fprintf('%17.13f %17.13f\n',ab(k,1),ab(k,2))
  else
    fprintf('%17.13f %18.14f\n',ab(k,1),ab(k,2))
  end
end
