This paper presents an explicit, computationally efficient, recursive formula for computing the normal forms, center manifolds and nonlinear transformations for general n-dimensional systems, associated with semisimple singularities. Based on the formula, we develop a Maple program, which is very convenient for an end-user who only needs to prepare an input file and then execute the program to "automatically" generate the result. Several examples are presented to demonstrate the computational efficiency of the algorithm. (C) 2013 Elsevier B. V. All rights reserved.