# 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