% TABLE2_2 Computation of Table 2.2.
%
f1='%5.0f %9.2e %9.2e %9.2e %9.2e %9.2e\n';
fprintf('\n')
disp('  Table 2.2 Condition number and error magnification for') 
disp('  Laguerre measure with parameter alpha=.5')
fprintf('\n')
disp('    n  (2.1.36)  (2.1.44)  (2.1.46)    maga      magb')
a=-.5; iopt=1;
for n=2:3:14
  ab0=zeros(2*n,2);
  mom=zeros(1,2*n);
  r=1:2*n; mom(r)=gamma(r+a);
  ab=r_laguerre(n,a);
  acl=acondGlow(n,ab);
  rc=rcondG(n,iopt,ab,ab0,mom);
  gc=gcondG(n,ab,mom);
  [abcheb,normsq]=chebyshev(n,mom);
  errab=abs(abcheb-ab)./moll(ab);
  maga=max(errab(:,1)/eps);
  magb=max(errab(:,2)/eps);
  fprintf(f1,n,acl,rc,gc,maga,magb)
end
