function quadtx_test ( )

%*****************************************************************************80
%
%% quadtx_test tests quadtx.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    08 January 2020
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'quadtx_test\n' );
  fprintf ( 1, '  quadtx estimates integral (a<=x<=b) f(x) dx.\n' );

  f = @(x)1.0 ./ sqrt ( 1.0 + x.^4 );
  a = 0.0;
  b = 1.0;

  value = quadtx ( f, a, b );

  fprintf ( 1, '\n' );
  fprintf ( 1, '  Estimate integral of sin(x) between %g and %g\n', a, b );
  fprintf ( 1, '  Integral estimate = %g\n', value );

  return
end
