A numerical method, based on the path-integral formalism, is presented to solve non-linear Fokker-Planck equations with natural boundary conditions. For one-dimensional stochastic processes, several specific examples possessing exact analytic solutions are evaluated numerically for purposes of comparison. Various discretization prescriptions are investigated and found to be equivalent as expected. The numerical method is shown to give accurate results provided the spatial discretization and the time step satisfy certain relationships determined by the drift and the diffusion functions of the nonlinear Fokker-Planck equations.

• Received 15 November 1982

DOI:https://doi.org/10.1103/PhysRevA.27.2663