Nouveau2 toy_document_fr.Rmd

module2/exo1/toy_document_fr.Rmd

@@ -17,6 +17,33 @@ Mon ordinateur m’indique que π vaut approximativement :
 pi
 #> [1] 3.141593
 
+**En utilisant la méthode des aiguilles de Buffon**
+set.seed(42)
+N <- 100000
+x <- runif(N)
+theta <- pi/2 * runif(N)
+2 / (mean(x + sin(theta) > 1))
+#> [1] 3.14327
+
+**Avec un argument « fréquentiel » de surface (Méthode de Monte Carlo)**
+set.seed(42)
+N  <- 1000
+df <- data.frame(X = runif(N), Y = runif(N))
+df$Accept <- (df$X^2 + df$Y^2 <= 1)
+
+# Visualisation
+library(ggplot2)
+ggplot(df, aes(x = X, y = Y, color = Accept)) +
+  geom_point(alpha = .2) +
+  coord_fixed() +
+  theme_bw()
+
+**d’obtenir une approximation (pas terrible) de π**
+4 * mean(df$Accept)
+#> [1] 3.156
+
+
+
 
 
 ```{r setup, include=FALSE}