For the study of relaxation processes in complex spin system, a general master equation, which can be used to simulate a vast range of pulse experiments, has been formulated using the Liouville representation of quantum mechanics. The state of a nonequilibrium spin system in magnetic field is described by a density vector in Liouville space and the time evolution of the system is followed by the application of a linear master operator to the density vector in this Liouville space. In this master equation the nuclear spin relaxation due to intramolecular dipolar interaction or randomly fluctuating field interaction is explicitly implemented as a relaxation supermatrix for a strong coupled two-spin (1/2) system. The whole dynamic information inherent in the spin system is thus contained in the density vector and the master operator. The radiofrequency pulses are applied in the same space by corresponding unitary rotational supertransformations of the density vector. If the resulting FID is analytically Fourier transformed, it is possible to represent the final nonstationary spectrum using a frequency dependent spectral vector and intensity determining shape vector. The overall algorithm including relaxation interactions is then translated into an ANSIFORTRAN computer program, which can simulate a variety of two dimensional spectra. Furthermore a new strategy is tested by simulation of multiple quantum signals to differentiate the two relaxation interaction types.