dadi <- function(n) {u <- 0; for (i in 1:n) u <- u+floor(runif(1)*6)+1; u}
# ovvero:     dadi <- function(n) sum(floor(runif(n)*6)+1)
# ovvero:     dadi <- function(n) sum(floor(runif(n,min=1,max=7)))
# runif(1) genera un numero tra 0 ed 1 a caso, runif(2) ne genera 2, ...
# runif(..., min=h, max=k) li genera in [h,k)
s <- 0; tot <- 0; n <- 1000; for (i in 1:n) {tot<-tot+1; if (dadi(2)>7) s <- s+1}; tot; s/tot
[1] 1000
[1] 0.414
# aggiungo 1000 lanci, in modo da averne 2000
n <- 1000; for (i in 1:n) {tot<-tot+1; if (dadi(2)>7) s <- s+1}; tot; s/tot
[1] 2000
[1] 0.428
n <- 3000; for (i in 1:n) {tot<-tot+1; if (dadi(2)>7) s <- s+1}; tot; s/tot
[1] 5000
[1] 0.409
n <- 5000; for (i in 1:n) {tot<-tot+1; if (dadi(2)>7) s <- s+1}; tot; s/tot
[1] 10000
[1] 0.4158
n <- 10000; for (i in 1:n) {tot<-tot+1; if (dadi(2)>7) s <- s+1}; tot; s/tot
[1] 20000
[1] 0.4173
n <- 30000; for (i in 1:n) {tot<-tot+1; if (dadi(2)>7) s <- s+1}; tot; s/tot
[1] 50000
[1] 0.42048
n <- 50000; for (i in 1:n) {tot<-tot+1; if (dadi(2)>7) s <- s+1}; tot; s/tot
[1] 1e+05
[1] 0.41842
n <- 100000; for (i in 1:n) {tot<-tot+1; if (dadi(2)>7) s <- s+1}; tot; s/tot
[1] 2e+05
[1] 0.417625
n <- 200000; for (i in 1:n) {tot<-tot+1; if (dadi(2)>7) s <- s+1}; tot; s/tot
[1] 4e+05
[1] 0.417535
n <- 600000; for (i in 1:n) {tot<-tot+1; if (dadi(2)>7) s <- s+1}; tot; s/tot
[1] 1e+06
[1] 0.416217