  We discuss the exact discretization of the classical harmonic oscillator
equation (including the inhomogeneous case and multidimensional
generalizations) with a special stress on the energy integral. We present and
suggest some numerical applications.
