dmsn function

Multivariate skew-normal distribution