%EXAMPLE5
%
f0='%10.0f %6.0f %21.15f\n';
f1='%16.0f %21.15f\n';
global r imod dig
dig=28; digits(dig);
N=12; 
ab0=sr_laguerre(dig,N);
disp('         r      n            I') 
for r=2:4
  z=vpa(2*r*pi*i);
%  z=vpa(-1);
  imod=1;
  ab=mod_d(N,r,z,ab0,@rho0_mu);
  if imag(subs(z))~=0
    imod=2;
    ab=mod_d(N,r,z,ab,@rho0_mu);
  end
  ab=real(subs(ab));
  for n=2:2:N
    xw=gauss(n,ab); x=xw(:,1); w=xw(:,2);
    I=sum((1/r)*w.*exp(-x/r).*f5(x,r));
    fprintf(f0,r,n,I)
  end
  fprintf('\n')
end
