The problem of extracting dipole couplings from a system of N spins I = 1/2 and one spin 5 by NMR techniques is analyzed. The resolution attainable using a variety of single quantum methods is reviewed. The theory of heteronuclear multiple quantum (HMQ) NMR is developed, with particular emphasis being placed on the superior resolution available in HMQ spectra. Several novel pulse sequences are introduced, including a two-step method for the excitation of HMQ coherence. Experiments on partially oriented [1 - 13C] benzene demonstrate the excitation of the necessary HMQ coherence and illustrate the calculation of relative line intensities. Spectra of high order HMQ coherence under several different effective Hamiltonians achievable by multiple pulse sequences are discussed. A new effective Hamiltonian, scalar heteronuclear recoupled interactions by multiple pulse (SHRIMP), achieved by the simultaneous irradiation of both spin species with the same multiple pulse sequence, is introduced. Experiments are described which allow heteronuclear couplings to be correlated with an S-spin spreading parameter in spectra free of inhomogeneous broadening.