Source code to solve numerically a classical harmonic oscillator