I have a function "f" as follows:
f <- function (X){ DD= 91.8372-10.3096*(as.numeric((X[1]<70.5)&(X[2]>69.5)&(X[2]<79.5)&(X[3]