Here, we assume an analytic expression of f is not available, and the outcome of f (X 1 ) can only be measured/recorded postexperimentation/simulation at X = X 1 . X may subject to known linear equality/inequality constraints. And the goal is to find the X * that minimizes f .