Periodic orbits based on orthogonal collocation

We compute Ntst time slices of a periodic orbit using orthogonal collocation. This is implemented in the structure BifurcationKit.PeriodicOrbitOCollProblem.

Large scale

The current implementation is not yet optimized for large scale problems. This will be improved in the future.

The general method is explained in BifurcationKit.jl.