diff --git a/module2/exo1/toy_document_fr.Rmd b/module2/exo1/toy_document_fr.Rmd index 2614355e9225700b1c143e110b5d52b9a5d16f9b..6cc552f23728df5b27822becac482c110580716e 100644 --- a/module2/exo1/toy_document_fr.Rmd +++ b/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}