%R_EXPLOG01

function ab=r_explog01(N) 
global mc mp iq idelta irout AB
mc=3; mp=0; iq=0; idelta=1; irout=1;
AB=[[0 .0001];[.0001 .01];[.01 1]];
eps0=.5e-14; Mmax=500;
[ab,Mcap,kount]=mcdis(N,eps0,@quadgp,Mmax);
[Mcap kount]
