From StatWiki
pdf(rpdf)
library(mvtnorm)
dmvnorm(c(0,0), c(0,0), diag(1,2))
x=seq(-4,4,length.out=100)
z=outer(x, x, function(x,y) dmvnorm(cbind(x,y), c(0,0), diag(1,2)) )
persp(x,x,z,col="lightblue", xlab="", zlab="", ylab="", box=F, border=NA, shade=0.75, theta=-60)
#