Continuation of time-periodic solutions
This should involve setting up an augmented system, that treats time as a spatial coordinate and has periodic boundary conditions. It is propably okay to use some finite-difference (forward/central?) scheme for time-discretization.
The time axis is normalized to 1 but scaled by some parameter tau, that is a free parameter. The initial guess comes from the imaginary part of the eigenvalues corresponding to the Hopf bifurcation.
I am unsure, whether this can be realized with a heavy usage of the Equation framework or whether we need a more intricate code.