% EXAMPLE3_44 Computations for Example 3.44.
%
disp('                              Example3.44')
fprintf('\n')
f1='%4.0f %18.14f %11.3e %18.14f %11.3e\n';
Nmax=15; ab=r_jacobi(Nmax+1);
x=1/3; exact=.97621001086991718;
rho0=log((1+x)/(1-x));
disp('   n       (3.1.127)         err          (3.1.132)         err')
for N=2:Nmax
  y1=cauchyPVI(N,x,'fex44','ddfex44',1,ab,rho0);
  y0=cauchyPVI(N,x,'fex44','ddfex44',0,ab,rho0);
  err1=abs(y1-exact); err0=abs(y0-exact);
  fprintf(f1,N,y1,err1,y0,err0)
end
