# 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.
#
prova <- 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]}
  barplot(carta2/prove*100,names.arg=c(1,2,3,4,5,6,7,"F")) }
#
prova(10000)  # Se vuoi copia questa riga