%R_ALGLOG1
function ab=r_alglog1(n)
global abjac ablag1 ableg mc mp iq idelta irout a xi
mc=3; mp=0; iq=1; idelta=2; irout=1; Mmax=400; eps0=.5e-12;
[ab,Mcap,err,kount]=mcdis(n,eps0,@quad_alglog1,Mmax);
[Mcap err]
