\begin{Schunk}
\begin{Sinput}
 mu <- c(0, 0, 0, 0)
 N <- 1000
 X <- mvrnorm(N, m = mu, S = Sigma) 
\end{Sinput}
\end{Schunk}
