We develop an efficient and practical three-dimensional design procedure based on discrete sensitivity analysis and capable of handling large numbers of design variables. The Euler equations are solved using a fully implicit unfactored algorithm. This new procedure is applied to the design of three-dimensional transport wings in transonic flow. We use a wing geometry model that is totally based on two- and three-dimensional Bezier-Bernstein parametrizations. Two wing design cases are presented; one case illustrates the procedure’s suitability to preliminary design, and the other demonstrates its ability to produce realistic realistic optimal shapes, even when starting from very elementary initial geometries.