%% QUAD_SCRIPT is a script to call QUAD_FUN.
%
%  Licensing:
%
%    This code is distributed under the GNU LGPL license.
%
%  Modified:
%
%    26 March 2010
%
%  Author:
%
%    John Burkardt
%
  n = 10000;  a = 0.0; b = 1.0;

  fprintf ( 1, '\n' );
  fprintf ( 1, 'QUAD_SCRIPT\n' );
  fprintf ( 1, '  Estimate integral using %d points\n', n );

  value = quad_fun( n, a, b );
  
  fprintf ( 1, '\n' );
  fprintf ( 1, '  Integral estimate is %f\n', value );
  fprintf ( 1, '  Exact integral    is %f\n', pi );
  fprintf ( 1, '  Error is             %e\n', abs ( value - pi ) );

