Fixed a bug with rmvnorm and rlmvnorm where an incorrect number
of variates were simulted when all boundaries were -Inf to
`Inf
Import from MASS; add rmassnorm and
rlmassnorm to simulate multivariate normal from the MASS
function
Added several tests