Let $X_t$ be a causal $AR(p)$ process. Compute a linear forecast $X_{t+2}$ based on $X_1, X_2, ..., X_t$ for $t \geq p+1$.
If $AR(p)$ is causal it means that it can be rewritten as a linear process: $X_t=\sum\limits_{i=0}^{\infty}\psi_iZ_{t-i}$ where $\sum\limits_{i=0}^{\infty}|\psi_i|<\infty$ and $Z_t$ is a white noise.
My prediction should be a linear combination of $X_1,...,X_t: \hat X_{t+2}=m+\sum\limits_{i=1}^t a_iX_{i}$ so it is a projection onto space spanned by $(X_1, ..., X_n)$. But how to proceed?