mThe energy eigenvalues of coupled oscillators in two dimensions with quartic and sextic couplings have been calculated to a high accuracy. For this purpose, unbounded domain of the wave function has been truncated and various combination of trigonometric functions are employed as the basis sets in a Rayleigh-Ritz variational method. The method is applicable to the multiwell oscillators as well. 0 1996 John Wiley & Sons, Inc.