A software package has been designed to simulate nuclear magnetic resonance spectra and images. Combining the product operator matrix with the non-linear Bloch equations, the software can efficiently simulate classical and quantum effects including scalar coupling, dipolar coupling, translational diffusion, chemical shift, radiation damping, transverse relaxation, and longitudinal relaxation.