Computing (un)stable manifolds with validated error bounds: non-resonant and resonant spectra (Q310798)

This paper concerns the computation of the local stable manifold of a hyperbolic fixed point of a nonlinear ODE. (Unstable manifolds can be computed by simply reversing time.) The approach used results in a posteriori analysis of truncation errors, which can be combined with management of round off errors to give a mathematically rigorous description of the manifold. The authors' method works even when eigenvalues are resonant, which happens generically if a parameter of the system is varied. They demonstrate their method in three examples: (i) the stable manifold of the origin for the Lorenz system with standard parameters, (ii) the same but for repeated eigenvalues, (iii) the calculation of a connecting orbit derived from a pattern-forming PDE.