Custom multivariate density via DensityDist

You might need to tell DensityDist that the support dimension is 1 (a vector) and, if you have any parameters, what the core dimensions of each are.

But in your very simple case is something not working?