A new algorithm for computing one-dimensional stable and unstable manifolds of maps

A new algorithm is presented to compute one-dimensional stable and unstable manifolds of fixed points for both two-dimensional and higher dimensional diffeomorphism maps. When computing the stable manifold, the algorithm does not require the explicit expression of the inverse map. The global manifold is grown from a local manifold and one point is added at each step. The new point is located with a "prediction and correction" scheme, which avoids searching the computed part of the manifold with a bisection method and accelerates the searching process. By using the fact that the Jacobian transports derivatives along the orbit of the manifold, the tangent component of the manifold is determined and a new accuracy criterion is proposed to check whether the new point that defines the manifold is acceptable. The performance of the algorithm is demonstrated with several numerical examples.