The transverse cylindrical projection has been used in Korea since 1910s when the nationwide geodetic network was firstly established. However, the projection has a number of different types of functions according to a way of its mathematical derivation as well as a section of its coefficients and terms, for instance Gauss- Schreiber(GS) and Gauss-Kruger(GK) types. Although the transverse cylindrical projection itself is assigned to a system, projected coordinates would be diverse with respect to the function used in the actual calculation. In order to investigate impact of functions used in the computation, five different equations (i.e., 2 GS and 3 GK) were implemented in this study by using MATLAB. A series of numerical simulation tests has been carried out to compare and characterize them in terms of projection accuracy, difference of projected coordinates and distortion. Furthermore, a comparison between GS and GK function was made under the Korean gridding system, consisting of four zones. Results from the numerical computations were qualitatively analyzed and summarized in this paper.