Thank you for an informative article. Do you have any suggestion for simulating multivariate normal data with an AR(p) correlation structure?