From the viewpoint of obtaining low viscous drag at large Reynolds numbers for axisymmetric bodies over which the boundary layer is fully turbulent, an iteration method is developed for solving the inverse problem of getting body shape from the velocity distribution by means of boundary integral equation for Stokes stream function. A calculation formula, an integral of a power of the potential-flow surface velocity, is used to get the viscous drag. Using an approximate steepest-descent-method in the optimum value problem the minimum drag shaping is conducted. The obtained value of the minimum drag and the body shape is nearly equal to ones known already. The presently obtained know-hows and means are expected to obtain the body shapes over which the boundary layer is partially laminar.
