My computer tells me that $\pi$ is /approximatively/
#+begin_src R :results output :session *R* :exports both
pi
#+end_src
...
...
@@ -21,7 +21,8 @@ pi
: [1] 3.141593
* Buffon's needle
Applying the method of [[https://en.wikipedia.org/wiki/Buffon%2527s_needle_problem][Buffon's needle]] we get the *approximation*
Applying the method of [[https://en.wikipedia.org/wiki/Buffon%2527s_needle_problem][Buffon's needle]], we get the *approximation*
#+begin_src R :results output :session *R* :exports results
set.seed(42)
N = 100000
...
...
@@ -34,10 +35,7 @@ theta = pi/2*runif(N)
: [1] 3.14327
* Using a surface fraction argument
A method that is easier to understand and does not make use of the
$\sin$ function is based on the fact that if $X\simU(0,1)$ and
$Y\simU(0,1)$, then $P[X^2+Y^2 \le1]=\pi/4$ (see [[https://en.wikipedia.org/wiki/Monte_Carlo_method]["Monte Carlo method" on
Wikipedia]]). The following code uses this approach:
A method that is easier to understand and does not make use of the $\sin$ function is based on the fact that if $X\simU(0,1)$ and $Y\simU(0,1)$, then $P[X^2+Y^2 \le1]=\pi/4$ (see [[https://en.wikipedia.org/wiki/Monte_Carlo_method]["Monte Carlo method" on Wikipedia]]). The following code uses this approach:
#+begin_src R :results output graphics :file figure_pi_mc1.png :exports both :width 600 :height 400 :session *R*