Custom Distribution - Multivariate SkewNormal

Yeah could be tricky, special because it looks like it may not have a closed form solution so you need to do some approximation