Numerical method for solving stochastic differential equations with dichotomous noise.
We propose a numerical method for solving stochastic differential equations with dichotomous Markov noise. The numerical scheme is formulated such that (i) the stochastic formula used follows the Stratonovich-Taylor form over the entire range of noise correlation times, including the Gaussian white noise limit; and (ii) the method is readily applicable to dynamical systems driven by arbitrary types of noise, provided there exists a way to describe the random increment of the stochastic process expressed in the Stratonovich-Taylor form. We further propose a simplified Taylor scheme that significantly reduces the computation time, while still satisfying the moment properties up to the required order. The accuracies and efficiencies of the proposed algorithms are validated by applying the schemes to two prototypical model systems that possess analytical solutions.