Numerical integration of the non-stationary Schrödinger equation with Duffing potential depending on two coordinates has been carried out. Oscillation types and the influence of coupling between two oscillators on frequency spectra are analyzed in detail.