Forced Duffing — dx/dt=y, dy/dt=f cos θ−β x−x³−d y, dθ/dt=ω. Classic f=27 ω=2.2 d=0.2 β=1. Distinct from duffing.html / ueda.html. RK4. Gold trail, cap 650. RM = settled silhouette.