First we construct the linear invariants of the system using equations of motion. These invariants can be renormalized to be transformed into ladder operators, which are also invariants. Next we construct the quadratic invariants in the same way as the linear ones. The quadratic invariant turns out to be proportional to the number operator, which makes possible to construct the harmonic oscillator algebra.