Buffon's needle problem

The problem in more mathematical terms is: Given a needle of length ${\displaystyle l}$ dropped on a plane ruled with parallel lines t units apart, what is the probability that the needle will lie across a line upon landing?

Let x be the distance from the center of the needle to the closest parallel line, and let θ be the acute angle between the needle and one of the parallel lines.

The uniform probability density function of x between 0 and t /2 is

${\displaystyle {\begin{cases}{\frac {2}{t}}&:\ 0\leq x\leq {\frac {t}{2}}\\0&:{\text{elsewhere.}}\end{cases}}}$

Here, x = 0 represents a needle that is centered directly on a line, and x = t/2 represents a needle that is perfectly centered between two lines. The uniform PDF assumes the needle is equally likely to fall anywhere in this range, but could not fall outside of it.

The uniform probability density function of θ between 0 and π/2 is

${\displaystyle {\begin{cases}{\frac {2}{\pi }}&:\ 0\leq \theta \leq {\frac {\pi }{2}}\\0&:{\text{elsewhere.}}\end{cases}}}$

Here, θ = 0 radians represents a needle that is parallel to the marked lines, and θ = π/2 radians represents a needle that is perpendicular to the marked lines. Any angle within this range is assumed an equally likely outcome.

The two random variables, x and θ, are independent, so the joint probability density function is the product

${\displaystyle {\begin{cases}{\frac {4}{t\pi }}&:\ 0\leq x\leq {\frac {t}{2}},\ 0\leq \theta \leq {\frac {\pi }{2}}\\0&:{\text{elsewhere.}}\end{cases}}}$

The needle crosses a line if

${\displaystyle x\leq {\frac {l}{2}}\sin \theta .}$

Now there are two cases.

The following solution for the "short needle" case, while equivalent to the one above, has a more visual flavor, and avoids iterated integrals.

We can calculate the probability ${\displaystyle P}$ as the product of 2 probabilities: ${\displaystyle P=P_{1}\cdot P_{2}}$, where ${\displaystyle P_{1}}$ is the probability that the center of the needle falls close enough to a line for the needle to possibly cross it, and ${\displaystyle P_{2}}$ is the probability that the needle actually crosses the line, given that the center is within reach.

Looking at the illustration in the above section, it is apparent that the needle can cross a line if the center of the needle is within ${\displaystyle l/2}$ units of either side of the strip. Adding ${\displaystyle {\frac {l}{2}}+{\frac {l}{2}}}$ from both sides and dividing by the whole width ${\displaystyle t}$, we obtain ${\displaystyle P_{1}={\frac {l}{t}}.}$

The red and blue needles are both centered at x. The red one falls within the gray area, contained by an angle of 2θ on each side, so it crosses the vertical line; the blue one does not. The proportion of the circle that is gray is what we integrate as the center x goes from 0 to 1

Now, we assume that the center is within reach of the edge of the strip, and calculate ${\displaystyle P_{2}}$. To simplify the calculation, we can assume that ${\displaystyle l=2}$.

Let x and θ be as in the illustration in this section. Placing a needle's center at x, the needle will cross the vertical axis if it falls within a range of 2θ radians, out of π radians of possible orientations. This represents the gray area to the left of x in the figure. For a fixed x, we can express θ as a function of x: ${\displaystyle \theta \left(x\right)=\cos ^{-1}\left(x\right)}$. Now we can let x move from 0 to 1, and integrate:

${\displaystyle P_{2}=\int _{0}^{1}{\frac {2\theta (x)}{\pi }}\,dx={\frac {2}{\pi }}\int _{0}^{1}\cos ^{-1}(x)\,dx={\frac {2}{\pi }}\cdot 1={\frac {2}{\pi }}.}$

Multiplying both results, we obtain ${\displaystyle P=P_{1}\cdot P_{2}={\frac {l}{t}}{\frac {2}{\pi }}={\frac {2l}{t\pi }}}$, as above.

There is an even more elegant and simple method of calculating the "short needle case". The end of the needle farthest away from any one of the two lines bordering its region must be located within a horizontal (perpendicular to the bordering lines) distance of ${\displaystyle l\cos \theta }$ (where ${\displaystyle \theta }$ is the angle between the needle and the horizontal) from this line in order for the needle to cross it. The farthest this end of the needle can move away from this line horizontally in its region is ${\displaystyle t}$. The probability that the farthest end of the needle is located no more than a distance ${\displaystyle l\cos \theta }$ away from the line (and thus that the needle crosses the line) out of the total distance ${\displaystyle t}$ it can move in its region for ${\displaystyle 0\leq \theta \leq \pi /2}$ is given by

${\displaystyle P={\frac {\int _{0}^{\frac {\pi }{2}}l\cos \theta d\theta }{\int _{0}^{\frac {\pi }{2}}td\theta }}={\frac {l}{t}}{\frac {\int _{0}^{\frac {\pi }{2}}\cos \theta d\theta }{\int _{0}^{\frac {\pi }{2}}d\theta }}={\frac {l}{t}}{\frac {1}{\frac {\pi }{2}}}={\frac {2l}{t\pi }}}$, as above.

The short-needle problem can also be solved without any integration, in a way that explains the formula for p from the geometric fact that a circle of diameter t will cross the distance t strips always (i.e. with probability 1) in exactly two spots. This solution was given by Joseph-Émile Barbier in 1860[3] and is also referred to as "Buffon's noodle".

An experiment to find π. Matches with the length of 9 squares have been thrown 17 times between rows with the width of 9 squares. 11 of the matches have landed at random across the drawn lines marked by the green points.
(2l·n)/(th)=(2·9·17)/(9·11)≈3.1≈π.

A Python 3 based simulation using Matplotlib to sketch Buffon's needle experiment with the parameters t = 5.0, l = 2.6. Observe the calculated value of Pi (y-axis) approaching 3.14 as the number of tosses (x-axis) approaches infinity.

In the first, simpler case above, the formula obtained for the probability ${\displaystyle P}$ can be rearranged to:

$${\displaystyle \pi ={\frac {2l}{tP}}}$$

Thus, if we conduct an experiment to estimate ${\displaystyle P}$, we will also have an estimate for π.

Suppose we drop n needles and find that h of those needles are crossing lines, so ${\displaystyle P}$ is approximated by the fraction ${\displaystyle h/n}$. This leads to the formula:

${\displaystyle \pi \approx {\frac {2l\cdot n}{th}}.}$

In 1901, Italian mathematician Mario Lazzarini performed Buffon's needle experiment. Tossing a needle 3408 times, he obtained the well-known approximation 355/113 for π, accurate to six significant digits.[4] Lazzarini's "experiment" is an example of confirmation bias, as it was set up to replicate the already well-known approximation of 355/113 (in fact, there is no better rational approximation with fewer than five digits in the numerator and denominator), yielding a more accurate "prediction" of π than would be expected from the number of trials, as follows: [5]

Lazzarini chose needles whose length was 5/6 of the width of the strips of wood. In this case, the probability that the needles will cross the lines is ${\textstyle {\frac {5}{3\pi }}}$. Thus if one were to drop n needles and get x crossings, one would estimate π as:

$${\displaystyle \pi \approx {\frac {5}{3}}\cdot {\frac {n}{x}}}$$

So if Lazzarini was aiming for the result 355/113, he needed n and x such that:

$${\displaystyle {\frac {355}{113}}={\frac {5}{3}}\cdot {\frac {n}{x}}}$$

or equivalently,

$${\displaystyle x={\frac {113n}{213}}.}$$

To do this, one should pick n as a multiple of 213, because then ${\textstyle {\frac {113n}{213}}}$ is an integer; one then drops n needles, and hopes for exactly ${\textstyle x={\frac {113n}{213}}}$ successes.

If one drops 213 needles and happens to get 113 successes, then one can triumphantly report an estimate of π accurate to six decimal places. If not, one can just do 213 more trials and hope for a total of 226 successes; if not, just repeat as necessary. Lazzarini performed 3408 = 213 · 16 trials, making it seem likely that this is the strategy he used to obtain his "estimate."

The above description of strategy might even be considered charitable to Lazzarini. A statistical analysis of intermediate results he reported for fewer tosses leads to a very low probability of achieving such close agreement to the expected value all through the experiment. This makes it very possible that the "experiment" itself was never physically performed, but based on numbers concocted from imagination to match statistical expectations, but too well, as it turns out.[5]

Dutch science journalist Hans van Maanen argues, however, that Lazzarini's article was never meant to be taken too seriously as it would have been pretty obvious for the readers of the magazine (aimed at school teachers) that the apparatus that Lazzarini said to have built cannot possibly work as described.[6]

• Buffon's noodle
• Bertrand paradox (probability)

• Badger, Lee (April 1994). "Lazzarini's Lucky Approximation of π". Mathematics Magazine. Mathematical Association of America. 67 (2): 83–91. doi:10.2307/2690682. JSTOR 2690682.
• Ramaley, J. F. (October 1969). "Buffon's Noodle Problem". The American Mathematical Monthly. Mathematical Association of America. 76 (8): 916–918. doi:10.2307/2317945. JSTOR 2317945.
• Mathai, A. M. (1999). An Introduction to Geometrical Probability. Newark: Gordon & Breach. p. 5. ISBN 978-90-5699-681-9.
• Dell, Zachary; Franklin, Scott V. (September 2009). "The Buffon-Laplace needle problem in three dimensions". Journal of Statistical Mechanics: Theory and Experiment. 2009 (9): 010. Bibcode:2009JSMTE..09..010D. doi:10.1088/1742-5468/2009/09/P09010.
• Schroeder, L. (1974). "Buffon's needle problem: An exciting application of many mathematical concepts". Mathematics Teacher, 67 (2), 183–6.

Wikimedia Commons has media related to Buffon's needle.

• Buffon's Needle Problem at cut-the-knot
• Math Surprises: Buffon's Noodle at cut-the-knot
• MSTE: Buffon's Needle
• Buffon's Needle Java Applet
• Estimating PI Visualization (Flash)
• Buffon's needle: fun and fundamentals (presentation) at slideshare
• Animations for the Simulation of Buffon's Needle by Yihui Xie using the R package animation
• 3D Physical Animation by Jeffrey Ventrella
• Padilla, Tony. "∏ Pi and Buffon's Needle". Numberphile. Brady Haran. Archived from the original on 2013-05-17. Retrieved 2013-04-09.