The presence of numerical noise inhibits gradient-based optimization and therefore limits the practicality of performing aircraft multidisciplinary design optimization (MDO). To address this issue, a procedure has been developed to create noise free algebraic models of subsonic and supersonic aerodynamic performance for use in the MDO of high-speed civil transport (HSCT) configurations. This procedure employs methods from statistical design of experiments theory to select a set of HSCT wing designs (fuselage/tail/engine geometry fixed) for which numerous detailed aerodynamic analyses are performed. Polynomial approximations (i.e., response surface models) are created from the aerodynamic data to provide analytical models relating aerodynamic quantities (e.g., wave drag and drag-due-to-lift) to the variables which define the HSCT wing configuration. A multidisciplinary design optimization of the HSCT is then performed using the response surface models in lieu of the traditional, local gradient based design methods. The use of response surface models makes possible the efficient and robust application of MDO to the design of an aircraft system. Results obtained from five variable and ten variable wing design problems presented here demonstrate the effectiveness of this response surface modeling method.