Numerical continuation of connecting orbits of maps in MATLAB

R. Khoshsiar Ghaziani, W. Govaerts, Iouri Aleksandrovitsj Kouznetsov, Hil Gaétan Ellart Meijer

    Research output: Contribution to journalArticleAcademicpeer-review

    18 Citations (Scopus)
    8 Downloads (Pure)


    We present new or improved methods to continue heteroclinic and homoclinic orbits to fixed points in iterated maps and to compute their fold bifurcation curves, corresponding to the tangency of the invariant manifolds. The proposed methods are applicable to general n-dimensional maps and are implemented in matlab. They are based on the continuation of invariant subspaces (CIS) algorithm, which is presented in a novel way. The systems of defining equations include the Riccati equations appearing in CIS for bases of the generalized stable and unstable eigenspaces. We use the bordering techniques to continue the folds, and provide full algorithmic details on how to treat the Jacobian matrix of the defining system as a sparse matrix in matlab. For a special - but important in applications - case n = 2 we describe the first matlab implementation of known algorithms to grow one-dimensional stable and unstable manifolds of the fixed points of noninvertible maps. The methods are applied to study heteroclinic and homoclinic connections in the generalized Hénon map.
    Original languageUndefined
    Article number10.1080/10236190802357677
    Pages (from-to)849-875
    Number of pages27
    JournalJournal of difference equations and applications
    Issue number8-9
    Publication statusPublished - 2009


    • IR-68366
    • METIS-263696
    • EWI-13858

    Cite this