% TABLE2_28 Computations for the second example in Example 2.61
%
fprintf('\n')
f='%8.0f %20.12f %19.12f\n';
N=40; x=0; y=.1; a=7.5;
ab0=r_hermite(N);
rho0=complex(0,-pi*erfcx(y));
[ab,nu]=chri6(N,ab0,y,eps0,nu0,numax,rho0,2);
disp('       n      modified Hermite        Hermite')
for n=5:5:40
  xw0=gauss(n,ab0);
  xw=gauss(n,ab);
  int0=sum(xw0(:,2).*bessel(0,a*xw0(:,1))./((xw0(:,1)-x).^2+y^2));
  int=sum(xw(:,2).*bessel(0,a*xw(:,1)));
  fprintf(f,n,real(int),real(int0))
end
