% TEST_IAL  Test of Ial
%
global abj abl
f0='%4.0f %15.12f %15.12f %10.2e %9.2e\n';
disp('   n      Ialg             Ilog           erra      errl')
for alpha=[-3/4 -1/2 0 1/2 3/4]
  abj=r_jacobi01(3,0,alpha); abl=r_jaclog(3,alpha);
  Ialgex=2/((alpha+1)*(alpha+2));
  Ilogex=-2*(2*alpha+3)/((alpha+1)^2*(alpha+2)^2);
  for n=1:3
    [Ialg,Ilog]=Ial(n);
    errIalg=abs((Ialg-Ialgex)/Ialgex);
    errIlog=abs((Ilog-Ilogex)/Ilogex);
    fprintf(f0,n,Ialg,Ilog,errIalg,errIlog)
  end
  fprintf('\n')
end
