Multivariate density estimation