I want to simulate three variables e1, e2 and e3 that are not independent of each other. The population mean of each variable must be 0 and each follows a normal distribution. Since they are not independent, the entries outside the main diagonal in the covariance matrix are not zeros. Is there a way to achieve this in R?