Numerical Taylor expansions for invariant manifolds

We consider numerical computation of Taylor expansions of invariant manifolds around equilibria of maps and flows. These expansions are obtained by writing the corresponding functional equation in a number of points, setting up a nonlinear system of equations and solving this system using a simplified Newton’s method. This approach will avoid symbolic or explicit numerical differentiation. The linear algebra issues of solving the resulting Sylvester equations are studied in detail.