An integral equation constrained optimization approach to finding minimum-drag shapes for solid bodies of revolution in Stokes flows subject to constraints on body's volume and shape has been developed. An axially symmetric Stokes flow problem has been formulated in the form of an integral equation (state equation), and finding an optimal shape has been reduced to an integral equation constrained optimization problem. The total variation of the Lagrangian functional of the problem has been reduced to the variation only with respect to the shape by the adjoint equation-based method, and the steepest descent direction for the coefficients in the function series representing the shape has been found analytically. It has been shown that solutions to the state and adjoint integral equations are related by a simple algebraic formula, which eliminates the need for solving the adjoint equation. The approach has been illustrated in finding the minimum-drag shape for a solid unit volume particle and in finding the ...