cmna (version 0.1.2)

mcint: Monte Carlo Integration

Description

Simple Monte Carlo Integraton

Usage

mcint(f, a, b, m = 1000, max.y = 1)
mcint2(f, xdom, ydom, m = 1000, max.z = 1)

Arguments

f

function to integrate

a

the lower-bound of integration

b

the upper-bound of integration

m

the number of subintervals to calculate

max.y

the largest expected value of the range

xdom

the domain on x of integration in two dimensions

ydom

the domain on y of integration in two dimensions

max.z

the largest expected value of the range in two dimensions

Value

the value of the integral

Details

The mcint function uses a simple Monte Carlo algorithm to estimate the value of an integral. The parameter n sets the total number of evaluation points. The parameter max.y is the maximum expected value of the range of function f. The mcint2 provides Monte Carlo integration in two dimensions.

See Also

Other integration: adaptint, gaussint, giniquintile, midpt, revolution-solid, romberg, simp38, simp, trap

Examples

f <- function(x) { sin(x)^2 + log(x)}
mcint(f, 0, 1)
mcint(f, 0, 1, m = 10e6)