Baseado nesta dica, resolvi experimentar uma Cobb-Douglas simples.

u <- function(x, y) (x^0.5)*(y^0.5)

x <- seq(.1, 40, by=.1)
y <- seq(.1, 40, by=.1)
a <- c(1,5,10,15, 20, 30)

persp(x, y, outer(x, y, u), ticktype=”detailed”)
contour(x, y, outer(x, y, u), levels=a)

Resultado? Olha que bonitinhas.