Chapter 3 — Itô's Formula
Guiding question
If f f f is a smooth function and W t W_t W t is Brownian motion, how does f ( W t ) f(W_t) f ( W t ) change over time? The
ordinary chain rule says d ( f ∘ W ) = f ′ d W d(f\circ W) = f'\,dW d ( f ∘ W ) = f ′ d W . Why is this wrong, and what replaces it?
Itô's formula is the chain rule of stochastic calculus. It is the single most important computational tool in the entire theory. Every SDE that can be solved explicitly is solved via Itô's formula. Every partial differential equation that arises in option pricing is derived via Itô's formula.
The reason the ordinary chain rule fails is the quadratic variation of Brownian motion: ( d W t ) 2 = d t (dW_t)^2 = dt ( d W t ) 2 = d t . This creates an additional drift term proportional to the second derivative f ′ ′ f'' f ′′ .
3.1 Why the ordinary chain rule fails
For a deterministic differentiable function g : [ 0 , T ] → R g:[0,T]\to\mathbb R g : [ 0 , T ] → R and a smooth function f : R → R f:\mathbb R\to\mathbb R f : R → R , the chain rule gives
f ( g ( T ) ) − f ( g ( 0 ) ) = ∫ 0 T f ′ ( g ( t ) ) g ′ ( t ) d t . f(g(T)) - f(g(0)) = \int_0^T f'(g(t))\,g'(t)\,dt. f ( g ( T )) − f ( g ( 0 )) = ∫ 0 T f ′ ( g ( t )) g ′ ( t ) d t .
In infinitesimal notation: d ( f ∘ g ) = f ′ ( g ) d g d(f\circ g) = f'(g)\,dg d ( f ∘ g ) = f ′ ( g ) d g . The second derivative plays no role because ( d g ) 2 ∼ ( g ′ ( t ) ) 2 ( d t ) 2 = O ( d t 2 ) (dg)^2\sim (g'(t))^2(dt)^2 = O(dt^2) ( d g ) 2 ∼ ( g ′ ( t ) ) 2 ( d t ) 2 = O ( d t 2 ) is negligible.
For stochastic processes, increments of order d t \sqrt{dt} d t change everything. Write the second-order Taylor expansion:
f ( W t + h ) − f ( W t ) = f ′ ( W t ) ( W t + h − W t ) + 1 2 f ′ ′ ( W t ) ( W t + h − W t ) 2 + O ( ∣ W t + h − W t ∣ 3 ) . f(W_{t+h}) - f(W_t)
= f'(W_t)(W_{t+h}-W_t)
+ \tfrac12 f''(W_t)(W_{t+h}-W_t)^2
+ O(|W_{t+h}-W_t|^3). f ( W t + h ) − f ( W t ) = f ′ ( W t ) ( W t + h − W t ) + 2 1 f ′′ ( W t ) ( W t + h − W t ) 2 + O ( ∣ W t + h − W t ∣ 3 ) .
Summing over a partition:
f ( W T ) − f ( W 0 ) ≈ ∑ k f ′ ( W t k ) Δ W k + 1 2 ∑ k f ′ ′ ( W t k ) ( Δ W k ) 2 . f(W_T)-f(W_0)
\approx \sum_k f'(W_{t_k})\Delta W_k
+ \tfrac12\sum_k f''(W_{t_k})(\Delta W_k)^2. f ( W T ) − f ( W 0 ) ≈ k ∑ f ′ ( W t k ) Δ W k + 2 1 k ∑ f ′′ ( W t k ) ( Δ W k ) 2 .
As the partition refines, the first sum converges to the Itô integral ∫ 0 T f ′ ( W t ) d W t \int_0^T f'(W_t)\,dW_t ∫ 0 T f ′ ( W t ) d W t . The second sum involves ( Δ W k ) 2 ≈ Δ t k (\Delta W_k)^2\approx\Delta t_k ( Δ W k ) 2 ≈ Δ t k , so it converges to 1 2 ∫ 0 T f ′ ′ ( W t ) d t \frac12\int_0^T f''(W_t)\,dt 2 1 ∫ 0 T f ′′ ( W t ) d t .
The third-order remainder is negligible because ( Δ W k ) 3 ∼ ( Δ t k ) 3 / 2 (\Delta W_k)^3 \sim (\Delta t_k)^{3/2} ( Δ W k ) 3 ∼ ( Δ t k ) 3/2 and ∑ ( Δ t k ) 3 / 2 ≤ ∣ Π ∣ 1 / 2 ⋅ T → 0 \sum(\Delta t_k)^{3/2} \le |\Pi|^{1/2}\cdot T\to 0 ∑ ( Δ t k ) 3/2 ≤ ∣Π ∣ 1/2 ⋅ T → 0 .
Why the second derivative survives
In classical calculus, the chain rule is first-order: only f ′ f' f ′ appears because ( d g ) 2 = O ( d t 2 ) → 0 (dg)^2 = O(dt^2)\to 0 ( d g ) 2 = O ( d t 2 ) → 0 .
For Brownian motion, ( d W ) 2 = O ( d t ) (dW)^2 = O(dt) ( d W ) 2 = O ( d t ) , which is the same order as d t dt d t itself. So the second derivative is not negligible: it contributes a term of order d t dt d t , exactly the same magnitude as the drift.
This is the geometric content of Itô's formula: Brownian paths are so rough that second derivatives must be included in the linearization.
Itô's formula (basic version)
Let f : R → R f:\mathbb R\to\mathbb R f : R → R be twice continuously differentiable. Then
f ( W T ) − f ( W 0 ) = ∫ 0 T f ′ ( W t ) d W t + 1 2 ∫ 0 T f ′ ′ ( W t ) d t . f(W_T) - f(W_0)
= \int_0^T f'(W_t)\,dW_t
+ \frac12\int_0^T f''(W_t)\,dt. f ( W T ) − f ( W 0 ) = ∫ 0 T f ′ ( W t ) d W t + 2 1 ∫ 0 T f ′′ ( W t ) d t . In differential notation:
d f ( W t ) = f ′ ( W t ) d W t + 1 2 f ′ ′ ( W t ) d t . df(W_t) = f'(W_t)\,dW_t + \tfrac12 f''(W_t)\,dt. df ( W t ) = f ′ ( W t ) d W t + 2 1 f ′′ ( W t ) d t .
Proof via partition and Taylor expansion
Let Π = { 0 = t 0 < t 1 < ⋯ < t n = T } \Pi = \{0=t_0<t_1<\cdots<t_n=T\} Π = { 0 = t 0 < t 1 < ⋯ < t n = T } . By the Taylor expansion to second order:
f ( W t k + 1 ) − f ( W t k ) = f ′ ( W t k ) Δ W k + 1 2 f ′ ′ ( W t k ) ( Δ W k ) 2 + R k , f(W_{t_{k+1}})-f(W_{t_k})
= f'(W_{t_k})\Delta W_k
+ \tfrac12 f''(W_{t_k})(\Delta W_k)^2
+ R_k, f ( W t k + 1 ) − f ( W t k ) = f ′ ( W t k ) Δ W k + 2 1 f ′′ ( W t k ) ( Δ W k ) 2 + R k , where R k = 1 2 [ f ′ ′ ( ξ k ) − f ′ ′ ( W t k ) ] ( Δ W k ) 2 R_k = \frac12[f''(\xi_k)-f''(W_{t_k})](\Delta W_k)^2 R k = 2 1 [ f ′′ ( ξ k ) − f ′′ ( W t k )] ( Δ W k ) 2 for some ξ k \xi_k ξ k between W t k W_{t_k} W t k and W t k + 1 W_{t_{k+1}} W t k + 1 .
Step 1. Summing and telescoping:
f ( W T ) − f ( W 0 ) = ∑ k f ′ ( W t k ) Δ W k + 1 2 ∑ k f ′ ′ ( W t k ) ( Δ W k ) 2 + ∑ k R k . f(W_T)-f(W_0)
= \sum_k f'(W_{t_k})\Delta W_k
+ \tfrac12\sum_k f''(W_{t_k})(\Delta W_k)^2
+ \sum_k R_k. f ( W T ) − f ( W 0 ) = k ∑ f ′ ( W t k ) Δ W k + 2 1 k ∑ f ′′ ( W t k ) ( Δ W k ) 2 + k ∑ R k . Step 2. The first sum converges in L 2 L^2 L 2 to ∫ 0 T f ′ ( W t ) d W t \int_0^T f'(W_t)\,dW_t ∫ 0 T f ′ ( W t ) d W t by definition of the Itô integral.
Step 3. For the second sum, write:
∑ k f ′ ′ ( W t k ) ( Δ W k ) 2 = ∑ k f ′ ′ ( W t k ) Δ t k + ∑ k f ′ ′ ( W t k ) [ ( Δ W k ) 2 − Δ t k ] . \sum_k f''(W_{t_k})(\Delta W_k)^2
= \sum_k f''(W_{t_k})\Delta t_k
+ \sum_k f''(W_{t_k})[(\Delta W_k)^2-\Delta t_k]. k ∑ f ′′ ( W t k ) ( Δ W k ) 2 = k ∑ f ′′ ( W t k ) Δ t k + k ∑ f ′′ ( W t k ) [( Δ W k ) 2 − Δ t k ] . The first part converges to ∫ 0 T f ′ ′ ( W t ) d t \int_0^T f''(W_t)\,dt ∫ 0 T f ′′ ( W t ) d t . For the second part, its L 2 L^2 L 2 norm vanishes as ∣ Π ∣ → 0 |\Pi|\to 0 ∣Π∣ → 0 (the increments ( Δ W k ) 2 − Δ t k (\Delta W_k)^2-\Delta t_k ( Δ W k ) 2 − Δ t k are centered and orthogonal, so their sum's variance is 2 ∑ k ( Δ t k ) 2 ≤ 2 ∣ Π ∣ ⋅ T → 0 2\sum_k(\Delta t_k)^2\le 2|\Pi|\cdot T\to 0 2 ∑ k ( Δ t k ) 2 ≤ 2∣Π∣ ⋅ T → 0 ).
Step 4. The remainder terms ∑ k R k → 0 \sum_k R_k\to 0 ∑ k R k → 0 in L 2 L^2 L 2 by continuity of f ′ ′ f'' f ′′ and the fact that ∑ k ( Δ W k ) 2 → T \sum_k(\Delta W_k)^2\to T ∑ k ( Δ W k ) 2 → T .
Combining, the result holds in L 2 L^2 L 2 .
Itô formula: f(W_t) decomposes into martingale + drift correction
Function: f(x) = x² f(x) = eˣ
0 1 1 Monte Carlo mean Itô formula: E[W_t²] = t Itô: d(W²) = 2W dW + dt. Taking expectations: E[W_t²] = t. Without the correction term dt, classical calculus would predict E[W_t²] = 0 (wrong).
3.3 The symbolic rule ( d W t ) 2 = d t (dW_t)^2 = dt ( d W t ) 2 = d t
Itô's formula is efficiently remembered via the multiplication table:
× \times × d W t dW_t d W t d t dt d t d W t dW_t d W t d t dt d t 0 0 0 d t dt d t 0 0 0 0 0 0
This says: ( d W t ) 2 = d t (dW_t)^2 = dt ( d W t ) 2 = d t , d W t ⋅ d t = 0 dW_t\cdot dt = 0 d W t ⋅ d t = 0 , ( d t ) 2 = 0 (dt)^2 = 0 ( d t ) 2 = 0 . The table encodes the order-of-magnitude argument:
d W t ∼ d t , ( d W t ) 2 ∼ d t , d W t ⋅ d t ∼ d t 3 / 2 → 0 , ( d t ) 2 ∼ d t 2 → 0. dW_t\sim\sqrt{dt},\quad
(dW_t)^2\sim dt,\quad
dW_t\cdot dt\sim dt^{3/2}\to 0,\quad
(dt)^2\sim dt^2\to 0. d W t ∼ d t , ( d W t ) 2 ∼ d t , d W t ⋅ d t ∼ d t 3/2 → 0 , ( d t ) 2 ∼ d t 2 → 0.
Using this table, the chain rule d f = f ′ d W + 1 2 f ′ ′ ( d W ) 2 = f ′ d W + 1 2 f ′ ′ d t df = f'dW + \frac12 f''(dW)^2 = f'dW + \frac12 f''dt df = f ′ d W + 2 1 f ′′ ( d W ) 2 = f ′ d W + 2 1 f ′′ d t follows immediately from the formal second-order Taylor expansion.
3.4 General Itô process
The formula extends to any Itô process, not just Brownian motion.
Itô's formula for Itô processes
Let X t X_t X t satisfy
d X t = b t d t + σ t d W t , dX_t = b_t\,dt + \sigma_t\,dW_t, d X t = b t d t + σ t d W t , where b b b and σ \sigma σ are adapted and satisfy integrability conditions. Let f : R → R f:\mathbb R\to\mathbb R f : R → R be C 2 C^2 C 2 . Then
d f ( X t ) = f ′ ( X t ) d X t + 1 2 f ′ ′ ( X t ) d [ X ] t , df(X_t) = f'(X_t)\,dX_t + \tfrac12 f''(X_t)\,d[X]_t, df ( X t ) = f ′ ( X t ) d X t + 2 1 f ′′ ( X t ) d [ X ] t , where d [ X ] t = σ t 2 d t d[X]_t = \sigma_t^2\,dt d [ X ] t = σ t 2 d t . Expanding:
d f ( X t ) = f ′ ( X t ) b t d t + f ′ ( X t ) σ t d W t + 1 2 f ′ ′ ( X t ) σ t 2 d t . df(X_t)
= f'(X_t)b_t\,dt
+ f'(X_t)\sigma_t\,dW_t
+ \tfrac12 f''(X_t)\sigma_t^2\,dt. df ( X t ) = f ′ ( X t ) b t d t + f ′ ( X t ) σ t d W t + 2 1 f ′′ ( X t ) σ t 2 d t .
The proof follows the same structure, using ( d X t ) 2 = σ t 2 d t (dX_t)^2 = \sigma_t^2\,dt ( d X t ) 2 = σ t 2 d t from the Itô multiplication table.
Itô's formula with time dependence
Let f = f ( t , x ) f = f(t,x) f = f ( t , x ) be C 1 , 2 C^{1,2} C 1 , 2 (once in t t t , twice in x x x ). Then
d f ( t , X t ) = ∂ t f d t + ∂ x f d X t + 1 2 ∂ x x f ( d X t ) 2 , df(t,X_t)
= \partial_t f\,dt
+ \partial_x f\,dX_t
+ \tfrac12\partial_{xx}f\,(dX_t)^2, df ( t , X t ) = ∂ t f d t + ∂ x f d X t + 2 1 ∂ xx f ( d X t ) 2 , which expands to
d f ( t , X t ) = [ ∂ t f + b t ∂ x f + 1 2 σ t 2 ∂ x x f ] d t + σ t ∂ x f d W t . df(t,X_t)
= \left[\partial_t f + b_t\,\partial_x f + \tfrac12\sigma_t^2\,\partial_{xx}f\right]dt
+ \sigma_t\,\partial_x f\,dW_t. df ( t , X t ) = [ ∂ t f + b t ∂ x f + 2 1 σ t 2 ∂ xx f ] d t + σ t ∂ x f d W t .
Three sources of change
The function f ( t , X t ) f(t,X_t) f ( t , X t ) changes due to three effects:
Time drift ∂ t f d t \partial_t f\,dt ∂ t f d t : direct time-dependence of f f f .
First-order X X X term ∂ x f ⋅ d X t \partial_x f\cdot dX_t ∂ x f ⋅ d X t : classical chain rule contribution.
Curvature correction 1 2 ∂ x x f ⋅ σ 2 d t \frac12\partial_{xx}f\cdot\sigma^2\,dt 2 1 ∂ xx f ⋅ σ 2 d t : Brownian quadratic variation creates an extra drift proportional to the second derivative. This is the purely stochastic term.
In the language of differential geometry, the second term is the gradient dotted with the velocity, and the third term is a Laplacian-like correction from the rough geometry of Brownian paths.
3.6 Key examples
Example 1: d ( W t 2 ) d(W_t^2) d ( W t 2 )
Take f ( x ) = x 2 f(x) = x^2 f ( x ) = x 2 , so f ′ = 2 x f' = 2x f ′ = 2 x , f ′ ′ = 2 f'' = 2 f ′′ = 2 :
d ( W t 2 ) = 2 W t d W t + 1 2 ⋅ 2 d t = 2 W t d W t + d t . d(W_t^2) = 2W_t\,dW_t + \tfrac12\cdot 2\,dt = 2W_t\,dW_t + dt. d ( W t 2 ) = 2 W t d W t + 2 1 ⋅ 2 d t = 2 W t d W t + d t .
Integrating:
W T 2 = 2 ∫ 0 T W t d W t + T . W_T^2 = 2\int_0^T W_t\,dW_t + T. W T 2 = 2 ∫ 0 T W t d W t + T .
This recovers the computation of Section 2.10 from a different angle.
Example 2: d ( e W t ) d(e^{W_t}) d ( e W t )
Take f ( x ) = e x f(x) = e^x f ( x ) = e x , so f ′ = f ′ ′ = e x f'=f''=e^x f ′ = f ′′ = e x :
d ( e W t ) = e W t d W t + 1 2 e W t d t . d(e^{W_t}) = e^{W_t}\,dW_t + \tfrac12 e^{W_t}\,dt. d ( e W t ) = e W t d W t + 2 1 e W t d t .
Taking expectations: d d t E [ e W t ] = 1 2 E [ e W t ] \frac{d}{dt}\mathbb E[e^{W_t}] = \tfrac12\mathbb E[e^{W_t}] d t d E [ e W t ] = 2 1 E [ e W t ] , which gives E [ e W t ] = e t / 2 \mathbb E[e^{W_t}] = e^{t/2} E [ e W t ] = e t /2 . This is the standard moment generating function of N ( 0 , t ) N(0,t) N ( 0 , t ) .
Example 3: d ( e θ W t − θ 2 t / 2 ) d(e^{\theta W_t - \theta^2 t/2}) d ( e θ W t − θ 2 t /2 ) — the exponential martingale
Take f ( t , x ) = e θ x − θ 2 t / 2 f(t,x) = e^{\theta x - \theta^2 t/2} f ( t , x ) = e θ x − θ 2 t /2 . Then
∂ t f = − 1 2 θ 2 f , ∂ x f = θ f , ∂ x x f = θ 2 f . \partial_t f = -\tfrac12\theta^2 f,\quad
\partial_x f = \theta f,\quad
\partial_{xx}f = \theta^2 f. ∂ t f = − 2 1 θ 2 f , ∂ x f = θ f , ∂ xx f = θ 2 f .
Itô's formula gives
d f = ( − 1 2 θ 2 f + 1 2 θ 2 f ) d t + θ f d W t = θ f d W t . df = \left(-\tfrac12\theta^2 f + \tfrac12\theta^2 f\right)dt + \theta f\,dW_t = \theta f\,dW_t. df = ( − 2 1 θ 2 f + 2 1 θ 2 f ) d t + θ f d W t = θ f d W t .
Since the drift vanishes, E t ( θ ) = e θ W t − θ 2 t / 2 \mathcal E_t(\theta) = e^{\theta W_t - \theta^2 t/2} E t ( θ ) = e θ W t − θ 2 t /2 is a local martingale (and under integrability conditions, a true martingale). This is the Doléans-Dade exponential or stochastic exponential , central to Girsanov's theorem.
3.7 Geometric Brownian Motion
The most important SDE in finance is:
d S t = μ S t d t + σ S t d W t , S 0 > 0. dS_t = \mu S_t\,dt + \sigma S_t\,dW_t,
\quad S_0 > 0. d S t = μ S t d t + σ S t d W t , S 0 > 0.
Strategy: apply Itô's formula to log S t \log S_t log S t .
Let f ( x ) = log x f(x) = \log x f ( x ) = log x (for x > 0 x > 0 x > 0 ). Then f ′ ( x ) = 1 / x f'(x) = 1/x f ′ ( x ) = 1/ x , f ′ ′ ( x ) = − 1 / x 2 f''(x) = -1/x^2 f ′′ ( x ) = − 1/ x 2 . Itô's formula gives:
d ( log S t ) = 1 S t d S t + 1 2 ( − 1 S t 2 ) ( d S t ) 2 . d(\log S_t)
= \frac{1}{S_t}\,dS_t + \tfrac12\left(-\frac{1}{S_t^2}\right)(dS_t)^2. d ( log S t ) = S t 1 d S t + 2 1 ( − S t 2 1 ) ( d S t ) 2 .
Now ( d S t ) 2 = σ 2 S t 2 d t (dS_t)^2 = \sigma^2 S_t^2\,dt ( d S t ) 2 = σ 2 S t 2 d t , so:
d ( log S t ) = μ S t d t + σ S t d W t S t − 1 2 σ 2 d t = ( μ − 1 2 σ 2 ) d t + σ d W t . d(\log S_t)
= \frac{\mu S_t\,dt + \sigma S_t\,dW_t}{S_t} - \tfrac12\sigma^2\,dt
= \left(\mu - \tfrac12\sigma^2\right)dt + \sigma\,dW_t. d ( log S t ) = S t μ S t d t + σ S t d W t − 2 1 σ 2 d t = ( μ − 2 1 σ 2 ) d t + σ d W t .
The right side is a Brownian motion with constant coefficients. Integrating:
log S t = log S 0 + ( μ − 1 2 σ 2 ) t + σ W t , \log S_t = \log S_0 + \left(\mu - \tfrac12\sigma^2\right)t + \sigma W_t, log S t = log S 0 + ( μ − 2 1 σ 2 ) t + σ W t ,
so
S t = S 0 exp [ ( μ − 1 2 σ 2 ) t + σ W t ] . \boxed{S_t = S_0\exp\!\left[\left(\mu - \tfrac12\sigma^2\right)t + \sigma W_t\right].} S t = S 0 exp [ ( μ − 2 1 σ 2 ) t + σ W t ] .
Log-normal returns and the volatility drag
The expected value is E [ S t ] = S 0 e μ t \mathbb E[S_t] = S_0 e^{\mu t} E [ S t ] = S 0 e μ t (the drift μ \mu μ drives the mean). But the
median of S t S_t S t is S 0 e ( μ − σ 2 / 2 ) t S_0 e^{(\mu-\sigma^2/2)t} S 0 e ( μ − σ 2 /2 ) t , which is smaller. The term − σ 2 / 2 -\sigma^2/2 − σ 2 /2 is the
volatility drag or Itô correction : randomness moves the median below the mean due to Jensen's
inequality applied to the convex exponential function. Higher volatility means slower median growth,
even with the same expected return.
Let X t = ( X t 1 , … , X t d ) \mathbf X_t = (X_t^1,\ldots,X_t^d) X t = ( X t 1 , … , X t d ) with d X t i = b t i d t + σ t i j d W t j dX_t^i = b_t^i\,dt + \sigma_t^{ij}\,dW_t^j d X t i = b t i d t + σ t ij d W t j (summed over j j j ). For f ∈ C 1 , 2 f\in C^{1,2} f ∈ C 1 , 2 ,
d f ( t , X t ) = ∂ t f d t + ∑ i ∂ i f d X t i + 1 2 ∑ i , j ∂ i j f d [ X i , X j ] t . df(t,\mathbf X_t)
= \partial_t f\,dt
+ \sum_i \partial_i f\,dX_t^i
+ \tfrac12\sum_{i,j}\partial_{ij}f\,d[X^i,X^j]_t. df ( t , X t ) = ∂ t f d t + i ∑ ∂ i f d X t i + 2 1 i , j ∑ ∂ ij f d [ X i , X j ] t .
where d [ X i , X j ] t = ∑ k σ t i k σ t j k d t d[X^i,X^j]_t = \sum_k\sigma_t^{ik}\sigma_t^{jk}\,dt d [ X i , X j ] t = ∑ k σ t ik σ t j k d t . The details are in Chapter 8.
3.9 Itô's product rule
As a special case, for two Itô processes X t X_t X t and Y t Y_t Y t :
d ( X t Y t ) = X t d Y t + Y t d X t + d [ X , Y ] t . d(X_t Y_t) = X_t\,dY_t + Y_t\,dX_t + d[X,Y]_t. d ( X t Y t ) = X t d Y t + Y t d X t + d [ X , Y ] t .
The last term is the stochastic correction. Compare with the deterministic product rule (Leibniz rule), which has no such correction.
Derivation from the bilinear identity
Apply Itô's formula to f ( x , y ) = x y f(x,y)=xy f ( x , y ) = x y : then ∂ x y f = 1 \partial_{xy}f = 1 ∂ x y f = 1 , all other second derivatives vanish. The formula gives d ( X Y ) = X d Y + Y d X + d [ X , Y ] d(XY) = X\,dY + Y\,dX + d[X,Y] d ( X Y ) = X d Y + Y d X + d [ X , Y ] .
3.10 Problems
Problem 1: d(W_t^3)
Use Itô's formula to compute d ( W t 3 ) d(W_t^3) d ( W t 3 ) and integrate to find ∫ 0 T W t 2 d W t \int_0^T W_t^2\,dW_t ∫ 0 T W t 2 d W t in terms of W T W_T W T and T T T .
Problem 2: Exponential process
Let X t = e a t + b W t X_t = e^{at + bW_t} X t = e a t + b W t for constants a , b a,b a , b . Find the SDE satisfied by X X X .
Problem 3: Itô's formula for $f(t,x) = x^2 - t$
Show directly that f ( t , W t ) = W t 2 − t f(t,W_t) = W_t^2 - t f ( t , W t ) = W t 2 − t is a local martingale. What is its SDE?
Problem 4: Integration by parts
Use the Itô product rule to show that
W T 2 = 2 ∫ 0 T W t d W t + T . W_T^2 = 2\int_0^T W_t\,dW_t + T. W T 2 = 2 ∫ 0 T W t d W t + T .
Problem 5: Lognormal moments
For GBM with S t = S 0 exp ( ( μ − σ 2 / 2 ) t + σ W t ) S_t = S_0\exp((\mu-\sigma^2/2)t + \sigma W_t) S t = S 0 exp (( μ − σ 2 /2 ) t + σ W t ) , compute E [ S t n ] \mathbb E[S_t^n] E [ S t n ] for integer n ≥ 1 n\ge 1 n ≥ 1 .
3.11 Selected solutions
Solution 1
With f ( x ) = x 3 f(x) = x^3 f ( x ) = x 3 , f ′ = 3 x 2 f'=3x^2 f ′ = 3 x 2 , f ′ ′ = 6 x f''=6x f ′′ = 6 x :
d ( W t 3 ) = 3 W t 2 d W t + 1 2 ⋅ 6 W t d t = 3 W t 2 d W t + 3 W t d t . d(W_t^3) = 3W_t^2\,dW_t + \tfrac12\cdot 6W_t\,dt = 3W_t^2\,dW_t + 3W_t\,dt. d ( W t 3 ) = 3 W t 2 d W t + 2 1 ⋅ 6 W t d t = 3 W t 2 d W t + 3 W t d t . Integrating: W T 3 = 3 ∫ 0 T W t 2 d W t + 3 ∫ 0 T W t d t W_T^3 = 3\int_0^T W_t^2\,dW_t + 3\int_0^T W_t\,dt W T 3 = 3 ∫ 0 T W t 2 d W t + 3 ∫ 0 T W t d t , so
∫ 0 T W t 2 d W t = 1 3 W T 3 − ∫ 0 T W t d t . \int_0^T W_t^2\,dW_t = \tfrac13 W_T^3 - \int_0^T W_t\,dt. ∫ 0 T W t 2 d W t = 3 1 W T 3 − ∫ 0 T W t d t . The remaining integral ∫ 0 T W t d t \int_0^T W_t\,dt ∫ 0 T W t d t is a Gaussian with mean 0 0 0 and variance T 3 / 3 T^3/3 T 3 /3 .
Solution 2
Let f ( t , x ) = e a t + b x f(t,x) = e^{at+bx} f ( t , x ) = e a t + b x . Then ∂ t f = a f \partial_t f = af ∂ t f = a f , ∂ x f = b f \partial_x f = bf ∂ x f = b f , ∂ x x f = b 2 f \partial_{xx}f = b^2f ∂ xx f = b 2 f . Itô's formula:
d X t = a f d t + b f d W t + 1 2 b 2 f d t = ( a + 1 2 b 2 ) X t d t + b X t d W t . dX_t = af\,dt + bf\,dW_t + \tfrac12 b^2 f\,dt = (a+\tfrac12 b^2)X_t\,dt + bX_t\,dW_t. d X t = a f d t + b f d W t + 2 1 b 2 f d t = ( a + 2 1 b 2 ) X t d t + b X t d W t .
Solution 3
With f ( t , x ) = x 2 − t f(t,x)=x^2-t f ( t , x ) = x 2 − t : ∂ t f = − 1 \partial_t f = -1 ∂ t f = − 1 , ∂ x f = 2 x \partial_x f=2x ∂ x f = 2 x , ∂ x x f = 2 \partial_{xx}f=2 ∂ xx f = 2 . Itô:
d ( W t 2 − t ) = ( − 1 + 1 2 ⋅ 2 ) d t + 2 W t d W t = 2 W t d W t . d(W_t^2-t) = (-1+\tfrac12\cdot 2)\,dt + 2W_t\,dW_t = 2W_t\,dW_t. d ( W t 2 − t ) = ( − 1 + 2 1 ⋅ 2 ) d t + 2 W t d W t = 2 W t d W t . The drift vanishes: W t 2 − t W_t^2-t W t 2 − t is a local martingale with zero mean.
Solution 5
Since S t n = S 0 n exp ( n ( μ − σ 2 / 2 ) t + n σ W t ) S_t^n = S_0^n\exp(n(\mu-\sigma^2/2)t + n\sigma W_t) S t n = S 0 n exp ( n ( μ − σ 2 /2 ) t + nσ W t ) , and e n σ W t e^{n\sigma W_t} e nσ W t has moment generating function e n 2 σ 2 t / 2 e^{n^2\sigma^2 t/2} e n 2 σ 2 t /2 :
E [ S t n ] = S 0 n exp [ ( n μ − n 2 σ 2 + n 2 2 σ 2 ) t ] = S 0 n exp [ ( n μ + n ( n − 1 ) 2 σ 2 ) t ] . \mathbb E[S_t^n] = S_0^n\exp\!\left[\left(n\mu - \tfrac{n}{2}\sigma^2 + \tfrac{n^2}{2}\sigma^2\right)t\right]
= S_0^n\exp\!\left[\left(n\mu + \tfrac{n(n-1)}{2}\sigma^2\right)t\right]. E [ S t n ] = S 0 n exp [ ( n μ − 2 n σ 2 + 2 n 2 σ 2 ) t ] = S 0 n exp [ ( n μ + 2 n ( n − 1 ) σ 2 ) t ] . For n = 1 n=1 n = 1 : E [ S t ] = S 0 e μ t \mathbb E[S_t] = S_0 e^{\mu t} E [ S t ] = S 0 e μ t . For n = 2 n=2 n = 2 : E [ S t 2 ] = S 0 2 e ( 2 μ + σ 2 ) t \mathbb E[S_t^2] = S_0^2 e^{(2\mu+\sigma^2)t} E [ S t 2 ] = S 0 2 e ( 2 μ + σ 2 ) t .
3.12 What this chapter gives you
Itô's formula is the engine of stochastic calculus:
d f ( t , X t ) = ∂ t f d t ⏟ time + ∂ x f d X t ⏟ first order + 1 2 ∂ x x f d [ X ] t ⏟ It o ˆ correction . df(t,X_t)
= \underbrace{\partial_t f\,dt}_{\text{time}}
+ \underbrace{\partial_x f\,dX_t}_{\text{first order}}
+ \underbrace{\tfrac12\partial_{xx}f\,d[X]_t}_{\text{Itô correction}}. df ( t , X t ) = time ∂ t f d t + first order ∂ x f d X t + It o ˆ correction 2 1 ∂ xx f d [ X ] t .
The third term is unique to stochastic calculus and arises from ( d W ) 2 = d t (dW)^2 = dt ( d W ) 2 = d t . It:
creates the volatility drag in GBM;
is the source of the 1 2 σ 2 \frac12\sigma^2 2 1 σ 2 term in the Black–Scholes PDE;
is what makes stochastic calculus fundamentally different from classical calculus.
The next chapter uses Itô's formula as the main tool for stating and solving stochastic differential equations.