Action variable of a shifted harmonic oscillator (source code)

= Action variable of a shifted harmonic oscillator

For $F=(p^2+W^2q^2+aq)/2$, shifting $Q=q+a/(2W^2)$ gives oscillator energy $E=F+a^2/(8W^2)$ and action $I=E/|W|$.