1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
|
% Generated by roxygen2: do not edit by hand
% Please edit documentation in R/buffon.needle.R
\name{buffon.needle}
\alias{buffon.needle}
\title{Simulation of Buffon's Needle}
\usage{
buffon.needle(
l = 0.8,
d = 1,
redraw = TRUE,
mat = matrix(c(1, 3, 2, 3), 2),
heights = c(3, 2),
col = c("lightgray", "red", "gray", "red", "blue", "black", "red"),
expand = 0.4,
type = "l",
...
)
}
\arguments{
\item{l}{numerical. length of the needle; shorter than \code{d}.}
\item{d}{numerical. distances between lines; it should be longer than
\code{l}.}
\item{redraw}{logical. redraw former `needles' or not for each drop.}
\item{mat, heights}{arguments passed to \code{\link{layout}} to set the layout
of the three graphs.}
\item{col}{a character vector of length 7 specifying the colors of:
background of the area between parallel lines, the needles, the sin curve,
points below / above the sin curve, estimated \eqn{\pi} values, and the
true \eqn{\pi} value.}
\item{expand}{a numerical value defining the expanding range of the y-axis
when plotting the estimated \eqn{\pi} values: the \code{ylim} will be
\code{(1 +/- expand) * pi}.}
\item{type}{an argument passed to \code{\link[graphics:plot.default]{plot}}
when plotting the estimated \eqn{\pi} values (default to be lines).}
\item{\dots}{other arguments passed to
\code{\link[graphics:plot.default]{plot}} when plotting the values of
estimated \eqn{\pi}.}
}
\value{
The values of estimated \eqn{\pi} are returned as a numerical vector
(of length \code{nmax}).
}
\description{
This function provides a simulation for the problem of Buffon's Needle, which
is one of the oldest problems in the field of geometrical probability.
}
\details{
This is quite an old problem in probability. For mathematical background,
please refer to \url{https://en.wikipedia.org/wiki/Buffon's_needle} or
\url{https://mste.illinois.edu/activity/buffon/}.
`Needles' are denoted by segments on the 2D plane, and dropped randomly to
check whether they cross the parallel lines. Through many times of `dropping'
needles, the approximate value of \eqn{\pi} can be calculated out.
There are three graphs made in each step: the top-left one is a simulation of
the scenario, the top-right one is to help us understand the connection
between dropping needles and the mathematical method to estimate \eqn{\pi},
and the bottom one is the result for each drop.
}
\note{
Note that \code{redraw} has great influence on the speed of the
simulation (animation) if the control argument \code{nmax} (in
\code{\link{ani.options}}) is quite large, so you'd better specify it as
\code{FALSE} when doing a large amount of simulations.
The maximum number of drops is specified in \code{ani.options('nmax')}.
}
\references{
Examples at \url{https://yihui.org/animation/example/buffon-needle/}
Ramaley, J. F. (Oct 1969). Buffon's Noodle Problem. \emph{The
American Mathematical Monthly} \bold{76} (8): 916-918.
}
\author{
Yihui Xie
}
|
