\begin{Schunk}
\begin{Sinput}
 mu <- rep(0, 3)
 Sigma <- matrix(c(1, .3, -.1, .3, 1, .2, -.1, .2, 1), nrow = 3)
 Sigma
\end{Sinput}
\begin{Soutput}
     [,1] [,2] [,3]
[1,]  1.0  0.3 -0.1
[2,]  0.3  1.0  0.2
[3,] -0.1  0.2  1.0
\end{Soutput}
\begin{Sinput}
 mvrnorm(5, mu = mu, Sigma = Sigma)
\end{Sinput}
\begin{Soutput}
           [,1]       [,2]       [,3]
[1,] -0.1796238 -1.5100143  1.0723613
[2,] -0.3950370  1.1846379  1.3880722
[3,]  0.8681192 -0.1118966 -0.2419762
[4,]  0.3432726 -1.4629586 -1.5008948
[5,]  0.5923470 -0.3757720 -1.5577252
\end{Soutput}
\end{Schunk}
