% EXAMPLE3_70  Computations for Example 3.70.
%
f1='%10.0f %8.3f %5.0f %20.14f\n';
f2='%10.0f %8.3f %5.0f %17.11f\n';
f3='%10.0f %8.3f %5.0f %12.6f\n';
fprintf('\n')
x=[.8 .9 .95 .99 .999 1.0];
eps0=10^2*eps; nu0=5; numax=100;
disp('         p      x      nu         R_p(x)')
for k=1:6
  z=x(k)^(-2);
  for p=2:3
    c=2^p*factorial(p-1);
    if p==2
      load -ascii absqm1log1;
      ab=absqm1log1;
    else
      load -ascii absqm1log2;
      ab=absqm1log2;
    end
    if k==5, eps0=10^5*eps; end
    if k==6, eps0=10^(10)*eps; end
    [rho,r,nu]=cauchy(0,ab,z,eps0,nu0,numax);
    rho(1)=rho(1)/(c*x(k));
    if k<5
      fprintf(f1,p,x(k),nu,rho(1))
    elseif k==5
      fprintf(f2,p,x(k),nu,rho(1))
    else
      fprintf(f3,p,x(k),nu,rho(1))
    end
  end
end
