Boundary layer solutions around body of revolution are analyzed throughout this study. Sphere, ellipsoids that have different fineness ratios and bodies which are obtained by the Karman-Trefftz airfoils that are revolved around their axis are used as axially symmetrical bodies. The velocity distributions that are needed in order to solve boundary layer equations are obtained by the solutions of potential flow. Therefore, an analytic solution of potential flow around ellipsoid and axial singularity distribution method are used. Thwaites-Walz method, which gives approximate solutions for all geometries, is applied to solve boundary layer equations. In order for the solutions to be tested, an analytic solution that gives an exact solution for sphere which is a special condition of ellipsoid is utilized. In order to obtain solutions in a quick way, a software program is developed in Visual Basic.