Quant Sandbox
DerivativesPortfolio & RiskFixed IncomeCorporate FinanceStochastic ModelsGuide
Guide
Derivatives pricing
01Forwards, calls and puts02Combining positions03Arbitrage and the forward price04Put–call parity05The one-step binomial model06Risk-neutral probabilities07Multi-step trees and backward induction08American options and early exercise09Delta and hedging10Geometric Brownian motion11The Black–Scholes formula12The Greeks: gamma, vega, theta and rho13Implied volatility and the volatility smile14Stochastic volatility and the Heston model15Static replication: Breeden–Litzenberger and Carr–Madan16Exotic options: barriers, Asians and digitals
Portfolio & risk
01Returns and volatility02Diversification and correlation03The efficient frontier and the capital market line04Risk aversion and the optimal complete portfolio05CAPM, beta and the security market line06Performance measures: Sharpe, Treynor, Jensen and information ratio07Factor models and systematic risk08Value at Risk and expected shortfall
Fixed income
01Bond prices and yields02Duration and convexity03Immunization and duration matching04Spot rates, forward rates and the yield curve05Bond futures and the cheapest-to-deliver06Interest rate swaps07Callable bonds and negative convexity08Credit risk: default probabilities and credit spreads
Corporate finance
01Interest, compounding and discounting02Annuities, perpetuities and loan amortization03Net present value and internal rate of return04The weighted average cost of capital05Capital structure and the Modigliani–Miller theorems06Discounted cash flow valuation07Relative valuation with multiples08Real options: the value of flexibility
Stochastic models
01Random walks and Brownian motion02Probability spaces, σ-algebras and random variables03Filtrations, adapted processes and stopping times04Conditional expectation and martingales05Properties of Brownian motion06The reflection principle and first passage times07The Brownian bridge08The Itô integral09Itô's lemma10Kolmogorov equations and the Feynman–Kac formula11Change of measure and Girsanov's theorem12Monte Carlo option pricing13Simulating SDEs: the Euler–Maruyama scheme14Mean reversion: Ornstein–Uhlenbeck and Vasicek15Jumps and the Merton jump-diffusion model
Knowledge graph

© 2026 Quant Sandbox · A study aid for quantitative finance

Privacy & dataNo cookies

This site does not use cookies and asks for no personal data. It keeps a few small settings in your browser's local storage — your chosen theme and language, the sidebar layout, and your progress and answers in the Guides — so it remembers them on your next visit. This data never leaves your device; clearing your browser storage removes it. Inputs you change can also be captured in the page URL, so you can share or bookmark a specific setup. The site uses Vercel Web Analytics, a privacy-friendly, cookieless service that records only anonymous, aggregated usage (such as page views) without identifying you or tracking you across other sites.

Guide/Stochastic models

The Itô integral

Lesson 8 of 15 · 1 min read

This lesson defines the stochastic integral ∫0tHs dWs\int_0^t H_s\,dW_s∫0t​Hs​dWs​, which gives meaning to the expression dX=a dt+b dWdX = a\,dt + b\,dWdX=adt+bdW used for asset prices. It builds on "Filtrations, adapted processes and stopping times" (adapted integrands), "Conditional expectation and martingales" (the martingale property) and "Properties of Brownian motion" (infinite total variation and quadratic variation [W]t=t[W]_t = t[W]t​=t). The next lesson, "Itô's lemma", is the chain rule for these integrals.

Why a Riemann–Stieltjes integral fails

For a deterministic function ggg, the Riemann–Stieltjes integral ∫0tf dg\int_0^t f\,dg∫0t​fdg is the limit of sums ∑if(si) (g(ti+1)−g(ti))\sum_i f(s_i)\,(g(t_{i+1}) - g(t_i))∑i​f(si​)(g(ti+1​)−g(ti​)) over partitions 0=t0<t1<⋯<tn=t0 = t_0 < t_1 < \dots < t_n = t0=t0​<t1​<⋯<tn​=t whose mesh max⁡i(ti+1−ti)\max_i (t_{i+1} - t_i)maxi​(ti+1​−ti​) tends to zero, with any evaluation points si∈[ti,ti+1]s_i \in [t_i, t_{i+1}]si​∈[ti​,ti+1​]. A classical theorem says that this limit exists for every continuous fff exactly when ggg has finite total variation. Brownian paths have infinite total variation on every interval, almost surely, so the pathwise construction is not available.

The failure is visible in the simplest case f=g=Wf = g = Wf=g=W. Write ΔWi=Wti+1−Wti\Delta W_i = W_{t_{i+1}} - W_{t_i}ΔWi​=Wti+1​​−Wti​​ and compare right-point and left-point sums:

∑iWti+1 ΔWi−∑iWti ΔWi=∑i(ΔWi)2⟶t.\sum_{i} W_{t_{i+1}}\,\Delta W_i - \sum_{i} W_{t_i}\,\Delta W_i = \sum_{i} (\Delta W_i)^2 \longrightarrow t.i∑​Wti+1​​ΔWi​−i∑​Wti​​ΔWi​=i∑​(ΔWi​)2⟶t.

The difference tends to the quadratic variation ttt, not to zero, so the answer depends on where in each interval the integrand is evaluated. A definition has to fix the evaluation point. Itô's choice is the left endpoint.

Simple integrands

Fix a probability space (Ω,F,P)(\Omega, \mathcal{F}, \mathbb{P})(Ω,F,P) with a filtration (Ft)t≥0(\mathcal{F}_t)_{t \ge 0}(Ft​)t≥0​ and a Brownian motion WWW with respect to it: WWW is adapted and each increment Wu−WsW_u - W_sWu​−Ws​, u>su > su>s, is independent of Fs\mathcal{F}_sFs​. A process HHH is simple if there are times 0=t0<t1<⋯<tn=T0 = t_0 < t_1 < \dots < t_n = T0=t0​<t1​<⋯<tn​=T and random variables ξ0,…,ξn−1\xi_0, \dots, \xi_{n-1}ξ0​,…,ξn−1​ with ξi\xi_iξi​ bounded and Fti\mathcal{F}_{t_i}Fti​​-measurable, such that

Hs=∑i=0n−1ξi 1(ti, ti+1](s).H_s = \sum_{i=0}^{n-1} \xi_i\,\mathbf{1}_{(t_i,\, t_{i+1}]}(s).Hs​=i=0∑n−1​ξi​1(ti​,ti+1​]​(s).

The value on (ti,ti+1](t_i, t_{i+1}](ti​,ti+1​] is fixed at the left endpoint, using only the information available there. For such HHH the Itô integral is defined directly as a finite sum, for tk≤t≤tk+1t_k \le t \le t_{k+1}tk​≤t≤tk+1​:

∫0tHs dWs=∑i=0k−1ξi (Wti+1−Wti)+ξk (Wt−Wtk).\int_0^t H_s\,dW_s = \sum_{i=0}^{k-1} \xi_i\,\big(W_{t_{i+1}} - W_{t_i}\big) + \xi_k\,\big(W_t - W_{t_k}\big).∫0t​Hs​dWs​=i=0∑k−1​ξi​(Wti+1​​−Wti​​)+ξk​(Wt​−Wtk​​).

In trading terms, ξi\xi_iξi​ is the number of shares held over (ti,ti+1](t_i, t_{i+1}](ti​,ti+1​], chosen at time tit_iti​, and the integral is the accumulated gain.

The Itô isometry

Theorem. For a simple HHH,

E ⁣[(∫0THs dWs)2]=E ⁣[∫0THs2 ds].\mathbb{E}\!\left[\Big(\int_0^T H_s\,dW_s\Big)^{2}\right] = \mathbb{E}\!\left[\int_0^T H_s^2\,ds\right].E[(∫0T​Hs​dWs​)2]=E[∫0T​Hs2​ds].

Proof. Let I=∑iξi ΔWiI = \sum_i \xi_i\,\Delta W_iI=∑i​ξi​ΔWi​ and expand the square into diagonal and cross terms. On the diagonal, ξi2\xi_i^2ξi2​ is Fti\mathcal{F}_{t_i}Fti​​-measurable and ΔWi\Delta W_iΔWi​ is independent of Fti\mathcal{F}_{t_i}Fti​​ with variance ti+1−tit_{i+1} - t_iti+1​−ti​, so by the tower property

E[ξi2(ΔWi)2]=E[ξi2 E[(ΔWi)2∣Fti]]=E[ξi2] (ti+1−ti).\mathbb{E}\big[\xi_i^2 (\Delta W_i)^2\big] = \mathbb{E}\Big[\xi_i^2\,\mathbb{E}\big[(\Delta W_i)^2 \mid \mathcal{F}_{t_i}\big]\Big] = \mathbb{E}[\xi_i^2]\,(t_{i+1} - t_i).E[ξi2​(ΔWi​)2]=E[ξi2​E[(ΔWi​)2∣Fti​​]]=E[ξi2​](ti+1​−ti​).

For a cross term with i<ji < ji<j, the product ξiξjΔWi\xi_i \xi_j \Delta W_iξi​ξj​ΔWi​ is Ftj\mathcal{F}_{t_j}Ftj​​-measurable, while E[ΔWj∣Ftj]=0\mathbb{E}[\Delta W_j \mid \mathcal{F}_{t_j}] = 0E[ΔWj​∣Ftj​​]=0; hence E[ξiξjΔWiΔWj]=0\mathbb{E}[\xi_i \xi_j \Delta W_i \Delta W_j] = 0E[ξi​ξj​ΔWi​ΔWj​]=0. Summing the diagonal gives E[I2]=∑iE[ξi2](ti+1−ti)=E[∫0THs2 ds]\mathbb{E}[I^2] = \sum_i \mathbb{E}[\xi_i^2](t_{i+1} - t_i) = \mathbb{E}\big[\int_0^T H_s^2\,ds\big]E[I2]=∑i​E[ξi2​](ti+1​−ti​)=E[∫0T​Hs2​ds]. □\square□

Adaptedness was used twice: to pull ξi\xi_iξi​ out of the conditional expectation, and to make the cross terms vanish. With an integrand that peeks into its own interval both steps break.

Itô isometry: E[(∫0TH dW)2]=E[∫0TH2 ds]\mathbb{E}\big[(\int_0^T H\,dW)^2\big] = \mathbb{E}\big[\int_0^T H^2\,ds\big]E[(∫0T​HdW)2]=E[∫0T​H2ds]. It turns a question about a random integral into an ordinary dsdsds-integral.

Extension to general integrands

Let LT2\mathcal{L}^2_TLT2​ be the set of progressively measurable processes (adapted and jointly measurable in (s,ω)(s, \omega)(s,ω)) with E[∫0THs2 ds]<∞\mathbb{E}\big[\int_0^T H_s^2\,ds\big] < \inftyE[∫0T​Hs2​ds]<∞. Every H∈LT2H \in \mathcal{L}^2_TH∈LT2​ is the limit of simple processes H(n)H^{(n)}H(n) in the norm ∥H∥2=E[∫0THs2 ds]\|H\|^2 = \mathbb{E}\big[\int_0^T H_s^2\,ds\big]∥H∥2=E[∫0T​Hs2​ds]. By the isometry, the integrals I(H(n))I(H^{(n)})I(H(n)) form a Cauchy sequence in L2(Ω)L^2(\Omega)L2(Ω), so they converge, and the limit does not depend on the approximating sequence. This limit is the Itô integral of HHH. The isometry therefore holds for all of LT2\mathcal{L}^2_TLT2​. Integrands with only ∫0THs2 ds<∞\int_0^T H_s^2\,ds < \infty∫0T​Hs2​ds<∞ almost surely can be handled by stopping (localisation); the result is then a local martingale rather than a martingale.

Properties

For H,K∈LT2H, K \in \mathcal{L}^2_TH,K∈LT2​ and constants a,ba, ba,b, the process It=∫0tHs dWsI_t = \int_0^t H_s\,dW_sIt​=∫0t​Hs​dWs​ has the following properties. Each one is checked on simple integrands and passes to the limit.

  • Linearity: ∫(aH+bK) dW=a∫H dW+b∫K dW\int (aH + bK)\,dW = a\int H\,dW + b\int K\,dW∫(aH+bK)dW=a∫HdW+b∫KdW.
  • Zero mean: E[It]=0\mathbb{E}[I_t] = 0E[It​]=0, since E[ξiΔWi]=E[ξi E[ΔWi∣Fti]]=0\mathbb{E}[\xi_i \Delta W_i] = \mathbb{E}\big[\xi_i\,\mathbb{E}[\Delta W_i \mid \mathcal{F}_{t_i}]\big] = 0E[ξi​ΔWi​]=E[ξi​E[ΔWi​∣Fti​​]]=0.
  • Martingale: E[It∣Fs]=Is\mathbb{E}[I_t \mid \mathcal{F}_s] = I_sE[It​∣Fs​]=Is​ for s≤ts \le ts≤t, and III has a version with continuous paths.
  • Isometry and covariance: E[It2]=E[∫0tHs2 ds]\mathbb{E}[I_t^2] = \mathbb{E}\big[\int_0^t H_s^2\,ds\big]E[It2​]=E[∫0t​Hs2​ds] and E[∫0tH dW∫0tK dW]=E[∫0tHsKs ds]\mathbb{E}\big[\int_0^t H\,dW \int_0^t K\,dW\big] = \mathbb{E}\big[\int_0^t H_s K_s\,ds\big]E[∫0t​HdW∫0t​KdW]=E[∫0t​Hs​Ks​ds].
  • Quadratic variation: [I]t=∫0tHs2 ds[I]_t = \int_0^t H_s^2\,ds[I]t​=∫0t​Hs2​ds, the integral analogue of [W]t=t[W]_t = t[W]t​=t.

For a deterministic integrand hhh the integral is Gaussian: ∫0th(s) dWs∼N(0,∫0th(s)2 ds)\int_0^t h(s)\,dW_s \sim N\big(0, \int_0^t h(s)^2\,ds\big)∫0t​h(s)dWs​∼N(0,∫0t​h(s)2ds). For example, ∫01s dWs∼N(0,1/3)\int_0^1 s\,dW_s \sim N(0, 1/3)∫01​sdWs​∼N(0,1/3).

Computing the integral of W against itself

Take H=WH = WH=W on [0,t][0, t][0,t] with the grid ti=it/nt_i = it/nti​=it/n. The left-point sum can be rewritten with the identity b(a−b)=12(a2−b2)−12(a−b)2b(a - b) = \tfrac12(a^2 - b^2) - \tfrac12(a - b)^2b(a−b)=21​(a2−b2)−21​(a−b)2, applied to a=Wti+1a = W_{t_{i+1}}a=Wti+1​​, b=Wtib = W_{t_i}b=Wti​​:

∑i=0n−1Wti ΔWi=12∑i=0n−1(Wti+12−Wti2)−12∑i=0n−1(ΔWi)2=12Wt2−12∑i=0n−1(ΔWi)2.\sum_{i=0}^{n-1} W_{t_i}\,\Delta W_i = \frac{1}{2}\sum_{i=0}^{n-1}\big(W_{t_{i+1}}^2 - W_{t_i}^2\big) - \frac{1}{2}\sum_{i=0}^{n-1}(\Delta W_i)^2 = \frac{1}{2}W_t^2 - \frac{1}{2}\sum_{i=0}^{n-1}(\Delta W_i)^2.i=0∑n−1​Wti​​ΔWi​=21​i=0∑n−1​(Wti+1​2​−Wti​2​)−21​i=0∑n−1​(ΔWi​)2=21​Wt2​−21​i=0∑n−1​(ΔWi​)2.

The first sum telescopes, and this is an exact identity for every nnn. The second sum is the only random approximation left. Each (ΔWi)2(\Delta W_i)^2(ΔWi​)2 has mean t/nt/nt/n and variance 2(t/n)22(t/n)^22(t/n)2, and they are independent, so

E[∑i(ΔWi)2]=t,Var⁡(∑i(ΔWi)2)=2t2n⟶0.\mathbb{E}\Big[\sum_i (\Delta W_i)^2\Big] = t, \qquad \operatorname{Var}\Big(\sum_i (\Delta W_i)^2\Big) = \frac{2t^2}{n} \longrightarrow 0.E[i∑​(ΔWi​)2]=t,Var(i∑​(ΔWi​)2)=n2t2​⟶0.

The sum converges to ttt in L2L^2L2, and therefore

∫0tWs dWs=Wt2−t2.\int_0^t W_s\,dW_s = \frac{W_t^2 - t}{2}.∫0t​Ws​dWs​=2Wt2​−t​.

The −t/2-t/2−t/2 is exactly half the quadratic variation. Ordinary calculus would give Wt2/2W_t^2/2Wt2​/2; the correction is there because squared increments of Brownian motion do not vanish in the limit. The result is consistent with the properties above: its mean is (t−t)/2=0(t - t)/2 = 0(t−t)/2=0, and its variance is Var⁡(Wt2)/4=2t2/4=t2/2\operatorname{Var}(W_t^2)/4 = 2t^2/4 = t^2/2Var(Wt2​)/4=2t2/4=t2/2, which matches the isometry ∫0tE[Ws2] ds=∫0ts ds=t2/2\int_0^t \mathbb{E}[W_s^2]\,ds = \int_0^t s\,ds = t^2/2∫0t​E[Ws2​]ds=∫0t​sds=t2/2.

Numbers: if W1=1.3W_1 = 1.3W1​=1.3, then ∫01W dW=(1.69−1)/2=0.345\int_0^1 W\,dW = (1.69 - 1)/2 = 0.345∫01​WdW=(1.69−1)/2=0.345, while W12/2=0.845W_1^2/2 = 0.845W12​/2=0.845. On a grid of nnn steps the left-point sum differs from (W12−1)/2(W_1^2 - 1)/2(W12​−1)/2 by −12(∑i(ΔWi)2−1)-\tfrac12\big(\sum_i(\Delta W_i)^2 - 1\big)−21​(∑i​(ΔWi​)2−1), whose standard deviation is 1/(2n)\sqrt{1/(2n)}1/(2n)​: 0.2240.2240.224 for n=10n = 10n=10, 0.07070.07070.0707 for n=100n = 100n=100 and 0.007070.007070.00707 for n=10 000n = 10\,000n=10000.

Left-point and midpoint evaluation

Evaluate instead at the midpoint mi=(ti+ti+1)/2m_i = (t_i + t_{i+1})/2mi​=(ti​+ti+1​)/2. Split ΔWi=Ai+Bi\Delta W_i = A_i + B_iΔWi​=Ai​+Bi​ with Ai=Wmi−WtiA_i = W_{m_i} - W_{t_i}Ai​=Wmi​​−Wti​​ and Bi=Wti+1−WmiB_i = W_{t_{i+1}} - W_{m_i}Bi​=Wti+1​​−Wmi​​, independent N(0,t/(2n))N(0, t/(2n))N(0,t/(2n)) variables. Since Wmi=Wti+AiW_{m_i} = W_{t_i} + A_iWmi​​=Wti​​+Ai​,

∑iWmi ΔWi=∑iWti ΔWi+∑iAi2+∑iAiBi⟶Wt2−t2+t2+0=Wt22.\sum_i W_{m_i}\,\Delta W_i = \sum_i W_{t_i}\,\Delta W_i + \sum_i A_i^2 + \sum_i A_i B_i \longrightarrow \frac{W_t^2 - t}{2} + \frac{t}{2} + 0 = \frac{W_t^2}{2}.i∑​Wmi​​ΔWi​=i∑​Wti​​ΔWi​+i∑​Ai2​+i∑​Ai​Bi​⟶2Wt2​−t​+2t​+0=2Wt2​​.

Here ∑iAi2→t/2\sum_i A_i^2 \to t/2∑i​Ai2​→t/2 by the same variance argument as before, and ∑iAiBi\sum_i A_i B_i∑i​Ai​Bi​ has mean zero and variance n (t/2n)2→0n\,(t/2n)^2 \to 0n(t/2n)2→0. The midpoint limit is the Stratonovich integral, written ∫0tWs∘dWs=Wt2/2\int_0^t W_s \circ dW_s = W_t^2/2∫0t​Ws​∘dWs​=Wt2​/2. In general ∫H∘dW=∫H dW+12[H,W]t\int H \circ dW = \int H\,dW + \tfrac12 [H, W]_t∫H∘dW=∫HdW+21​[H,W]t​, where [H,W][H, W][H,W] is the quadratic covariation; for H=WH = WH=W this is t/2t/2t/2.

The figure draws one Brownian path on [0,1][0, 1][0,1] and, for the chosen number of steps nnn (10, 100, 1000 or 10 000), plots against ttt the running left-point sum ∑WtiΔWi\sum W_{t_i} \Delta W_i∑Wti​​ΔWi​, the running midpoint sum ∑W(ti+ti+1)/2ΔWi\sum W_{(t_i + t_{i+1})/2} \Delta W_i∑W(ti​+ti+1​)/2​ΔWi​, and the curves (Wt2−t)/2(W_t^2 - t)/2(Wt2​−t)/2 and Wt2/2W_t^2/2Wt2​/2. Increase nnn and note that the left sum settles on (Wt2−t)/2(W_t^2 - t)/2(Wt2​−t)/2 and the midpoint sum on Wt2/2W_t^2/2Wt2​/2. The readouts at t=1t = 1t=1 compare the left sum with (W12−1)/2(W_1^2 - 1)/2(W12​−1)/2, the midpoint sum with W12/2W_1^2/2W12​/2, and show the difference of the two sums, which approaches 0.5=t/20.5 = t/20.5=t/2 as nnn grows. Press "New path" to repeat with another path; the gap of t/2t/2t/2 stays.

Left-point and midpoint sumsOpen in Itô Process →
(W² − t)/2Left-point sum (Itô)W²/2Midpoint sum (Stratonovich)
Left sum · (W₁² − 1)/2
-0.500 · -0.496
Midpoint sum · W₁²/2
-0.005 · 0.004
Difference (→ 0.5)
0.495

The two integrals have complementary strengths. Stratonovich integrals obey the ordinary chain rule, which is why they are common in physics, where white noise models a smooth but fast physical signal. Itô integrals are martingales: E[∫0tW dW]=0\mathbb{E}\big[\int_0^t W\,dW\big] = 0E[∫0t​WdW]=0, while E[∫0tW∘dW]=t/2\mathbb{E}\big[\int_0^t W \circ dW\big] = t/2E[∫0t​W∘dW]=t/2. Finance uses Itô for a modelling reason. A trader chooses the position ξi\xi_iξi​ at the start of the period, using the information Fti\mathcal{F}_{t_i}Fti​​, and then the price moves. The gain ξi (Sti+1−Sti)\xi_i\,(S_{t_{i+1}} - S_{t_i})ξi​(Sti+1​​−Sti​​) is a left-point sum. Evaluating at the midpoint would let the position depend on half of the price move it is paid for, which is a form of trading on future information. Under a measure where discounted prices are martingales, the Itô integral keeps the gains from any such strategy a martingale too, which is the continuous-time statement that no strategy turns a fair game into a sure profit.

Summary

  • Brownian paths have infinite variation, so ∫H dW\int H\,dW∫HdW cannot be defined path by path; the value of the sums depends on the evaluation point.
  • For simple adapted HHH the Itô integral is ∑iξi ΔWi\sum_i \xi_i\,\Delta W_i∑i​ξi​ΔWi​ with ξi\xi_iξi​ known at tit_iti​; the Itô isometry E[(∫H dW)2]=E[∫H2 ds]\mathbb{E}[(\int H\,dW)^2] = \mathbb{E}[\int H^2\,ds]E[(∫HdW)2]=E[∫H2ds] extends it to all of LT2\mathcal{L}^2_TLT2​ by L2L^2L2 limits.
  • The Itô integral is linear, has mean zero, is a martingale and has quadratic variation ∫H2 ds\int H^2\,ds∫H2ds.
  • ∫0tW dW=(Wt2−t)/2\int_0^t W\,dW = (W_t^2 - t)/2∫0t​WdW=(Wt2​−t)/2: the −t/2-t/2−t/2 comes from ∑(ΔWi)2→t\sum (\Delta W_i)^2 \to t∑(ΔWi​)2→t. Midpoint evaluation gives the Stratonovich value Wt2/2W_t^2/2Wt2​/2.
  • Finance uses Itô because positions are chosen before the price moves they are paid on.

The next lesson, "Itô's lemma", turns the computation for W2W^2W2 into a general rule: d(Wt2)=2Wt dWt+dtd(W_t^2) = 2W_t\,dW_t + dtd(Wt2​)=2Wt​dWt​+dt is the differential form of ∫0tW dW=(Wt2−t)/2\int_0^t W\,dW = (W_t^2 - t)/2∫0t​WdW=(Wt2​−t)/2.

Try itSee the Itô correction at work on simulated log prices in the Itô Process lab →
3 questions on this lesson at the end of the series →
← PreviousThe Brownian bridgeNext →Itô's lemma
Connections
The Itô integral
Open the full graph →
Labs for this lesson
Itô ProcessBrownian Motion