TW 251

K. Engelborghs, K. Lust, D. Roose
A Newton-Picard method for accurate computation of period doubling bifurcation points of large-scale systems of ODEs

Abstract

Periodic solutions of certain large scale systems of ODEs can be computed efficiently using a hybrid Newton-Picard scheme, especially in a continuation context. In this paper we describe and analyse a numerical method to accurately compute period doubling bifurcation points using the Newton-Picard approach. The method avoids the computation of the full monodromy matrix, which is present in the determining system for period doubling bifurcation points. Test results are presented that demonstrate the numerical properties.

report.pdf / mailto: K. Engelborghs