Skip to content

Numerical approximations of strong (un)stable manifolds

Research output: Working paperWorking paper and Preprints

  • Hinke M Osinga
  • Reza Rokni
  • S Townley
Original languageEnglish
StatePublished - 2002


It is well-known how to compute global one-dimensional stable or unstable manifolds of a hyperbolic equilibrium of a smooth vector field. Such manifolds consist only of two trajectories and arbitrarily large pieces can be generated using an initial point close to the equilibrium on the linear approximation of the manifold. The attraction properties (in forward or backward time) of the local manifolds ensure that the computation error, which depends on the arclength of the computed piece, remains bounded. This paper discusses how these error bounds change as the equilibrium loses its hyperbolicity, or when the one-dimensional, say, stable manifold is in fact a strong stable manifold that is contained in a higher-dimensional stable manifold. For these cases, the local manifolds are not locally attracting (either in forward or in backward time) and the standard error bound does not work. The theoretical analysis is illustrated with numerical computations.

Additional information

Additional information: Preprint of a paper later published by Taylor and Francis (2004), Dynamical Systems, 19(3), pp.195-215, ISSN 1468-9367

View research connections

Related faculties, schools or groups