clear

scalar sd = sqrt(6)
matrix sdm = (sd,sd)
matrix m = (4,4)
drawnorm z1 z2, n(15) means(m) sds(sdm)
summarize 

clear
matrix sd = sqrt(6)
drawnorm z3, n(200) means(4) sds(sd) 
summarize z3
histogram z3, normal