\begin{Schunk}
\begin{Sinput}
    library(mvtnorm)
    set.seed(12345)
    X0 <- mvtnorm::rmvnorm(n=10, mean = c(0,0,0), sigma = diag(3))
    ## create a smaller data set, starting at same position
    set.seed(12345)
    X1 <- mvtnorm::rmvnorm(n=5,  mean = c(0,0,0), sigma = diag(3))
    X0
\end{Sinput}
\begin{Soutput}
            [,1]       [,2]       [,3]
 [1,]  0.5855288  0.7094660 -0.1093033
 [2,] -0.4534972  0.6058875 -1.8179560
 [3,]  0.6300986 -0.2761841 -0.2841597
 [4,] -0.9193220 -0.1162478  1.8173120
 [5,]  0.3706279  0.5202165 -0.7505320
 [6,]  0.8168998 -0.8863575 -0.3315776
 [7,]  1.1207127  0.2987237  0.7796219
 [8,]  1.4557851 -0.6443284 -1.5531374
 [9,] -1.5977095  1.8050975 -0.4816474
[10,]  0.6203798  0.6121235 -0.1623110
\end{Soutput}
\begin{Sinput}
    X1
\end{Sinput}
\begin{Soutput}
           [,1]       [,2]       [,3]
[1,]  0.5855288  0.7094660 -0.1093033
[2,] -0.4534972  0.6058875 -1.8179560
[3,]  0.6300986 -0.2761841 -0.2841597
[4,] -0.9193220 -0.1162478  1.8173120
[5,]  0.3706279  0.5202165 -0.7505320
\end{Soutput}
\end{Schunk}
