Data605discussion10

M = matrix(c(1, 0, 0,
             .3, .4, .3,
             .2, .1, .7),
           nrow = 3, byrow = T
)
colnames(M) = c("H", "Y", "D")
rownames(M) = c("H", "Y", "D")

M %^% 2

##      H    Y    D
## H 1.00 0.00 0.00
## Y 0.48 0.19 0.33
## D 0.37 0.11 0.52

## Probability is 1