Application of the resonant normal form to high-order nonlinear ODEs using Mathematica