A spin system with unperturbed Hamiltonian ${\mathcal{H}}_{0}=\frac{1}{2}\ensuremath{\Sigma}\ensuremath{\Sigma}{B}_{\mathrm{jk}}{{\ensuremath{\sigma}}_{x}}^{j}{{\ensuremath{\sigma}}_{x}}^{k}\ensuremath{-}\ensuremath{\gamma}H\ensuremath{\Sigma}{{\ensuremath{\sigma}}_{x}}^{j}$ relaxing via the spin-lattice coupling $G=\frac{1}{2}\ensuremath{\Sigma}\ensuremath{\Sigma}{C}_{\mathrm{jk}}{{\ensuremath{\sigma}}_{x}}^{j}({{\ensuremath{\sigma}}_{y}}^{k}+{{\ensuremath{\sigma}}_{z}}^{k})$ is studied by means of the general density-matrix theory of magnetic relaxation. By making some assumptions about the magnitude and time constants of the lattice correlation functions
〈Cij(t)Ckl(0)〉, a master equation is obtained. It agrees at high temperatures with a master equation previously suggested by Glauber for the one-dimensional nearest-neighbor case. At high temperatures the magnetic moment relaxes with a single relaxation time, and the spin pair-correlation functions satisfy a closed set of equations. At low temperatures, however, the equations for the magnetization and the correlation functions are coupled to higher-order moments.