\begin{Schunk}
\begin{Sinput}
 library(rockchalk)
 myR <- lazyCor(X = 0.3, d = 5)
 myR
\end{Sinput}
\begin{Soutput}
     [,1] [,2] [,3] [,4] [,5]
[1,]  1.0  0.3  0.3  0.3  0.3
[2,]  0.3  1.0  0.3  0.3  0.3
[3,]  0.3  0.3  1.0  0.3  0.3
[4,]  0.3  0.3  0.3  1.0  0.3
[5,]  0.3  0.3  0.3  0.3  1.0
\end{Soutput}
\begin{Sinput}
 mySD <- c(0.5, 0.5, 0.5, 1.5, 1.5)
 myCov <- lazyCov(Rho = myR, Sd = mySD)
 myCov
\end{Sinput}
\begin{Soutput}
      [,1]  [,2]  [,3]  [,4]  [,5]
[1,] 0.250 0.075 0.075 0.225 0.225
[2,] 0.075 0.250 0.075 0.225 0.225
[3,] 0.075 0.075 0.250 0.225 0.225
[4,] 0.225 0.225 0.225 2.250 0.675
[5,] 0.225 0.225 0.225 0.675 2.250
\end{Soutput}
\begin{Sinput}
 myMu <- c(1.1, 2.0, 1.1, -0.2, 0)
 ## Draw one to see what that does
 set.seed(123123)
 mvrnorm(1, mu = myMu, myCov)
\end{Sinput}
\begin{Soutput}
[1]  1.2809807  2.1131828  0.9241272 -0.2839534 -0.4944292
\end{Soutput}
\end{Schunk}
