In fitting a model, there always exists a discrepancy between the fitted model and the true functional relationship. In measuring this discrepancy, Box and Drapper (1959) used the criterion dividing the discrepancy into two parts which are the bias error part and the variance error one in single response case. In this paper, an optimum design which makes these two types of errors as small as possible is found by extending the Box and Drapper criterion to multiple response situation. Especially, a design is found to meat rotatability conditions when we fit a quadratic model to each response fearing cubic bias. Using the central composite design, an application of general results to a specific case is shown to help understanding the material.