# Studio della prob. che a "7 e 1/2", avendo "1 quadri" e "3 fiori", # pescata una carta si sballi, ossia si estragga una carta che valga 4 o più # Numeriamo le carte secondo: Come Quando Fuori Piove (da 1 a 10 cuori, da 11 a # 20 quadri, da 21 a 30 fiori, da 31 a 40 picche); 1Q e 3F sono le carte 11 e 33. # Vengono fatte 10 mila prove. Puoi ripeterle o cambiare numero di prove. # Il programma "test" genera la distribuzione di n uscite. test <- function(n) { no <- 0; prove <- n; carta <- vector(length=40); carta1 <- vector(length=40) for (k in 1:prove) { for (i in 1:40) carta[i] <- 1; carta[11] <- 0; carta[33] <- 0; ko <- 0; k <- 0; while(k==0) {c <- floor(runif(1)*40)+1; if(carta[c]==1) k <- 1}; carta1[c] <- carta1[c]+1; if((3 < c & c < 8)|(13 < c & c < 18)|(23 < c & c < 28)|(33 < c & c < 38)) ko <- 1; no <- no+ko } no/prove*100 carta2 <- vector(length=8) for (c in 1:40) {k <- c-trunc(c/10)*10; if (k > 7 | k==0) k <- 8; carta2[k] <- carta2[k]+carta1[c]} carta2/prove*100 } # # Un esempio di uscite: U=test(10000); BarNames = c("1","2","3","4","5","6","7","F"); Bar(U); sum(U[4:7]) # yellow,cyan,... % 8.07 10.84 8.17 10.76 10.31 10.24 10.88 30.73 # [1] 42.19 # Il valore teorico: 16/38 # [1] 0.4210526