% EXAMPLE3_19 Routine used in Example 3.19.
%
f1='%9.0f %10.2e %10.1e %16.8e %9.2e\n';
fprintf('\n')
a=3.1;
fprintf('                           alpha= %3.3f\n',a)
disp('        n   min tauK  interlacing    max tauK     min lam')
for N=1:125
  ab=r_jaclog(ceil(3*N/2+1),a);
  xwg=gauss(N,ab);
  xwk=kronrod(N,ab);
  fprintf(f1,N,xwk(1,1),max(abs(xwg(:,1)-xwk(2:2:2*N,1))),...
    xwk(2*N+1,1),min(xwk(:,2)))
end
