second-order polynomial models are
developed because first-order models often give
lack-of-fit [24]. The selected approximation
model is then used to find optimal setting of
input (design) variables that maximize (or
minimize) the mean value of the response.
However, when the variance is not constant,
classical RSM can be misleading [25].