Bernal Martínez, Francisco Manuel2021-03-102021-03-102019-05Bernal, F. (2018). An Implementation of Milstein’s Method for General Bounded Diffusions. Journal of Scientific Computing, 79(2), 867–890.0885-7474https://hdl.handle.net/10016/32092Despite its generality and powerful convergence properties, Milstein’s method for functionals of spatially bounded stochastic differential equations is widely regarded as difficult to implement. This has likely prevented it from being utilised in applications. In this paper, we design and analyse in detail one such implementation. The presented method turns out to be on par with other, popular schemes in terms of computational cost—but with a (nearly) linear weak convergence rate under the usual smoothness requirements on coefficients and boundary. Two byproducts of theoretical interest are a new, non-standard rank-one update formula, and a connection between numerics of bounded diffusions and Eikonal equations. Three examples are worked out, confirming the accuracy and robustness of the method.24eng© 2018, Springer Science Business Media, LLC, part of Springer NatureStopped diffusionReflected diffusionFeynman-kac formulaStochastic numericsRank-one updates to matrix decompositionsEikonal equationresearch articleMatemáticashttps://doi.org/10.1007/s10915-018-0884-6open access8672890Journal of Scientific Computing79AR/0000026502