Let be independent random variables with the standard normal distribution. The unconstrained Milstein method step is
The noise is additive, so its Milstein correction vanishes. The diagonal, independent noise fields commute, so no cross iterated stochastic integral is required.
Implement each reflecting boundary condition for a diffusion by folding the proposed coordinate back into . One formula that also handles multiple overshoots is
The complete update is
Under the usual smoothness assumptions, the Milstein discretization has strong order of convergence one and weak order of convergence one. The reflection enforces the boundary pathwise.

Articles by others on the same topic (0)

There are currently no matching articles.