\begin{Schunk}
\begin{Sinput}
 Sigma.R <- chol(Sigma)
 Sigma.R
\end{Sinput}
\begin{Soutput}
     [,1]     [,2]      [,3]      [,4]
[1,]    1 1.000000  4.200000  2.800000
[2,]    0 1.732051 -3.233162 -1.154701
[3,]    0 0.000000  4.572381 -2.163716
[4,]    0 0.000000  0.000000  1.464582
\end{Soutput}
\end{Schunk}
