curretal88exp <- function(xx) { ########################################################################## # # CURRIN ET AL. (1988) EXPONENTIAL FUNCTION # # Authors: Sonja Surjanovic, Simon Fraser University # Derek Bingham, Simon Fraser University # Questions/Comments: Please email Derek Bingham at dbingham@stat.sfu.ca. # # Copyright 2013. Derek Bingham, Simon Fraser University. # # THERE IS NO WARRANTY, EXPRESS OR IMPLIED. WE DO NOT ASSUME ANY LIABILITY # FOR THE USE OF THIS SOFTWARE. If software is modified to produce # derivative works, such modified software should be clearly marked. # Additionally, this program is free software; you can redistribute it # and/or modify it under the terms of the GNU General Public License as # published by the Free Software Foundation; version 2.0 of the License. # Accordingly, this program is distributed in the hope that it will be # useful, but WITHOUT ANY WARRANTY; without even the implied warranty # of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU # General Public License for more details. # # For function details and reference information, see: # http://www.sfu.ca/~ssurjano/ # ########################################################################## # # INPUT: # xx = c(x1, x2) # ######################################################################### x1 <- xx[1] x2 <- xx[2] fact1 <- 1 - exp(-1/(2*x2)) fact2 <- 2300*x1^3 + 1900*x1^2 + 2092*x1 + 60 fact3 <- 100*x1^3 + 500*x1^2 + 4*x1 + 20 y <- fact1 * fact2/fact3 return(y) }