Recent research suggests that any set of first order linear vector ODEs can be converted to a set of specific vector ODEs adhering to what we have called "Quantum Harmonical Form (QHF)". QHF has been developed using a virtual quantum multi harmonic oscillator system where mass and force constants are considered to be time variant and the Hamiltonian is defined as a conic structure over positions and momenta to conserve the Hermiticity. As described in previous works, the conversion to QHF requires the matrix coefficient of the first set of ODEs to be a normal matrix. In this paper, this limitation is circumvented using a space extension approach expanding the potential applicability of this method. Overall, conversion to QHF allows the investigation of a set of ODES using mathematical tools available to the investigation of the physical concepts underlying quantum harmonic oscillators. The utility of QHF in the context of dynamical systems and dynamical causal modeling in behavioral and cognitive neuroscience is briefly discussed.