% TABLE2_23 Complex exponential integral
%
fprintf('\n')
f='%13.2f %9.4f %18.12f %17.12f %6.0f\n';
N=0; eps0=.5*10^(-12);
numax=620; 
ab=r_laguerre(numax);
disp('          r       phi          real(E1)          imag(E1)        nu')
for r=[.5 1 2 5]
  for phi=[0 pi/4 pi/2 3*pi/4]
    z=r*exp(complex(0,phi));
    nu0=10;
    [rho,rm1,nu]=cauchy(N,ab,-z,eps0,nu0,numax);
    E1=-exp(-z)*rm1;
    fprintf(f,r,phi,real(E1),imag(E1),nu)
  end
  fprintf('\n')
end
