Multiple-shooting algorithm

What follows is the description of the multiple-shooting algorithm. Two functions are implemented, depending whther the orbit period is known a priori (i.e. for non-autonomous systems) or not (i.e. most of the autonomous systems).

Both methods are adaptations of the algorithms described by Lust [cite here].

Two methods to compute the Jacobian (or Monodromy matrix) are also described.

• Numerical scheme
• Numerical scheme (unknown period)
• Jacobian computation
  • Analytical computation
  • Numerical computation
• Multiple-shooting solver
• Floquet multipliers