This function returns a matrix of draws from a multivariate normal distribution truncated on the interval [lb, ub].

rtmvnorm(n, mu, sigma, lb, ub)

See also