Exercises 12.1.7: Give the domain and range of the multi-variable function.

\[f(x,~y)= x^2+y^2 + 2\]

fun1 <- function(x,y) { x^2 +y^2 + 2 }

x <- runif(100000,-1000,1000)
y <- runif(100000,-1000,1000)
result <- fun1(x,y)

hist(result, breaks = 100)