From latitude, longitude (φ, λ) to UTM coordinates (E, N)

First we compute some intermediate values:

${\displaystyle t=\sinh \left(\tanh ^{-1}\left(\sin \varphi \right)-{\frac {2{\sqrt {n}}}{1+n}}\tanh ^{-1}\left({\frac {2{\sqrt {n}}}{1+n}}\sin \varphi \right)\right),}$

${\displaystyle \xi '=\tan ^{-1}\left({\frac {t}{\cos(\lambda -\lambda _{0})}}\right),\,\,\,\eta '=\tanh ^{-1}\left({\frac {\sin(\lambda -\lambda _{0})}{\sqrt {1+t^{2}}}}\right),}$

${\displaystyle \sigma =1+\sum _{j=1}^{3}2j\alpha _{j}\cos(2j\xi ')\cosh(2j\eta '),\,\,\,\tau =\sum _{j=1}^{3}2j\alpha _{j}\sin(2j\xi ')\sinh(2j\eta ').}$

The final formulae are:

${\displaystyle E=E_{0}+k_{0}A\left(\eta '+\sum _{j=1}^{3}\alpha _{j}\cos(2j\xi ')\sinh(2j\eta ')\right),}$

${\displaystyle N=N_{0}+k_{0}A\left(\xi '+\sum _{j=1}^{3}\alpha _{j}\sin(2j\xi ')\cosh(2j\eta ')\right),}$

${\displaystyle k={\frac {k_{0}A}{a}}{\sqrt {\left\{1+\left({\frac {1-n}{1+n}}\tan \varphi \right)^{2}\right\}{\frac {\sigma ^{2}+\tau ^{2}}{t^{2}+\cos ^{2}(\lambda -\lambda _{0})}}}},}$

${\displaystyle \gamma =\tan ^{-1}\left({\frac {\tau {\sqrt {1+t^{2}}}+\sigma t\tan(\lambda -\lambda _{0})}{\sigma {\sqrt {1+t^{2}}}-\tau t\tan(\lambda -\lambda _{0})}}\right).}$

where ${\displaystyle E}$ is Easting, ${\displaystyle N}$ is Northing, ${\displaystyle k}$ is the Scale Factor, and ${\displaystyle \gamma }$ is the Grid Convergence.

From UTM coordinates (E, N, Zone, Hemi) to latitude, longitude (φ, λ)

Note: Hemi = +1 for Northern, Hemi = −1 for Southern

First let's compute some intermediate values:

${\displaystyle \xi ={\frac {N-N_{0}}{k_{0}A}},\,\,\,\eta ={\frac {E-E_{0}}{k_{0}A}},}$

${\displaystyle \xi '=\xi -\sum _{j=1}^{3}\beta _{j}\sin \left(2j\xi \right)\cosh \left(2j\eta \right),\,\,\,\eta '=\eta -\sum _{j=1}^{3}\beta _{j}\cos \left(2j\xi \right)\sinh \left(2j\eta \right),}$

${\displaystyle \sigma '=1-\sum _{j=1}^{3}2j\beta _{j}\cos \left(2j\xi \right)\cosh \left(2j\eta \right),\,\,\,\tau '=\sum _{j=1}^{3}2j\beta _{j}\sin \left(2j\xi \right)\sinh \left(2j\eta \right),}$

${\displaystyle \chi =\sin ^{-1}\left({\frac {\sin \xi '}{\cosh \eta '}}\right).}$

The final formulae are:

${\displaystyle \varphi =\chi +\sum _{j=1}^{3}\delta _{j}\sin \left(2j\chi \right),}$

${\displaystyle \lambda _{0}=\mathrm {Z} \mathrm {o} \mathrm {n} \mathrm {e} \times 6^{\circ }-183^{\circ }\,}$

${\displaystyle \lambda =\lambda _{0}+\tan ^{-1}\left({\frac {\sinh \eta '}{\cos \xi '}}\right),}$

${\displaystyle k={\frac {k_{0}A}{a}}{\sqrt {\left\{1+\left({\frac {1-n}{1+n}}\tan \varphi \right)^{2}\right\}{\frac {\cos ^{2}\xi '+\sinh ^{2}\eta '}{\sigma '^{2}+\tau '^{2}}}}},}$

${\displaystyle \gamma =\mathrm {H} \mathrm {e} \mathrm {m} \mathrm {i} \times \tan ^{-1}\left({\frac {\tau '+\sigma '\tan \xi '\tanh \eta '}{\sigma '-\tau '\tan \xi '\tanh \eta '}}\right).}$