Skip to content

Instantly share code, notes, and snippets.

@eliardocosta

eliardocosta/pi_por_simulacao.R

Last active June 18, 2021 13:47
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save eliardocosta/94854893198814d62fe62470380bf675 to your computer and use it in GitHub Desktop.
Save eliardocosta/94854893198814d62fe62470380bf675 to your computer and use it in GitHub Desktop.
Download ZIP
Obtendo a constante Pi por simulação de pontos no (0,1)x(0,1) e calculando a proporção de pontos dentro da circunferência.
Raw
pi_por_simulacao.R
r <- 0.5 # raio da circunferencia
n <- 1E5
x <- runif(n)
y <- runif(n)
pdentro <- numeric(n)
for (i in seq_len(n)) { # contando os pontos dentro da circunferencia
pdentro[i] <- ifelse((x[i] - 0.5)^2 + (y[i] - 0.5)^2 <= r^2, 1, 0)
}
4*sum(pdentro)/n # estimador de Pi
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment