Numerical method for solving stochastic differential equations with Poissonian white shot noise.

We propose a numerical integration scheme to solve stochastic differential equations driven by Poissonian white shot noise. Our formula, which is based on an integral equation, which is equivalent to the stochastic differential equation, utilizes a discrete time approximation with fixed integration time step. We show that our integration formula approaches the Euler formula if the Poissonian noise approaches the Gaussian white noise. The accuracy and efficiency of the proposed algorithm are examined by studying the dynamics of an overdamped particle driven by Poissonian white shot noise in a spatially periodic potential. We find that the accuracy of the proposed algorithm only weakly depends on the parameters characterizing the Poissonian white shot noise; this holds true even if the limit of Gaussian white noise is approached.