An iterative procedure for the calculation of the thick axisymmetric boundary layer and wake near the stern of a body of revolution is presented. Procedure consists of the potential flow calculation by a method of the integral equation of first kind and the calculation of boundary layer and wake by a differential me¬thod of the boundary layer theory. Additionally, higher order terms are included in the conventional momentum equations and continuity equation for the consider¬ation of the characteristics of axisymmetric flow different from the one of two dimentional flow and the thick boundary layer. These solutions are matched at the edge of boundary layer and wake. The results obtained by the present me¬thod are compared with the experimental data and it is found that the nominal wake distribution at the propeller plane of a axisymmetric body is in good agree¬ment with the experiment.