This paper describes EXPRESS (EXchange Program for RElaxing Spin Systems), a computer program that simulates the effects of Markovian jump dynamics for a wide variety of solid state nuclear magnetic resonance experiments. A graphical interface is described that facilitates the definition of rotational jumps around non-commuting axes that may occur at arbitrary, different rates. Solid state deuteron NMR is widely used to investigate such processes, and EXPRESS allows simulations of deuteron quadrupole echo and magic angle spinning line (MAS) shapes, as well as partially relaxed line shapes for measurements of anisotropic relaxation of Zeeman and quadrupolar order. Facilities are included for chemical shift tensors (at user-defined orientations relative to the quadrupole coupling tensors), so that EXPRESS is potentially useful for studies of paramagnetic systems where these interactions are of comparable magnitude. Many of the same techniques used for deuterons can be extended to half-integer quadrupolar nuclei. The same interface that specifies rotational "sites" for deuteron NMR studies is usable in EXPRESS to simulate static line shapes, MAS line shapes, and multiple pulse Carr-Purcell-Meiboom-Gill (CPMG) line shapes for the central transition of high spin quadrupoles with second order quadrupole coupling and chemical shift anisotropy. Representative simulations are displayed that show effects of slow libration on deuteron quadrupole echo line shapes and relaxation time anisotropies. EXPRESS is also used to investigate fundamental differences in the mechanism of echo formation in deuteron MAS and quadrupole CPMG experiments, and to illustrate significant differences between these techniques in the context of high spin quadrupolar nuclei. The program is modular and platform-independent, with facilities for users to add routines for experiments not yet envisioned.