% TABLE2_29 Computations for Table 2.26.
%
f1='%6.1f %15.10f %15.2e %24.10f\n';
fprintf('\n')
disp('           The behavior of beta_j^k for large gamma')
fprintf('\n')
n=21; g=100;
mom=zeros(2,2*n);
mom(1,1)=2; mom(2,1)=2*g;
abm=r_jacobi(2*n-1);
abjac=r_jacobi(n-2,1,1);
[B,normsq]=chebyshev_sob(n,mom,abm);
disp('    k       beta_1^k    max_{j>=2}|beta_j^k|    beta_{k-2}^{(1,1)}')
for k=5:5:20
  j=2:k;
  m=max(abs(B(j+1,k+1)));
  fprintf(f1,k,B(2,k+1),m,abjac(k-1,2))
end
