%R_EXPLOG

function ab=r_explog(N) 
global mc mp iq idelta irout AB
mc=1; mp=0; iq=0; idelta=1; irout=1; AB=[0 Inf];
eps0=.5e-10; Mmax=500;
[ab,Mcap,kount]=mcdis(N,eps0,@quadgp,Mmax);
[Mcap kount]
