Direct quadrature method of moments solution of the Fokker–Planck equation

The direct quadrature method of moments is presented as an efficient and accurate means of numerically computing solutions of the Fokker–Planck equation corresponding to stochastic nonlinear dynamical systems. The theoretical details of the solution procedure are first presented. The method is then used to solve Fokker–Planck equations for both 1D and 2D (noisy van der Pol oscillator) processes which possess nonlinear stochastic differential equations. Higher-order moments of the stationary solutions are computed and prove to be very accurate when compared to analytic (1D process) and Monte Carlo (2D process) solutions.