#*****************************************************************************80
#
## adaptint_test tests adaptint.
#
#  Licensing:
#
#    Copyright 2016 James P. Howard, II
#
#    The computer code and data files on this web page are distributed under
#    https://opensource.org/licenses/BSD-2-Clause, the BSD-2-Clause license.
#
#  Modified:
#
#    07 March 2020
#
#  Author:
#
#    John Burkardt
#
cat ( date ( ), "\n" )
cat ( "\n" )
cat ( "adaptint_test\n" )
cat ( " ", version$version.string, "\n" )
cat ( "  adaptint estimates an integral using an adaptive midpoint method.\n" )

source ( "/home/burkardt/public_html/r_src/adaptint/adaptint.R" )  

f <- function ( x ) { ( sin ( x ) )^2 + log ( x ) }
cat ( "\n" )
cat ( "  f(x) = ( sin ( x ) )^2 + log ( x )\n" )
cat ( "\n" )

a = 1.0
b = 10.0
n = 4
value <- adaptint ( f, a, b, n = n )
cat ( "  adaptint ( f, ", a, ",", b, ", n = ", n, " ) = ", value, "\n" )

a = 1.0
b = 10.0
n = 5
value <- adaptint ( f, a, b, n = n )
cat ( "  adaptint ( f, ", a, ",", b, ", n = ", n, " ) = ", value, "\n" )

a = 1.0
b = 10.0
n = 6
value <- adaptint ( f, a, b, n = n )
cat ( "  adaptint ( f, ", a, ",", b, ", n = ", n, " ) = ", value, "\n" )
#
#  Terminate.
#
cat ( "\n" )
cat ( "adaptint_test\n" )
cat ( "  Normal end of execution.\n" )
cat ( date ( ), "\n" )

quit ( )