Solution (source code)

= Solution

Spatial <derivatives> vanish for $u=y(t)$, giving $\ddot y+\gamma\dot y+\omega^2y=f(t)$.