Given the following Hamiltonian:
$$\hat H = -\frac{\hbar^2}{2m}\frac{\partial^2}{\partial x_1^2} -\frac{\hbar^2}{2m}\frac{\partial^2}{\partial x_2^2} + \frac{1}{2}m\omega_1^2x_1^2 + \frac{1}{2}m\omega_2^2x_2^2 \, .$$
Where there is no coupled term $x_1x_2$, is the correct solution to assume a separation ansatz $\Psi_n = \psi_1(x_1)\psi_2(x_2)$, substitute in, divide through by $\psi_1(x_1)\psi_2(x_2)$ and get two equations of the form
$$\frac{\hbar^2}{2m}\frac{\partial^2\psi_i}{\partial x_i^2} + \frac{1}{2}m\omega_1^2x_i^2\psi_i = E_i\psi_i $$
where the solutions are the standard solutions to the quantum harmonic oscillator?
Following this the full solution $\Psi_n$ is the product of two independent solutions $\psi_1\psi_2$ with $E_{tot}=E_1+E_2$, so the total energy is a sum of energies but the discrete wavefunction is a product of discrete wavefunctions?