#*****************************************************************************80
#
## mcint2_test tests mcint2.
#
#  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:
#
#    09 March 2020
#
#  Author:
#
#    John Burkardt
#
cat ( date ( ), "\n" )
cat ( "\n" )
cat ( "mcint2_test\n" )
cat ( " ", version$version.string, "\n" )
cat ( "  mcint2 estimates an integral over a rectangle using a Monte Carlo approach.\n" )

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

f <- function ( x, y ) { x^2 * y }
cat ( "\n" )
cat ( "  f(x) = x^2 * y\n" )
cat ( "\n" )

xdom = c ( 0.0, 1.0 )
ydom = c ( 0.0, 1.0 )
m = 1000
value <- mcint2 ( f, xdom, ydom, m = m )
cat ( "  mcint2 ( f, c(0,1), c(0,0), m = ", m, " ) = ", value, "\n" )

xdom = c ( 0.0, 1.0 )
ydom = c ( 0.0, 1.0 )
m = 10000
value <- mcint2 ( f, xdom, ydom, m = m )
cat ( "  mcint2 ( f, c(0,1), c(0,0), m = ", m, " ) = ", value, "\n" )

xdom = c ( 0.0, 1.0 )
ydom = c ( 0.0, 1.0 )
m = 100000
value <- mcint2 ( f, xdom, ydom, m = m )
cat ( "  mcint2 ( f, c(0,1), c(0,0), m = ", m, " ) = ", value, "\n" )

xdom = c ( 0.0, 1.0 )
ydom = c ( 0.0, 1.0 )
m = 1000000
value <- mcint2 ( f, xdom, ydom, m = m )
cat ( "  mcint2 ( f, c(0,1), c(0,0), m = ", m, " ) = ", value, "\n" )
#
#  Terminate.
#
cat ( "\n" )
cat ( "mcint2_test\n" )
cat ( "  Normal end of execution.\n" )
cat ( date ( ), "\n" )

quit ( )