Numerical treatment in Prandtl lifting-line theory of the nonlinearity associated with a 2-D lift curve, when the local incidence is larger than the incidence of maximum lift, is proposed. It is shown that the use of an artificial viscosity term makes the solution unique and allows the iterative method to converge to a physically meaningful solution, that is in agreement with the exact solution for the test case. The design and analysis of winglets is presented. The winglets considered are small fences placed upward at the tip of the wing to improve the wing efficiency by decreasing the induced drag. The effect of yaw on a wing equipped with such optimal winglets indicates that they provide weathercock stability.