% EX1_IAL  Example 1 for Ialg and Ilog 
%
global abj abl
f0='%8.0f %18.14f %18.14f\n';
disp('       n        Ialg               Ilog')
alpha=-1/2;
abj=r_jacobi01(15,0,alpha); 
abl=r_jacobilog(15,0,alpha);
for n=1:15
  [Ialg,Ilog]=Ial(n);
  fprintf(f0,n,Ialg,Ilog)
end
