function ab=sr_log1pe(nofdig,N)
%SR_LOG1PE Recurrence coefficients for a Binet-like weight
%function.
%   AB=SR_LOG1PE(NOFDIG,N) generates in DIG-digit arithmetic
%   the Nx2 array AB of the first N recurrence coefficients for
%   the weight function w(x)=log(1+exp(-abs(x))) on [-Inf,Inf]
%   using a 2-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');
mc=2; mp=0; iq=1; idelta=2; irout=1; Mmax=2000;
eps0=.5*10^(-nofdig); ablag=sr_laguerre(dig,Mmax);
ab=smcdis(dig,N,eps0,@squad_log1pe,Mmax);
ab=vpa(ab,nofdig);
toc
