The MTGLE theory is formulated for spin systems using the exact low temperature limit solution of the one dimensional ferromagnetic Heisenberg model as a guide. The final equations are solved to reproduce the exact time correlation function and GLE response function for the spin 1/2 XY model in the high temperature limit [Sur et al., Phys. Rev. B 12, 3845 (1975)]. The advantage of the method is that it can be applied to complex systems with relative ease and in the case of spin systems provides an avenue for investigating effects of anisotropic dipolar spin interactions. One such application to exchange narrowing of the ESR line shape is discussed.