  We show that, by using the quantum orthogonal functions invariant, we are
able to solve a coupled of time dependent harmonic oscillators where all the
time dependent frequencies are arbitrary. We do so, by transforming the time
dependent Hamiltonian of the interaction by a set of unitary operators. In
passing, we show that $N$ time dependent and coupled oscillators have a
generalized orthogonal functions invariant from which we can write a
Ermakov-Lewis invariant.
