The homonuclear dipolar coupling is the internal spin interaction that contributes the most to the line shapes in magic-angle-spinning (MAS) H-1 NMR spectra of solids, and linewidths typically extend over several hundred Hertz, limiting the H-1 resolution. Understanding and reducing this contribution could provide rich structural information for organic solids. Here, we use average Hamiltonian theory to study two- and three-spin systems in the fast MAS regime. Specifically, we develop analytical expressions to third order in the case of two and three inequivalent spins (I = 1/2). The results show that the full third-order expression of the Hamiltonian, without secular approximations or truncation to second order, is the description that agrees the best, by far, with full numerical calculations. We determine the effect on the NMR spectrum of the different Hamiltonian terms, which are shown to produce both residual shifts and splittings in the three-spin systems. Both the shifts and splittings have a fairly complex dependence on the spinning rate with the eigenstates having a polynomial.r dependence. The effect on powder line shapes is also shown, and we find that the anisotropic residual shift does not have zero average so that the powder line shape is broadened and shifted from the isotropic position. This suggests that in H-1 MAS spectra, even at the fastest MAS rates attainable today, the positions observed are not exactly the isotropic shifts. (c) 2021 Author(s).