This paper develops a new computation method for obtaining a significant refinement of the normal forms for high dimensional nonlinear systems. Several explicit formulae are derived herein, which can be used to compute the coefficients of the simplest normal form and the associated nonlinear transformation. Through the illustrative example, the feasibility and merit of this novel method are also demonstrated and discussed.