% FIGURE2_3 Computations to generate Figure 2.3.
%
N=20;
numax=200; eps0=1e-5; 
ab=r_jacobi(numax);
x=0:.1:10;
y=.1:.1:10;
[X,Y]=meshgrid(x,y);
for ix=1:100
  for jy=1:101
    xx=X(ix,jy); yy=Y(ix,jy);
    z=complex(xx,yy);
    nu0=nu0jac(N,z,eps0);
    ker=cauchy(N,ab,z,eps0,nu0,numax);
    Z(ix,jy)=abs(ker(11));
  end
end
v=[1e-4 1e-6 1e-8 1e-10 1e-12 1e-14 1e-16];
contour(x,y,Z,v)
axis([0 10 0 10])
text(7.3,5.5,'10^{-16}','FontSize',12)
text(4.7,4,'10^{-14}','FontSize',12)
text(3.1,2.65,'10^{-12}','FontSize',12)
text(2.25,1.65,'10^{-10}','FontSize',12)
text(1.65,.9,'10^{-8}','FontSize',12)