In this work we study 2-and 3-body oscillators with quadratic and sextic pairwise potentials which depend on relative distances, |r i − r j |, between particles. The two-body harmonic oscillator is two-parametric and can be reduced to a one-dimensional radial Jacobi oscillator, while in the 3-body case such a reduction is not possible in general. Our study is restricted to solutions in the space of relative motion which are functions of mutual (relative) distances only (S-states). We pay special attention to the cases where the masses of the particles and spring constants are unequal as well as to the atomic, where one mass is infinite, and molecular, where two masses are infinite, limits. In general, three-body harmonic oscillator is 7-parametric depending on 3 masses and 3 spring constants, and frequency. In particular, the first and second order integrals of the 3-body oscillator for unequal masses are searched: it is shown that for certain relations involving masses and spring constants the system becomes maximally (minimally) superintegrable in the case of two (one) relations.(J) 1 is always proportional to the vector of relative distance between particles 1 and 2. Evidently, the variables in T r are separated and the kinetic energy of relative motion is the sum of kinetic energies in the Jacobi coordinate directions.By adding to T r the harmonic oscillator potentials in each Jacobi coordinate direction we