diff --git a/Rstudio b/Rstudio new file mode 100644 index 0000000000000000000000000000000000000000..50ba6d9b578ed79caaa9d926568e117e3315f632 --- /dev/null +++ b/Rstudio @@ -0,0 +1,27 @@ +title: "Estimations de π avec R" +author: "Votre nom" +date: "`r format(Sys.Date(), '%d %B %Y')`" +output: html_document +--- + +```{r setup, include=FALSE} +knitr::opts_chunk$set(echo = TRUE, message = FALSE, warning = FALSE) +library(ggplot2) +piset.seed(42) +N <- 100000 +x <- runif(N) +theta <- pi / 2 * runif(N) +pi_buffon <- 2 / mean(x + sin(theta) > 1) +pi_buffon +set.seed(42) +N <- 1000 +df <- data.frame(X = runif(N), Y = runif(N)) +df$Accept <- (df$X^2 + df$Y^2 <= 1) +ggplot(df, aes(x = X, y = Y, color = Accept)) + + geom_point(alpha = 0.2) + + coord_fixed() + + theme_bw() + + +--- +