simulate <- function(){
  aantal <- 1000
  means <- rep(NA,aantal)
  sds <- rep(NA,aantal)

  for (i in seq(1,aantal)) {
    result <- rep(NA,300)
    weken <- rmultinom(300, 1, c(dhyper(0,6,39,6), dhyper(1,6,39,6), dhyper(2,6,39,6), dhyper(3,6,39,6), dhyper(4,6,39,6), dhyper(5,6,39,6), dhyper(6,6,39,6)))

    winst <- 0
    for (i in seq(1,300)){
      winst <- weken[[3, i]] * 1.25 + weken[[4, i]] * 12.5 + weken[[5, i]] * 125 + weken[[6, i]] * 5000 + weken[[7, i]] * 100000
      result[i] <- winst
    }
    sds[i] <- sd(result)
    means[i] <- mean(result)

  }
  View(means)

  gem <- mean(means)
  sdd <- mean(sds)

  print(gem)
  print(sdd)


}

simulate()