function showquads ( )

%% showquads() uses quad to estimate the integral of the
% humps function from 0 to 1.
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'showquads()\n' );

    clc
	home
    disp('Tolerance       quad       N-quad')
    disp(' ')
    for tol = logspace(-2,-6,5)
      [Q1,count1] = quad('humps',0,1,tol);
      s1 = sprintf('  %12.7f  %5.0f',Q1,count1);
      disp([sprintf('%8.6f ',tol) s1 ])  
   end

  return
end

