## Not run:
# f <- function(xy) {
# x <- xy[1]; y <- xy[2]
# 3*(1-x)^2 * exp(-(x^2) - (y+1)^2) -
# 10*(x/5 - x^3 - y^5) * exp(-x^2 - y^2) -
# 1/3 * exp(-(x+1)^2 - y^2)
# }
# ezcontour(f, col = "navy")
# ezcontour(f, filled = TRUE)
# ezmesh(f)
# ezmesh(f, col="lightblue", theta = -15, phi = 30)
# ## End(Not run)
Run the code above in your browser using DataLab