[R] Plot bivariate density with densities margins

[Eduardo García Portugués]

Dear R users,

I would like to plot the bivariate density surface with its marginal
densities on the sides of the 3D box, just like in the picture I attach. I
tried to found information about how to do it but did not find anything.

Does anyone know how to do it?

Thanks in advance,

Eduardo.