On the computation of unique normal forms and quadratic convergence

An algorithm is given which computes the transformation that transforms the initial system to final normalized system to a given order, while reducing the order to which intermediate results need to be computed. The algorithm will convergence quadratically when one is only interested in obtaining the classical normal form and somewhere in between linearly and quadratically in the general case.