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