Combining Response Surface Methodology with Numerical Models for Optimization of Class-Based Queueing Systems