P = matrix(c(.5, .25, .25,
.5, .25, .25,
.5, .25, .25), nrow = 3, ncol = 3, byrow=T)
P
## [,1] [,2] [,3]
## [1,] 0.5 0.25 0.25
## [2,] 0.5 0.25 0.25
## [3,] 0.5 0.25 0.25
P2 = P %*% P
P2
## [,1] [,2] [,3]
## [1,] 0.5 0.25 0.25
## [2,] 0.5 0.25 0.25
## [3,] 0.5 0.25 0.25
P3 = P2 %*% P
P3
## [,1] [,2] [,3]
## [1,] 0.5 0.25 0.25
## [2,] 0.5 0.25 0.25
## [3,] 0.5 0.25 0.25
Suppose n = 4 then \(P^n\):
P4 = P3 %*% P
P4
## [,1] [,2] [,3]
## [1,] 0.5 0.25 0.25
## [2,] 0.5 0.25 0.25
## [3,] 0.5 0.25 0.25