ÿ + cẏ + ky = forcing(r) with c = 0.2, k = 1.0, integrated together with ṙ = E(r) as a single state. The forcing reads the physical external state through `to_physical_external`.