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()


---