Time plot left, phase portrait right. Orbits integrated with RK4 through `evaluate(sys.first_order_dynamics, r)`.