Eleven different matrix forms of the third-order nonlinear susceptibility for all crystal classes are obtained by contracting subcript method for the first time.Meanwhile,the simplified forms under Kleinman symmetry condition are also presented.Calculation methods of the third-order effective nonlinear susceptibility in crystals when four-wave interactions occur are studied theoretically.