/SDE-higham

Some basic algorithms for stochastic differential equations in Python/NumPy

Primary LanguagePython

Some basic algorithms for stochastic differential equations in NumPy

Overview

Modified from the MATLAB versions in Higham “An Algorithmic Introduction to Numerical Simulation of Stochastic Differential Equations”, SIAM Review, Vol. 43, No. 3, 2001. http://www.caam.rice.edu/~cox/stoch/dhigham.pdf

Details about the algorithms can be found in the paper.

List of algorithms

FilenameDescription
bpath2.pySimulation of a Brownian path
bpath3.pyFunction along a Brownain path
stint.pyApproximate stochastic integrals
em.pyEuler-Maruyama method on linear SDE
emstrong.pyTest strong convergence of Euler-Maruyama
emweak.pyTest weak convergence of Euler-Maruyama
milstrong.pyTest strong convergence of Milstein method
stab.pyMean-square and asymptotic stability test for E-M
chain.pyTest stochastic chain rule
milstein-3d.pyMilstein’s method applied to a 3D SDE (not part of SIAM paper)