N <- 1000000
X <- sample.int(6, N, TRUE)
mean(X)

X <- runif(N, -1, +1)
Y <- runif(N, -1, +1)
Z <- sqrt(X^2+Y^2) <= 1
4*sum(Z)/N