\hat{H}=-\frac{\hbar^2}{2m}(\frac{\partial^2}{\partial x^2}+\frac{\partial^2}{\partial y^2})+\frac{1}{2}m\omega^2x^2