x <- seq(-3, 3, length = 50)
y <- x
otu <- function(x, y) 1/8 * (6 * exp(-((2/3 * abs(x) - 1)^2 + (2/3 *
y)^2) - 1/3 * (2/3 * y + 1/2)^3) + 2/3 * exp(-2.818^11 * ((abs(2/3 * x) -
1)^2 + (2/3 * y)^2)^2) + 2/3 * y - (2/3 * x)^4)
z <- outer(x, y, otu)
persp(x, y, z, theta = 30, phi = 30, expand = 0.5, col = "green")
x <- seq(-10, 10, length = 50)
y <- x
z <- outer(x, y, function(x, y) sin(sqrt(x^2 + y^2))/sqrt(x^2 + y^2))
persp(x, y, z, theta = 30, phi = 30, expand = 0.5, col = "green")
x <- seq(-2, 2, length = 50)
y <- x
z <- outer(x, y, function(x, y) (exp(-(x^2 + y^2)/2) * (x^2 + y^2) *
(2 * pi)))
persp(x, y, z, theta = 30, phi = 30, expand = 0.5, col = "green")
