% TABLE2_15 Computations for Table 2.15.
%
f1=' %4.0f %18.10e %18.10e %10.2e %9.2e\n';
f2=' %53.2e %9.2e\n';
fprintf('\n')
disp('                     Results for Example 2.41')
N=40; Mmax=400; eps0=1000*eps;
load -ascii abhrhermite;
abex=abhrhermite(1:N,:); abex(1,2)=sqrt(pi)/2;
ab=r_hrhermite(N,Mmax,eps0); 
fprintf('\n')
disp('    n       alpha_n             beta_n          erra      errb')
errabm=max(abs(ab-abex)./moll(abex));
n=[1 2 7 16 27 40];
errab=abs((ab(n,:)-abex(n,:))./moll(abex(n,:)));
abe=[(n-1)' abex(n,:) errab];
fprintf(f1,abe')
fprintf(f2,errabm')
