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