R Code for Discrete Probability Example 1

result <- rep(0,10000)# 10000 simulated rolls
for(i in 1:10000) result[i] <- sample(1:6,1)+sample(1:6,1)
table(result)
prob1 <- table(result)/ # simulated probabilities
prob1
outcomes<-c(2:12)
prob2<-rep(0,11)
for(i in 1:6) for(j in 1:6){
prob2[outcomes==(i+j)] <- prob2[outcomes==(i+j)]+(1/6)*(1/6)}
prob2