This function generates a random draw from the Wishart distribution.
rwish(v, S)
v
S
A (p x p) matrix.
dwish
A <- rwish(3, matrix(c(1,.3,.3,1),2,2))