# -----------------------------------------------------------------------
# G I O C O   D E L   C R A P S
# Il Banco lancia due dadi e si fa la somma dei punti.
# Il Giocatore vince se esce 7 o 11. Perde se esce 2 o 3 o 12.
# Se esce un altro numero, questo è detto PUNTO.
# Il gioco continua finchè‚ il giocatore vince se esce il PUNTO
# o perde se esce il 7. Che probabilità di vincere ha il giocatore?
# La probabilità teorica è circa 0.493

dado <- c(1:6)

prove <- 100000
num <- 0; den <- 0
num <- rep(NA,prove); den <- num

load("rs_craps_01.rdata")

runif(1)
.Random.seed <- rs

# -----------------------------------------------------------------------
for (i in 1:prove) {
    out <- sample(dado,size=2,replace=TRUE)
    den[i] <- 1
    x <- sum(out)
    
    swt <- 3
    if (x==7 | x==11)        swt <- 1
    if (x==2 | x==3 | x==12) swt <- 2
    
    switch(swt,
           num[i] <- +1,
           num[i] <- -1,
           {punto <- x
            repeat {
                    out <- sample(dado,size=2,replace=TRUE)
                    x <- sum(out)
                    if (x==punto) {
                       num[i] <- +1
                       break
                    }
                    if (x==7) {
                       num[i] <- -1
                       break
                    }   
            }
           }
    )
}

x <- cumsum(den)
y <- cumsum(num)

plot(x,y,log="x",type="l",xlab="Logaritmo del numero della partita",ylab="Guadagno totale")   # guadagno "totale" fino alla puntata x-esima
abline(h=0,lty=2)            

win.graph()
plot(x,y/x,log="x",type="l",xlab="Logaritmo del numero della partita",ylab="Guadagno medio") # guadagno medio ...
abline(h=0,lty=2)            
# -----------------------------------------------------------------------