2.2 #8

reps <- 10000

cars <- c(rep("S", 3), rep("D", 3), rep("I", 3))

SportsAdjacent <- \() {
  arrangement <- sample(cars, 9, replace = FALSE)
  sport_positions <- which(arrangement == "S")
  (max(sport_positions) - min(sport_positions)) == 2
}

sim2 <- replicate(reps, SportsAdjacent())

simProb2 <- mean(sim2)
simProb2

[1] 0.0849