Chapter 2
Monte Carlo Integration
This chapter gives an introduction to Monte Carlo integration. The main goals are to review some basic concepts of probability theory, to de®ne the notation and terminology that we will be using, and to summarize the variance reduction techniques that have proven most useful in computer graphics. Good references on Monte Carlo methods include Kalos & Whitlock [1986], Hammer- sley & Handscomb [1964], and Rubinstein [1981]. Sobol' [1994] is a good starting point for those with little background in probability and statistics. Spanier & Gelbard [1969] is the classic reference for Monte Carlo applications to neutron transport problems; Lewis & Miller [1984] is a good source of background information in this area. For quasi-Monte Carlo methods, see Niederreiter [1992], Beck & Chen [1987], and Kuipers & Niederreiter [1974].
2.1 A brief history
Monte Carlo methods originated at the Los Alamos National Laboratory in the early years after World War II. The ®rst electronic computer in the United States had just been com- pleted (the ENIAC), and the scientists at Los Alamos were considering how to use it for the design of thermonuclear weapons (the H-bomb). In late 1946 Stanislaw Ulam suggested the use of random sampling to simulate the ¯ight paths of neutrons, and John von Neumann
29 30 CHAPTER 2. MONTE CARLO INTEGRATION
developed a detailed proposal in early 1947. This led to small-scale simulations whose re- sults were indispensable in completing the project. Metropolis & Ulam [1949] published a paper in 1949 describing their ideas, which sparked to a great deal of research in the 1950's [Meyer 1956]. The name of the Monte Carlo method comes from a city in Monaco, famous for its casinos (as suggested by Nick Metropolis, another Monte Carlo pioneer). In isolated instances, random sampling had been used much earlier to solve numerical problems [Kalos & Whitlock 1986]. For example, in 1777 the Comte de Buffon performed
an experiment in which a needle was dropped many times onto a board marked with equidis-
d L tant parallel lines. Letting L be the length of the needle and be the distance between
the lines, he showed that the probability of the needle intersecting a line is