Skip to content

M
mooc-rr
Commit 413d5b1a authored by fce93069b4bce444fe8565ce6d49f908's avatar fce93069b4bce444fe8565ce6d49f908
Browse files
Options

Update toy_document_fr.Rmd

parent b7b7820a
Hide whitespace changes
Inline Side-by-side
Showing with 5 additions and 3 deletions
module2/exo1/toy_document_fr.Rmd
View file @ 413d5b1a
...@@ -9,6 +9,7 @@ output: html_document ...@@ -9,6 +9,7 @@ output: html_document
```{r setup, include=FALSE} ```{r setup, include=FALSE}
knitr::opts_chunk$set(echo = TRUE) knitr::opts_chunk$set(echo = TRUE)
``` ```
## En demandant à la lib maths ## En demandant à la lib maths
Mon ordinateur m'indique que $\pi$ vaut *approximativement* Mon ordinateur m'indique que $\pi$ vaut *approximativement*
...@@ -22,7 +23,7 @@ Mais calculé avec la __méthode__ des [aiguilles de Buffon](https://fr.wikipedi ...@@ -22,7 +23,7 @@ Mais calculé avec la __méthode__ des [aiguilles de Buffon](https://fr.wikipedi
```{r} ```{r}
set.seed(42) set.seed(42)
N=100000 N=100000
x =runif(N) x = runif(N)
theta = pi/2*runif(N) theta = pi/2*runif(N)
2/(mean(x+sin(theta)>1)) 2/(mean(x+sin(theta)>1))
``` ```
...@@ -32,7 +33,7 @@ Sinon, une méthode plus simple à comprendre et ne faisant pas intervenir d'app ...@@ -32,7 +33,7 @@ Sinon, une méthode plus simple à comprendre et ne faisant pas intervenir d'app
```{r} ```{r}
set.seed(42) set.seed(42)
N = 1000 N=1000
df = data.frame(X = runif(N), Y = runif(N)) df = data.frame(X = runif(N), Y = runif(N))
df$Accept = (df$X**2 + df$Y**2 <=1) df$Accept = (df$X**2 + df$Y**2 <=1)
library(ggplot2) library(ggplot2)
...@@ -44,4 +45,5 @@ Il est alors aisé d'obtenir une approximation (pas terrible) de $\pi$ en compta ...@@ -44,4 +45,5 @@ Il est alors aisé d'obtenir une approximation (pas terrible) de $\pi$ en compta
```{r} ```{r}
4*mean(df$Accept) 4*mean(df$Accept)
``` ```
\ No newline at end of file
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment