When pricing path-dependent derivatives in multi-factor models (e.g. a basket option under two correlated GBMs, or a Heston model combining a stock SDE and a variance SDE), Monte Carlo simulations require discretising the system of SDEs. The Euler-Maruyama scheme uses only the first-order stochastic Taylor term. The Milstein scheme adds the second-order terms, achieving a higher convergence rate — but in dimension , the second-order terms involve the Lévy area: an iterated stochastic integral that cannot be computed from the BM increments and alone.
This creates a fundamental obstacle: in 1D, the Itô correction is deterministic, so second-order terms reduce to . In , the off-diagonal second-order terms introduce random correction terms (the Lévy areas) that require either simulation (expensive) or bounding (introducing error). Understanding this is essential for anyone implementing or assessing multi-factor Monte Carlo engines.
Learning objectives. By the end of this note you should be able to:
- Derive the first- and second-order stochastic Taylor expansion of an Itô process via iterated application of Itô's formula.
- Define the iterated stochastic integral and the Lévy area .
- Explain why for is not determined by and in dimension .
- State the Euler-Maruyama and Milstein schemes and compare their strong convergence orders in 1D and .
- Identify when Lévy areas matter in practical Monte Carlo.
In this lesson
- The Stochastic Taylor Expansion
- The Lévy Area
- Numerical Implications
- When Lévy Areas Matter (and When They Don't)
Before you start
- Covariation, the Product Rule, and Multiple Brownian Motions
- Derivation — Quadratic Variation of Brownian Motion
The full lesson, with worked examples and exercise solutions, is for members.
Read the full lessonNew to SelfMadeQuant.com? See membership options or explore the curriculum.