%SR_EXPLOG

function ab=sr_explog(dig,N) 
global mc mp iq idelta irout AB 
digits(dig); 
mc=5; mp=0; iq=0; idelta=1; irout=1;
AB=[[vpa(0) vpa(1/10000)];[vpa(1/10000) vpa(1/1000)]; ...
[vpa(1/1000) vpa(1/100)];[vpa(1/100) vpa(1/10)];[vpa(1/10) vpa(1)]];
eps0=.5e-32; Mmax=5000;
[ab,Mcap,kount]=smcdis(dig,N,eps0,@squadgp,Mmax);
[Mcap kount]
