> d <- 0.85
> m <- matrix(c(0,1,1,0,0,1,1,0,0),3)
> pr <- rep(1,3)
> m
     [,1] [,2] [,3]
[1,]    0    0    1
[2,]    1    0    0
[3,]    1    1    0
> pr
[1] 1 1 1
> (pr <- (1-d)/3 + (d * (m %*% (pr/colSums(m)))))
      [,1]
[1,] 0.900
[2,] 0.475
[3,] 1.325
> (pr <- (1-d)/3 + (d * (m %*% (pr/colSums(m)))))
        [,1]
[1,] 1.17625
[2,] 0.43250
[3,] 0.83625
> (pr <- (1-d)/3 + (d * (m %*% (pr/colSums(m)))))
          [,1]
[1,] 0.7608125
[2,] 0.5499063
[3,] 0.9175312
>