% TABLE2_22  Computations for Table 2.22.
f1='                  n          beta_n          m=%2.0f\n';
f2='%19.0f %23.15e\n';
fprintf('\n')
disp('                   Induced Legendre polynomials') 
fprintf('\n')
N=20; M=11;
ab0=r_jacobi(N+M);
for m=[0 2 6 11]
  fprintf(f1,m)
  ab=indop(N,m,ab0);
  n=[1 2 7 13 20];
  fprintf(f2,[n'-1 ab(n,2)]') 
  fprintf('\n')
end
