require(MASS) #simulate from a multivariate normal distribution
for (i in 1:n)
{
class[i] = sample(1:k, 1, prob=pi)
X[i,] = mvrnorm(1, meanX, covX)
require(MASS) #simulate from a multivariate normal distribution
for (i in 1:n)
{
class[i] = sample(1:k, 1, prob=pi)
X[i,] = mvrnorm(1, meanX, covX)