An improved adjoint operator method was proposed to compute the third order normal form of six dimensional nonlinear dynamical systems and the associated nonlinear transformation for the first time.