###################################################
### code chunk number 44: Cs40_Vinfty
###################################################
m <- nrow(B)
vecV <- solve(diag(m * m) - kronecker(B, B)) %*% as.vector(Q)
V_inf <- matrix(vecV, nrow = m, ncol = m)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.