# MATA271 Tietokoneharjoitus 1
# Lasse Leskelä 2012-09-24
# Simuloidaan useita otoksia Billin pelikassa 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+1);
Z <- numeric(Otokset);

for (i in 1:Otokset) {

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

  # Lasketaan pelaajan kassavirta otoksen funktiona.
  # V[t] on kassan arvo juuri ennen t:ttä pelikierrosta
  V[1] <- V0;
  for (t in 1:n) {
    if ((omega[t] >= 1) & (omega[t] <= 18)) {
      V[t+1] <- V[t] + 1;
    } else {
      V[t+1] <- V[t] - 1;
    }
  }

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

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);