Consider the model $$x_i=\rho_1x_{i-1}+\rho_2x_{i-2}+\dots+\rho_px_{i-p}+\omega_i,\:\mathbf{\omega}\sim N(\vec0,\mathbf{I}\sigma^2)$$
In the case of order 1 autocorrelation (i.e. where $\rho_2$ and up are 0), I am able to make a covariance matrix of $\mathbf{x}$ quite easily (see here). But I need to be able to compute the variances and covariances for any order of autocorrelation, and I have not been able to find any literature that covers the math behind it. Does anyone know where I can find out how to compute this matrix?