d²φ/dt² + (1/Q)dφ/dt + sin(φ) = I_dc + I_ac·cos(ωt)
Poincaré Section (stroboscopic)
Phase velocity dφ/dt vs time