We introduce an algebraic model of coupled anharmonic oscillators capable of describing the complete spectrum and intensities of X-Y stretching vibrations of medium size and large polyatomic molecules. We apply this model to the calculation of C-H and C-D stretches of C6H6 and C6D6. The model appears to describe the data accurately.