I'm trying to find the autocorrelation function for a discrete parameter Markov chain, $\{X(k)\}_{k=0}^{\infty}$ with a transition probability matrix given by $$ P = \begin{bmatrix} \frac{1}{3} & \frac{1}{3} & \frac{1}{3}\\ \frac{1}{3} & \frac{1}{3} & \frac{1}{3}\\ \frac{1}{3} & \frac{1}{3} & \frac{1}{3} \end{bmatrix} $$
and initial-state probabilities $\textbf{p}(0) = [\frac{1}{3}, \frac{1}{3}, \frac{1}{3}]$ and state space $E = \{0,1,2\}$
So the auto-correlation will be given, for $n\ge 1$, by $\mathbf{E}[X(k)X(k+n)]$. However I am not quite sure how to calculate this. For something like $\mathbf{E}[X(k)]$ one could use the definition of expectation and get $\sum_{k=0}^{3}X(k)\mathbf{P}(X(k)=k)$. However i am not quite sure how to deal with the product.