Optimization of Highly Non-Linear Simulation Functions Using Improved Adaptive Response Surface Methodology