Abstract
To describe turbulent wake behind a submarine, it is very important to know turbulent kinetic energy distributions in the wake. To get the distribution is to solve the turbulent kinetic energy equation, and to solve the equation, it is needed both information of ${\lambda}$ and ${\sigma}$ which define physical characteristics of the wake. This paper gives analytical solution of the equation, which is driven from $8^{th}$ order polynomial fitting, as a function of given ${\lambda}$, even though there is no information of ${\sigma}$. In comparison between numerical solution(i.e. exact solution) and analytical solution, the relative errors between them are less than to 5% in the range of 0 < ${\xi}$ < 0.95 in most given ${\lambda}$.