9887

Adaptive Monte Carlo Integration

This Demonstration compares adaptive and nonadaptive Monte Carlo integration for two different functions, and . The plot shows the places on the interval where sample points are added as the number of sample points is increased. The actual values of the integrals to six significant figures are 2.66667 and 5.65685. The adaptive technique generally gets better estimates with the same number of sample points by subdividing the subinterval with the highest error estimate. Normally this process would be repeated until some error criterion is satisfied, but in this Demonstration only four subdivisions are shown.

SNAPSHOTS

  • [Snapshot]
  • [Snapshot]
  • [Snapshot]

DETAILS

In one dimension, if is integrable on the interval , then the average value of is . From this, the nonadaptive (basic) Monte Carlo method estimates the integral by , where is the mean of the function evaluated at randomly sampled points. The corresponding error estimate is , where is the variance of the function evaluations. Thus the error can be reduced by increasing the number of sample points or reducing the variance. The key idea of adaptive Monte Carlo is to reduce the variance by recursive subdivision. Parts of this Demonstration are modifications of examples from the Wolfram Mathematica Documentation Center tutorial "Advanced Numerical Integration in Mathematica" (see link below). Students should try to explain why Monte Carlo integration has more difficulty for the function with the singularity at 0. (Hint: look at the formulas for the integral estimate and error estimate.)
    • Share:

Embed Interactive Demonstration New!

Just copy and paste this snippet of JavaScript code into your website or blog to put the live Demonstration on your site. More details »

Files require Wolfram CDF Player or Mathematica.









 
RELATED RESOURCES
Mathematica »
The #1 tool for creating Demonstrations
and anything technical.
Wolfram|Alpha »
Explore anything with the first
computational knowledge engine.
MathWorld »
The web's most extensive
mathematics resource.
Course Assistant Apps »
An app for every course—
right in the palm of your hand.
Wolfram Blog »
Read our views on math,
science, and technology.
Computable Document Format »
The format that makes Demonstrations
(and any information) easy to share and
interact with.
STEM Initiative »
Programs & resources for
educators, schools & students.
Computerbasedmath.org »
Join the initiative for modernizing
math education.
Step-by-step Solutions »
Walk through homework problems one step at a time, with hints to help along the way.
Wolfram Problem Generator »
Unlimited random practice problems and answers with built-in Step-by-step solutions. Practice online or make a printable study sheet.
Wolfram Language »
Knowledge-based programming for everyone.
Powered by Wolfram Mathematica © 2014 Wolfram Demonstrations Project & Contributors  |  Terms of Use  |  Privacy Policy  |  RSS Give us your feedback
Note: To run this Demonstration you need Mathematica 7+ or the free Mathematica Player 7EX
Download or upgrade to Mathematica Player 7EX
I already have Mathematica Player or Mathematica 7+