Skip to content

Latest commit

 

History

History
88 lines (86 loc) · 2.96 KB

README.md

File metadata and controls

88 lines (86 loc) · 2.96 KB
Raw

Computation of periodic orbits for non-autonomous systems

These routines compute the initial conditions for a periodic orbit of a non-autonomous with periodic dependence on time:

$\displaystyle \dot{x}=f(x,t),\quad f(x,t+T)=f(x,t).$    

That is, the goal is to compute a point $ x^*$ such that, for a given $ t_0$,

$\displaystyle \varphi(t_0+qT;x)=x.$    

The method used consists on performing a Newton method to solve the equation

$\displaystyle S_{t_0}(x)-x=0,$    

where $ S_{t_0}$ is the stroboscopic map (time $ qT$ return map):

$\displaystyle S_{t_0}(x)=\varphi(t_0+qT;x)$    

The program needs the gsl libraries, and can be compiled using
g++ -o main main.c rk78.c lu.c memory.c -lgsl -lgslcblas