\begin{Schunk}
\begin{Sinput}
 library(MASS)
 mu <- c(0, 0, 0)
 Sigma <- diag(c(1, 1, 1))
 Sigma
\end{Sinput}
\begin{Soutput}
     [,1] [,2] [,3]
[1,]    1    0    0
[2,]    0    1    0
[3,]    0    0    1
\end{Soutput}
\begin{Sinput}
 mvrnorm(5, mu, Sigma)
\end{Sinput}
\begin{Soutput}
           [,1]       [,2]       [,3]
[1,] -0.1162478 -1.8179560  0.5855288
[2,]  1.8173120  0.6300986  0.7094660
[3,]  0.3706279 -0.2761841 -0.1093033
[4,]  0.5202165 -0.2841597 -0.4534972
[5,] -0.7505320 -0.9193220  0.6058875
\end{Soutput}
\end{Schunk}
