Time plot left, phase portrait right. Every driver is built by `externalsystem_from_symbolics` and integrated through `evaluate(sys.first_order_dynamics, r)` — no right-hand side is written twice.