These routines compute the initial conditions for a periodic orbit of a non-autonomous with periodic dependence on time:
That is, the goal is to compute a point such that, for a given ,
The method used consists on performing a Newton method to solve the equation
where is the stroboscopic map (time return map):
The program needs the gsl libraries, and can be compiled using