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