fun <- function(X){
D <- length(X)
f <- abs(sum(X^2) - D)^(1/4) + (0.5 * sum(X^2) + sum(X))/D + 0.5
return(f)
}
n <- 30
k <- 2
tolerance <- 0.01
lower <- c(-5,-5)
upper <- c(5,5)
res <- ExplicitExploration(fun, lower = lower,
upper = upper,n = n,
maxiter = 20,
k = k)
z <- outer(X = seq(-5, 5, 0.05), Y = seq(-5, 5, 0.05),
FUN = Vectorize(function(X, Y) { fun(c(X, Y)) }))
contour(seq(-5, 5, 0.05),seq(-5, 5, 0.05),z,
nlevels = 20, cex.axis = .8)
points(res$par_historical[,1],res$par_historical[,2],
col = "blue")
points(res$par[,1],res$par[,2], col = "red",
pch = 19)
Run the code above in your browser using DataLab