# Section 4; Question 1:
#
n_mc = 100000
data = rep( NA, n_mc ) 
for( ii in 1:n_mc ){
    d = sum( sample( 1:5, 2 ) )
    data[ii] = d
}
P = table(data)/n_mc
print(P)
print(sum(P[3:5]))