I need to make a 2D plot for 3D data (for example a multivariate distribution and its density values) in a way that the third dimension (i.e. density values) is shown by a shaded area. any idea?

library(mvtnorm)

library(Compositional)

pairedData = rmvnorm(n=20,mean = c(0,0),sigma = diag(2))

z = mkde(pairedData, h=0.1)