Function to generate an AR(1) variance-covariance matrix with parameter rho s.t. |rho| < 1.
ARcov(rho, n)
Useful links