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.