Using the compact matrix forms of third order susceptibility, we have derived the expressions of effective third order nonlinear coefficients for third harmonic generation (THG) processes in uniaxial crystals. These expressions are significant in finding the optimum incident directions for THG in uniaxial crystals.