% EXAMPLE3_59 Computation for Example 3.59.
%
f1='%12.0f %4.0f %12.4e\n';
fprintf('\n')
disp('           m    n      err')
ta=zeros(80,2);
for m=1:3
  ab=r_laguerre(80,m+1); ab(1,2)=ab(1,2)/prod(1:m);
  for n=[5 10 20 40 80]
    xw=gauss(n,ab);
    ta(1:n,1)=xw(:,1); ta(1:n,2)=xw(:,2)./(xw(:,1).^(m+1));
    x=linspace(0,ta(n,1));
    errmax=max(abs(splinefunc(n,m,x,ta)-exp(-x)));
    fprintf(f1,m,n,errmax)
  end 
  fprintf('\n')
end
