Abstract
This paper describes the construction of a numerical code which is based solely on the product operator formalism and is aimed at simulating the full dynamical behaviour of nuclear spin systems under the influence of time-dependent radiofrequency fields, relaxation and exchange processes. To demonstrate the capabilities of the program, application of a sequence of selective shaped pulses to an AMX spin system is presented and analysed.