% EXAMPLE2_64 Routine used in Example 2.64.
%
N=12; s=1; same=0; eps0=1000*eps; numax=1000; nd=[N N];
alpha=.5;
for xi=[.0914 .0913 .09 .07 .03]
  xi
  ab0=r_jacobi(numax,alpha);
  z=complex(0,xi);
  nu0=nu0jac(0,z,eps0); iopt=1; rho0=0;
  ab1=chri6(N,ab0,xi,eps0,nu0,numax,rho0,iopt); 
  zw1=gauss(N,ab1);
  ab=r_jacobi(N,alpha-1); zw=gauss(N,ab);
  xw=[zw(:,1) zw1(:,1) zw(:,2) zw1(:,2)];
  a0=ab(1,1);
  B=stieltjes_sob(N,s,nd,xw,a0,same);
  z=sobzeros(N,N,B)
end
