Let $$P=\begin{bmatrix}\frac{1}{3}&0&\frac{2}{3}\\\frac{1}{3}&\frac{2}{3}&0\\\frac{1}{3}&\frac{1}{3}&\frac{1}{3}\end{bmatrix}$$ find $p_{11}^{(n)},p_{12}^{(n)},p_{13}^{(n)}$
Since the characteristic function is $$-\mu^3+\frac{4}{3}\mu^2-\frac{1}{3}\mu=0$$ the eigenvalues are $$\mu_0=1,\mu_1=\frac{1}{3},\mu_2=0$$ Here the problem starts, from what I understood from the book, we need to found three constants such that $$p_{ij}^{(n)}=\mu_0^nA+\mu_1^n+\mu_2^nC=A+\left(\frac{1}{3}\right)^nB$$
But how can I find the values of A and B?
Hmm. I don't know how you can directly find $A$, $B$, $C$ but I can suggest a way to go about the problem. Now that you have found the eigenvalues $\mu_i$, solve for the corresponding normalized eigenvectors $e_i$. Then $$P^n = E^{-1}\Lambda^n E$$ where $\Lambda$ is the diagonal matrix of eigenvalues and $E = \sum_{i=0}^2 {e_ie_i^T}$. Then you can read the values off of independent entries from this result.