A method which generates a multiphonon basis through the construction and the iterative solution of a set of equations of motion is adopted to investigate the negative parity spectrum and the structure of the giant dipole resonance in 16O.