Numerical path integration of a non-homogeneous Markov process

The numerical path integration method, based on Gauss-Legendre integration scheme, is applied to a Duffing oscillator subject to both sinusoidal and white noise excitations. The response of the system is a Markov process with one of the drift coefficients being periodic. It is a non-homogeneous Markov process that does not have a stationary probability distribution. When applying the numerical procedure, the values of transition probability density at the Gaussian-Legendre quadrature points need only be calculated for time steps of the first period of the sinusoidal excitation, and they can be saved for use in all subsequent periods. The numerical procedure is capable of capturing the evolution of the probability density from an initial distribution to one that is changing and rotating periodically in the phase space.