From OrganicDesign Wiki
# Simulation 1
out <- c()
for(i in 1:1000) {
x <- sample(1:100, 10000, replace=TRUE)
out[i] <- mean(x)
}
hist(out)
# Empirical
summary(out)
# Expect
(1 + 100)/2
# Simulation 2
# rbeta generation
x <- round((rbeta(100000, 3,3))*100 - 1)
table(x)
table(x%%10)
hist(x%%10)
[[:Talk:I Ching / Divination/R]]