% TABLE2_27 Computations for the first example in Example 2.61
%
fprintf('\n')
f1='%8.0f %20.12f %19.12f %8.0f\n';
f2='%8.0f %20.12f %19.12f\n';
N=40; eps0=.5*10^(-8);
x=2; y=1; z=complex(x,y);
a=7.5; nu0=N+1; numax=200;
ab0=r_hermite(numax);
[ab,nu]=chri5(N,ab0,z,eps0,nu0,numax,rho0,1);
ab=real(ab);
disp('       n      modified Hermite        Hermite          nu')
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)));
  if n==5
    fprintf(f1,n,real(int),real(int0),nu)
  else
    fprintf(f2,n,real(int),real(int0))
  end
end
