--- title: "A propos de pi" author: "Berger Audric" date: "10/02/2022" output: html_document --- #Formule mathématique ```{r setup, include=FALSE} y=2 x=sin(y*y) print(x) ``` ```{r setup, include=FALSE} p=$\pi$ print(p) ``` #Lien aiguille de buffon [Aiguille de Buffon](https://fr.wikipedia.org/wiki/Aiguille_de_Buffon) ```{r} set.seed(42) N = 100000 x = runif(N) theta = pi/2*runif(N) R=2/(mean(x+sin(theta)>1)) plot=R library(ggplot2) ggplot(R) ``` >ggplot ne fonctionne pas problème avec FORTITUDE? ```{r} set.seed(42) N = 1000 df = data.frame(X = runif(N), Y = runif(N)) df$Accept = (df$X**2 + df$Y**2 <=1) library(ggplot2) ggplot(df, aes(x=X,y=Y,color=Accept)) + geom_point(alpha=.2) + coord_fixed() + theme_bw() ``` >Ce code fonctionne..