%R_ALGLOG2
function ab=r_alglog2(n)
global abjac ablag1 ablag2 ableg mc mp iq idelta irout a xi
mc=4; mp=0; iq=1; idelta=2; irout = 1; Mmax=500; eps0=.5e-12;
[ab,Mcap,err,kount]=mcdis(n,eps0,@quad_alglog2,Mmax);
[Mcap err]
