which included eight tests for factorial design, six tests for axial points and 6 tests for replication of the central point (Table 3). A second order polynomial response surface was fitted using the methods and tools described by Lenth (Lenth, 2009). The response surface model has the form: