% TABLE2_3 Computation of Table 2.3.
%
f1='%5.0f     ---     %9.2e %9.2e %9.2e %9.2e\n';
fprintf('\n')
disp('  Table 2.3 Condition number and error magnification for') 
disp('  Hermite measure')
fprintf('\n')
disp('    n  (2.1.36)    (2.1.44)  (2.1.46)    maga      magb')
iopt=1;
for n=2:3:14
  ab0=zeros(2*n,2);
  mom=zeros(1,2*n);
  r=1:2:2*n-1; mom(r)=gamma(r/2);
  ab=r_hermite(n);
  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,rc,gc,maga,magb)
end
