Task 4

# Task 4
# assume 9 number cards, 3 figure cards, 1 ace => 13 card types => probability = 1/13
p <- 1/13
cardDistr <- rgeom(n = 200, prob = p)

quantile(x = cardDistr, 0.05)
quantile(x = cardDistr, 0.95)

hist(cardDistr)
abline(v = quantile(x = cardDistr, 0.05), lwd = 3, lty = 3, col = "green")
abline(v = quantile(x = cardDistr, 0.95), lwd = 3, lty = 3, col = "green")

# the theoretical distribution of P(X = 10)
dgeom(x = 10, p = p)