A range of physical systems, particularly of chemical nature involving reactions, perform self-excited oscillations coupled by diffusion. The role of diffusion is not trivial so that initial differences in the phase of the oscillations between different points in space do not necessarily disappear as time goes; they may self-sustain. The dynamics of the phase depend on the values of the controlling parameters of the system. We consider a 6th-order nonlinear partial differential equation resulting in such dynamics. The equation is solved using central finite-difference discretization in space. The resulting system of ordinary differential equations is integrated in time using a Matlab solver. The numerical code is tested using forced versions of the equation, which admit exact analytical solutions. The comparison of the exact and numerical solutions demonstrates satisfactory agreement.