Duffing Oscillator:  ẍ + δẋ + αx + βx³ = γcos(ωt)

0.50
0.10
0.50
1.00
Left: phase portrait (x vs ẋ). Right: time series x(t). Adjust parameters to explore hardening/softening springs, limit cycles, subharmonic response, and the period-doubling route to chaos.