Optimum airfoils in the sense of maximum lift coefficient are obtained for incompressible fluid flow at large Reynolds number. The maximum lift coefficient is achieved by requiring that the turbulent skin friction be zero in the pressure rise region on the airfoil upper surface. Under this constraint, the pressure distribution is optimized. The optimum pressure distribution is a function of Reynolds number and the trailing edge velocity. Geometries of those airfoils which will generate these optimum pressure distributions are obtained using a direct-iterative method which is developed in this study. This method can be used to design airfoils consisting of any number of elements. Numerical examples of one- and two-element airfoils are given. The maximum lift coefficients obtained range from 2 to 2.5.