Origins and basic principle
The Monte Carlo method began in 1946 when Stanisław Ulam, a mathematician at Los Alamos recovering from illness, tried to work out the odds of winning a game of solitaire. If the calculation was intractable, he reasoned, one could simply play many times and count the successes. John von Neumann applied the idea to neutron diffusion, and Nicholas Metropolis named it after the casino Ulam's uncle frequented. Metropolis and Ulam (1949) published the first paper. The classic illustration estimates π: scatter points at random in a unit square, and four times the fraction landing inside the quarter circle approaches π. The average converges to the true value as samples grow, thanks to the law of large numbers.
Accuracy and variance reduction
By the central limit theorem, Monte Carlo error shrinks roughly in proportion to one over the square root of the number of samples. Ten times more accuracy requires a hundred times more samples, but this rate is nearly independent of dimension, so Monte Carlo beats grid-based numerical methods for high-dimensional integrals. Variance-reduction techniques such as importance sampling, control variates, stratified sampling (the same principle as stratified probability sampling in surveys) and quasi-Monte Carlo improve efficiency. Quality also depends on the pseudo-random number generator; Matsumoto and Nishimura's (1998) Mersenne Twister is widely used. Hammersley and Handscomb (1964), Robert and Casella (2004) and Owen (2013) are standard texts.
Applications
The algorithm of Metropolis et al. (1953) and Hastings's (1970) generalization grew into MCMC, which uses Markov chains to sample from complex distributions, and after Gelfand and Smith (1990) it became the standard computational tool of Bayesian statistics. Efron's (1979) bootstrap, which resamples data to assess the uncertainty of estimates, is also a Monte Carlo approach. In finance, Boyle (1977) introduced Monte Carlo to option pricing, and it is used in scenario analysis for risk management. Physics uses it for particle transport and statistical mechanics, and computer graphics for path-traced rendering. In AI, Monte Carlo tree search, rooted in Kocsis and Szepesvári's (2006) UCT algorithm, was central to the AlphaGo system of Silver et al. (2016), an important branch of AI search algorithms.
FAQ
- Is Monte Carlo simulation always accurate?
- No. Its results are probabilistic approximations whose error depends on sample size, random-number quality and model assumptions, which is why they are usually reported with confidence intervals or standard errors.
- How does the Monte Carlo method differ from MCMC?
- Plain Monte Carlo draws independent samples directly, while MCMC uses a Markov chain to generate a sequence of dependent samples from distributions that are hard to sample directly.
Sources reviewed: 2026-10-26
References
- Metropolis, N., & Ulam, S. (1949). The Monte Carlo method. Journal of the American Statistical Association, 44(247), 335–341.
- Metropolis, N., Rosenbluth, A. W., Rosenbluth, M. N., Teller, A. H., & Teller, E. (1953). Equation of state calculations by fast computing machines. Journal of Chemical Physics, 21(6), 1087–1092.
- Hastings, W. K. (1970). Monte Carlo sampling methods using Markov chains and their applications. Biometrika, 57(1), 97–109.
- Eckhardt, R. (1987). Stan Ulam, John von Neumann, and the Monte Carlo method. Los Alamos Science, 15, 131–137.
- Hammersley, J. M., & Handscomb, D. C. (1964). Monte Carlo Methods. London: Methuen.
- Boyle, P. P. (1977). Options: A Monte Carlo approach. Journal of Financial Economics, 4(3), 323–338.
- Efron, B. (1979). Bootstrap methods: Another look at the jackknife. Annals of Statistics, 7(1), 1–26.
- Gelfand, A. E., & Smith, A. F. M. (1990). Sampling-based approaches to calculating marginal densities. Journal of the American Statistical Association, 85(410), 398–409.
- Kocsis, L., & Szepesvári, C. (2006). Bandit based Monte-Carlo planning. In Machine Learning: ECML 2006 (pp. 282–293). Berlin: Springer.
- Silver, D., et al. (2016). Mastering the game of Go with deep neural networks and tree search. Nature, 529(7587), 484–489.
- Robert, C. P., & Casella, G. (2004). Monte Carlo Statistical Methods (2nd ed.). New York: Springer.
- Owen, A. B. (2013). Monte Carlo Theory, Methods and Examples. Stanford, CA: Stanford University.
- Matsumoto, M., & Nishimura, T. (1998). Mersenne twister: A 623-dimensionally equidistributed uniform pseudo-random number generator. ACM Transactions on Modeling and Computer Simulation, 8(1), 3–30.
