Itô's Lemma and Stochastic Calculus › Worked Examples

Michael BrenndoerferNovember 23, 202510 min read
Quantitative FinanceAll chapters
Part 5 of 15Stochastic Calculus and Pricing Theory
Chapter 27Itô's Lemma and Stochastic Calculus › Worked Examples
Reading section 5 of 11 · Worked Examples · 10 minSection 5 of 11 · 10 min

Book contents

Section 5 of 11 · 10 min read

Itô's Lemma and Stochastic Calculus › Worked ExamplesOne section · Sections in this chapter5 / 11 ↓
Reading tools Definitions & inline links

Choose your expertise level to adjust how many terms are explained. Beginners see more tooltips, experts see fewer to maintain reading flow. Hover over underlined terms for instant definitions.

Article links

Make inline references clickable

Let's apply Itô's Lemma to several important examples that arise frequently in financial modeling. These examples illustrate both the mechanics of applying the formula and the surprising results that emerge from stochastic calculus.

Example 1: Squared Brownian Motion

Consider f(Wt)=Wt2f(W_t) = W_t^2 where WtW_t is standard Brownian motion. Here dWt=dWtdW_t = dW_t (the process is Brownian motion itself with μ=0\mu = 0 and σ=1\sigma = 1). This simple example clearly demonstrates how Itô's Lemma differs from ordinary calculus.

Computing the partial derivatives:

∂f∂x=2Wt,∂2f∂x2=2,∂f∂t=0\frac{\partial f}{\partial x} = 2W_t, \quad \frac{\partial^2 f}{\partial x^2} = 2, \quad \frac{\partial f}{\partial t} = 0

where:

  • f(x,t)=x2f(x, t) = x^2: the function being differentiated

  • WtW_t: the stochastic process (Brownian motion)

Applying Itô's Lemma with μ=0\mu = 0 and σ=1\sigma = 1:

d(Wt2)=(0+0+12⋅12⋅2)dt+1⋅2Wt dWt(substitute partials into Itoˆ formula)=dt+2Wt dWt(simplify)\begin{aligned} d(W_t^2) &= \left( 0 + 0 + \frac{1}{2} \cdot 1^2 \cdot 2 \right) dt + 1 \cdot 2W_t \, dW_t && \text{(substitute partials into Itô formula)} \\ &= dt + 2W_t \, dW_t && \text{(simplify)} \end{aligned}

where:

  • d(Wt2)d(W_t^2): differential of the squared Brownian motion

  • dtdt: contribution from the Itô correction term

  • 2Wt dWt2W_t \, dW_t: contribution from the standard chain rule

Notice that the ordinary chain rule would give d(Wt2)=2Wt dWtd(W_t^2) = 2W_t \, dW_t, missing the dtdt term entirely. The extra dtdt term comes from the Itô correction and reflects the fact that Brownian motion is constantly fluctuating. Even when Wt=0W_t = 0, the squared process is increasing on average. This may seem counterintuitive at first, but it makes sense: when the Brownian motion passes through zero, it doesn't stay there. The constant oscillations around zero still contribute positive values to the squared process.

We can verify this result and gain further insight by integrating both sides:

WT2=∫0Tdt+2∫0TWt dWt=T+2∫0TWt dWt\begin{aligned} W_T^2 &= \int_0^T dt + 2\int_0^T W_t \, dW_t \\ &= T + 2\int_0^T W_t \, dW_t \end{aligned}

where:

  • WT2W_T^2: terminal squared value of Brownian motion

  • TT: total time horizon (integral of the deterministic drift dtdt)

  • ∫0TWt dWt\int_0^T W_t \, dW_t: stochastic integral term

  • WTW_T: value of Brownian motion at time TT

Rearranging gives us an explicit formula for the stochastic integral:

∫0TWt dWt=12(WT2−T)\int_0^T W_t \, dW_t = \frac{1}{2}(W_T^2 - T)

where:

  • ∫0TWt dWt\int_0^T W_t \, dW_t: stochastic integral of Brownian motion with respect to itself

  • WT2W_T^2: squared terminal value of Brownian motion

Compare this to ordinary calculus, where ∫0Tx dx=12T2\int_0^T x \, dx = \frac{1}{2}T^2. The −T/2-T/2 correction term appears because stochastic integrals behave differently from ordinary integrals. This example showcases how the Itô correction manifests in integral form.

Out[5]:
Visualization
Line chart of 50 muted squared Brownian motion paths with an emphasized mean path and dashed theoretical expectation line growing linearly from 0 to 2.
Sample paths of squared Brownian motion $W_t^2$ and their mean behavior. The muted lines represent 50 individual realizations, while the emphasized mean tracks the theoretical expectation $t$ (dashed line). This linear growth in expectation, despite $W_t$ averaging to zero, illustrates the Itô correction term $dt$.

Example 2: Log of Geometric Brownian Motion

This is the most important example for finance because it connects the geometric Brownian motion model of asset prices to the lognormal distribution. Suppose StS_t follows geometric Brownian motion:

dSt=μSt dt+σSt dWtdS_t = \mu S_t \, dt + \sigma S_t \, dW_t

where:

  • StS_t: asset price

  • μ\mu: drift parameter

  • σ\sigma: volatility parameter

  • dWtdW_t: Brownian increment

Let f(St)=ln⁡(St)f(S_t) = \ln(S_t). We want to find the dynamics of the log-price, which will reveal important properties about the distribution of future prices. Computing derivatives with respect to SS:

∂f∂S=1St,∂2f∂S2=−1St2,∂f∂t=0\frac{\partial f}{\partial S} = \frac{1}{S_t}, \quad \frac{\partial^2 f}{\partial S^2} = -\frac{1}{S_t^2}, \quad \frac{\partial f}{\partial t} = 0

where:

  • f(S)f(S): natural logarithm function

  • StS_t: geometric Brownian motion process

Notice that the second derivative is negative, indicating that the logarithm function is concave. This concavity will cause the Itô correction to subtract from the drift rather than add to it.

Applying Itô's Lemma with μS=μSt\mu_S = \mu S_t and σS=σSt\sigma_S = \sigma S_t:

d(ln⁡St)=(0+μStSt+12σ2St2⋅(−1St2))dt+σStSt dWt(apply Itoˆ’s Lemma)=(μ−12σ2)dt+σ dWt(simplify terms)\begin{aligned} d(\ln S_t) &= \left( 0 + \frac{\mu S_t}{S_t} + \frac{1}{2} \sigma^2 S_t^2 \cdot \left( -\frac{1}{S_t^2} \right) \right) dt + \frac{\sigma S_t}{S_t} \, dW_t && \text{(apply Itô's Lemma)} \\ &= \left( \mu - \frac{1}{2}\sigma^2 \right) dt + \sigma \, dW_t && \text{(simplify terms)} \end{aligned}

where:

  • d(ln⁡St)d(\ln S_t): differential of the log-price process

  • μ−12σ2\mu - \frac{1}{2}\sigma^2: drift of the log-price

  • σ\sigma: volatility of the log-price

This is an main result. Even though StS_t follows geometric Brownian motion (a complicated multiplicative process with state-dependent coefficients), its logarithm follows simple arithmetic Brownian motion with constant drift μ−12σ2\mu - \frac{1}{2}\sigma^2 and constant volatility σ\sigma. This step to additive dynamics makes GBM tractable.

Integrating from 00 to TT:

ln⁡ST−ln⁡S0=(μ−12σ2)T+σWTln⁡STS0=(μ−12σ2)T+σWT\begin{aligned} \ln S_T - \ln S_0 &= \left( \mu - \frac{1}{2}\sigma^2 \right) T + \sigma W_T \\ \ln \frac{S_T}{S_0} &= \left( \mu - \frac{1}{2}\sigma^2 \right) T + \sigma W_T \end{aligned}

where:

  • ST,S0S_T, S_0: asset prices at times TT and 00

  • μ,σ\mu, \sigma: drift and volatility parameters

  • WTW_T: Brownian motion value at TT

Since WT∼N(0,T)W_T \sim \mathcal{N}(0, T), we can immediately determine the distribution of log-returns:

ln⁡STS0∼N((μ−12σ2)T,σ2T)\ln \frac{S_T}{S_0} \sim \mathcal{N}\left( \left( \mu - \frac{1}{2}\sigma^2 \right) T, \sigma^2 T \right)

where:

  • N(⋅)\mathcal{N}(\cdot): normal distribution

  • ln⁡(ST/S0)\ln(S_T/S_0): log-return over the period [0,T][0, T]

This proves that log-returns are normally distributed under the GBM model. Taking the exponential of both sides yields the explicit solution to the GBM SDE:

ST=S0exp⁡((μ−12σ2)T+σWT)S_T = S_0 \exp\left( \left( \mu - \frac{1}{2}\sigma^2 \right) T + \sigma W_T \right)

where:

  • STS_T: asset price at time TT

  • S0S_0: initial asset price

  • μ\mu: expected return parameter

  • σ\sigma: volatility parameter

  • WTW_T: value of Brownian motion at time TT

The term −12σ2-\frac{1}{2}\sigma^2 is called the Itô drift correction. It explains why the expected log-return (μ−12σ2)(\mu - \frac{1}{2}\sigma^2) is less than the expected percentage return μ\mu. This correction arises from the concavity of the logarithm function combined with the volatility of the underlying process. It carries into portfolio theory and long-term wealth accumulation, as we shall see in later chapters.

Example 3: The Exponential Martingale

Consider Yt=eσWt−12σ2tY_t = e^{\sigma W_t - \frac{1}{2}\sigma^2 t}. This exponential form appears frequently in risk-neutral pricing and change of measure arguments. The specific form of the exponent, with its −12σ2t-\frac{1}{2}\sigma^2 t term, is carefully chosen to produce a special property.

Let g(Wt,t)=eσWt−12σ2tg(W_t, t) = e^{\sigma W_t - \frac{1}{2}\sigma^2 t}, or equivalently, g(w,t)=eσw−12σ2tg(w, t) = e^{\sigma w - \frac{1}{2}\sigma^2 t}.

Computing derivatives:

∂g∂t=−12σ2eσw−12σ2t=−12σ2Yt∂g∂w=σeσw−12σ2t=σYt∂2g∂w2=σ2eσw−12σ2t=σ2Yt\begin{aligned} \frac{\partial g}{\partial t} &= -\frac{1}{2}\sigma^2 e^{\sigma w - \frac{1}{2}\sigma^2 t} = -\frac{1}{2}\sigma^2 Y_t \\ \frac{\partial g}{\partial w} &= \sigma e^{\sigma w - \frac{1}{2}\sigma^2 t} = \sigma Y_t \\ \frac{\partial^2 g}{\partial w^2} &= \sigma^2 e^{\sigma w - \frac{1}{2}\sigma^2 t} = \sigma^2 Y_t \end{aligned}

where:

  • g(w,t)g(w, t): the function defining the martingale

  • YtY_t: value of the exponential martingale

  • σ\sigma: volatility parameter

Applying Itô's Lemma (with dWt=dWtdW_t = dW_t, so μ=0\mu = 0, σW=1\sigma_W = 1):

dYt=(−12σ2Yt+0+12⋅12⋅σ2Yt)dt+1⋅σYt dWt(apply Itoˆ’s Lemma)=σYt dWt(drift terms cancel)\begin{aligned} dY_t &= \left( -\frac{1}{2}\sigma^2 Y_t + 0 + \frac{1}{2} \cdot 1^2 \cdot \sigma^2 Y_t \right) dt + 1 \cdot \sigma Y_t \, dW_t && \text{(apply Itô's Lemma)} \\ &= \sigma Y_t \, dW_t && \text{(drift terms cancel)} \end{aligned}

where:

  • dYtdY_t: change in the martingale process

  • σYt\sigma Y_t: volatility term proportional to the process level

The process YtY_t has zero drift! This means YtY_t is a martingale, which is a process whose expected future value equals its current value. The careful choice of −12σ2t-\frac{1}{2}\sigma^2 t in the exponent exactly cancels the Itô correction that would otherwise appear from the convexity of the exponential function. This cancellation is not coincidental; it is precisely engineered to eliminate the drift.

In the next chapter, this exponential martingale provides the mathematical mechanism for changing from the real-world probability measure to the risk-neutral measure, which is essential for derivative pricing.

Example 4: Product of Two Itô Processes

Sometimes we need the dynamics of a product XtYtX_t Y_t, for example when computing the dynamics of wealth or hedged portfolio value. The product rule for Itô processes differs from the ordinary product rule by an additional term:

d(XtYt)=XtdYt+YtdXt+dXtdYtd(X_t Y_t) = X_t dY_t + Y_t dX_t + dX_t dY_t

where:

  • d(XtYt)d(X_t Y_t): differential of the product

  • Xt,YtX_t, Y_t: two Itô processes

  • dXt dYtdX_t \, dY_t: quadratic covariation term

The extra term dXt dYtdX_t \, dY_t has no analog in ordinary calculus. It arises from the cross-term in the Taylor expansion and survives because of the non-zero quadratic variation of Brownian motion.

Using the multiplication rules that encode the behavior of infinitesimal products:

  • dt⋅dt=0dt \cdot dt = 0 (second-order in time, negligible)

  • dt⋅dWt=0dt \cdot dW_t = 0 (mixed term, order (Δt)3/2(\Delta t)^{3/2}, negligible)

  • dWt⋅dWt=dtdW_t \cdot dW_t = dt (quadratic variation of Brownian motion)

If dXt=μXdt+σXdWtdX_t = \mu_X dt + \sigma_X dW_t and dYt=μYdt+σYdWtdY_t = \mu_Y dt + \sigma_Y dW_t (driven by the same Brownian motion), then:

dXt⋅dYt=σXσYdtdX_t \cdot dY_t = \sigma_X \sigma_Y dt

where:

  • σX,σY\sigma_X, \sigma_Y: diffusion coefficients of XX and YY

The full product rule becomes:

d(XtYt)=Xt(μYdt+σYdWt)+Yt(μXdt+σXdWt)+σXσYdt(substitute differentials)=(μXYt+μYXt+σXσY)dt+(σXYt+σYXt)dWt(group drift and diffusion terms)\begin{aligned} d(X_t Y_t) &= X_t (\mu_Y dt + \sigma_Y dW_t) + Y_t (\mu_X dt + \sigma_X dW_t) + \sigma_X \sigma_Y dt && \text{(substitute differentials)} \\ &= (\mu_X Y_t + \mu_Y X_t + \sigma_X \sigma_Y) dt + (\sigma_X Y_t + \sigma_Y X_t) dW_t && \text{(group drift and diffusion terms)} \end{aligned}

where:

  • d(XtYt)d(X_t Y_t): differential of the product process

  • (σXYt+σYXt)dWt(\sigma_X Y_t + \sigma_Y X_t) dW_t: interaction term from covariance term σXσY\sigma_X \sigma_Y in the drift represents the correlation between the two processes. When both processes are driven by the same Brownian motion and both have positive diffusion coefficients, this term adds to the drift of the product. This effect captures the idea that correlated random movements reinforce each other when considering the product.