# MATA280 Tietokoneharjoitus
# Lasse Leskelä 2013-11-17
# Simuloidaan useita otoksia Valtterin pelikassan päätösarvosta.

V0 <- 300;       # Alkupääoma
n  <- 300;       # Pelikierrosten lkm
Otokset <- 500;  # Otosten lkm

# Alustetaan vektorit (eli luodaan sopivan kokoiset nollavektorit).
V <- numeric(n);
Z <- numeric(Otokset);

for (i in 1:Otokset) {

  # Arvotaan i:s otos.
  omega <- sample(0:36, n, replace=TRUE);

  # Lasketaan i:nnettä otosta vastaava pelikassa ajan funktiona.
  V <- numeric(n)
  V_after <- V0
  for (t in 1:n) {
    V_before <- V_after
    if ((omega[t] >= 1) & (omega[t] <= 18)) {
      V_after <- V_before + 1
    } else {
      V_after <- V_before - 1
    }
    V[t] <- V_after
  }

  # Tallennetaan pelikassan päätösarvo i:nnessä otoksessa.
  Z[i] <- V[n];
}

mean_emp <- mean(Z)
pr_win_10_emp  <- length(which(Z>1.1*V0))/length(Z);
pr_lose_10_emp <- length(which(Z<.9*V0))/length(Z);

hist(Z);