function ab=sr_usrbinet_old(nofdig,N,Ngauss)
%SR_USRBINET_OLD Recurrence coefficients for the upper subrange
%Binet weight function.
%   AB=SR_USRBINET_OLD(NOFDIG,N,Ngauss) generates in DIG-digit 
%   arithmetic the Nx2 array AB of the first N recurrence coefficients for
%   for the weight function w(x)=-log(1-exp(-x)) on [1,Inf] using
%   a 1-component discretization procedure. The variable DIG must
%   be entered when prompted. 

global dig eps0 ablag mc mp iq idelta irout
tic
dig=input('enter the value of dig');
Mmax=2000; ablag=sr_laguerre(dig,Mmax);
mc=1; mp=0; iq=1; idelta=2; irout=1; 
eps0=.5*10^(-nofdig); 
[ab,Mcap,kount]=smcdis(dig,N,eps0,@squad_usrbinet,Mmax); [Mcap kount]
toc
