changes made according to Dirk's comments

However, I don't see the relationship to multivariate Normal draws.
Sure, both could be put into an actual package, but I don't see
relationship between the codes.