%SR_EXLOG

function ab=sr_exlog(dig,nofdig,N,a)
global mc mp iq idelta irout dig a
dig=input('enter the value of dig');
a=input('enter the value of a');
abalglog=loadvpa('coeff_lalgm1halflog1.txt',9,32,100,2);
abj=sr_jacobi01(dig,Mmax,0,a);
ablag=sr_laguerre(dig,Mmax);
mc=3; mp=0; iq=1; idelta=2; irout=1; Mmax=200;
eps0=,5*10^(-nofdig);
[ab,Mcap,kount]=smcdis(dig,N,eps0,@quad_exlog,Mmax);
