Started calculating Floquet multipliers. Jacobian calculation becomes unstable for long orbits (after 2-cycles), I think I will survive by dividing it into subintervals. Will do that tomorrow.