\begin{Schunk}
\begin{Sinput}
 library(mvtnorm)
 meanV <- c(0,0,0,0)
 varV <- 3*diag(4)
 X <- rmvnorm(10, mean = meanV, sigma = varV)
 W <- t(X) %*% X
 W
\end{Sinput}
\end{Schunk}
