Skip to content
homeaboutworkprojectsthesiswritingresume
Loading
~/blog/stochastic-differential-equations0%dark
  1. home/
  2. writing/
  3. Stochastic Differential Equations

04 June 2026 · 7 min read · updated 13 June 2026

Stochastic Differential Equations

A stochastic differential equation specifies a process through its drift and diffusion, and under a Lipschitz condition it has a unique strong solution. We define the strong solution, prove Gronwall's inequality, prove uniqueness by bounding the mean-square distance between two solutions and closing it with Gronwall, and prove existence by Picard iteration, where the Ito isometry and Doob's inequality produce a factorial decay that makes the successive approximations converge. The Ornstein-Uhlenbeck process and the diffusion models of finance are the solutions this theorem guarantees.

  • 7 equations
  • 7 results
  • 10 connections
  • stochastic-processes
  • stochastic-calculus
  • sde
On this page▾
  • Strong solutions
  • Gronwall's inequality
  • Uniqueness
  • Existence

7 min left

  • Strong solutions1m
  • Gronwall's inequality1m
  • Uniqueness1m
  • Existence3m

A stochastic differential equation describes a process by how it moves in the next instant, a deterministic drift plus a random push proportional to a Brownian increment. The Ito integral gives the equation a rigorous meaning, and the question is whether the description determines a process. Under a Lipschitz condition it does, uniquely, and the proof is the stochastic version of the Picard-Lindelof argument for ordinary differential equations, with the Ito isometry controlling the random part. This post proves existence and uniqueness [1], [2]. Here WWW is a Brownian motion and the coefficients b,σ:[0,T]×R→Rb,\sigma:[0,T]\times\R\to\Rb,σ:[0,T]×R→R are measurable.

#Strong solutions

Definition1

A strong solution of the stochastic differential equation

dXt=b(t,Xt) dt+σ(t,Xt) dWt,X0=x0,(1)dX_t=b(t,X_t)\,dt+\sigma(t,X_t)\,dW_t,\qquad X_0=x_0, \tag{1}dXt​=b(t,Xt​)dt+σ(t,Xt​)dWt​,X0​=x0​,(1)

is an adapted continuous process XXX with ∫0TE[b(t,Xt)2+σ(t,Xt)2] dt<∞\int_0^T\E[b(t,X_t)^2+\sigma(t,X_t)^2]\,dt<\infty∫0T​E[b(t,Xt​)2+σ(t,Xt​)2]dt<∞ satisfying the integral equation Xt=x0+∫0tb(s,Xs) ds+∫0tσ(s,Xs) dWsX_t=x_0+\int_0^t b(s,X_s)\,ds+\int_0^t\sigma(s,X_s)\,dW_sXt​=x0​+∫0t​b(s,Xs​)ds+∫0t​σ(s,Xs​)dWs​ for all t≤Tt\le Tt≤T almost surely.

The equation is the integral identity, the differential being shorthand. The coefficients are assumed Lipschitz and of linear growth, meaning there is a constant KKK with

∣b(t,x)−b(t,y)∣+∣σ(t,x)−σ(t,y)∣≤K∣x−y∣,b(t,x)2+σ(t,x)2≤K2(1+x2),(2)\abs{b(t,x)-b(t,y)}+\abs{\sigma(t,x)-\sigma(t,y)}\le K\abs{x-y},\qquad b(t,x)^2+\sigma(t,x)^2\le K^2(1+x^2), \tag{2}∣b(t,x)−b(t,y)∣+∣σ(t,x)−σ(t,y)∣≤K∣x−y∣,b(t,x)2+σ(t,x)2≤K2(1+x2),(2)

for all t,x,yt,x,yt,x,y. Uniqueness uses the Lipschitz bound; existence additionally uses the linear-growth bound.

#Gronwall's inequality

The analytic engine is a comparison principle that turns a self-referential integral bound into an exponential one.

Lemma2

If a nonnegative integrable ggg satisfies g(t)≤a+C∫0tg(s) dsg(t)\le a+C\int_0^t g(s)\,dsg(t)≤a+C∫0t​g(s)ds for all t≤Tt\le Tt≤T with constants a,C≥0a,C\ge 0a,C≥0, then g(t)≤a eCtg(t)\le a\,e^{Ct}g(t)≤aeCt.

Proof

If C=0C=0C=0 the hypothesis is g(t)≤a=a e0⋅tg(t)\le a=a\,e^{0\cdot t}g(t)≤a=ae0⋅t and there is nothing to prove, so take C>0C>0C>0. Let G(t)=∫0tg(s) dsG(t)=\int_0^t g(s)\,dsG(t)=∫0t​g(s)ds, which is absolutely continuous with G′(t)=g(t)≤a+CG(t)G'(t)=g(t)\le a+CG(t)G′(t)=g(t)≤a+CG(t) for almost every ttt, so ddt(e−CtG(t))=e−Ct(G′−CG)≤ae−Ct\frac{d}{dt}\big(e^{-Ct}G(t)\big)=e^{-Ct}(G'-CG)\le ae^{-Ct}dtd​(e−CtG(t))=e−Ct(G′−CG)≤ae−Ct almost everywhere. Integrating from 000 to ttt with G(0)=0G(0)=0G(0)=0 gives e−CtG(t)≤aC(1−e−Ct)e^{-Ct}G(t)\le\frac aC(1-e^{-Ct})e−CtG(t)≤Ca​(1−e−Ct), that is CG(t)≤a(eCt−1)CG(t)\le a(e^{Ct}-1)CG(t)≤a(eCt−1). Substituting back, g(t)≤a+CG(t)≤a+a(eCt−1)=a eCtg(t)\le a+CG(t)\le a+a(e^{Ct}-1)= a\,e^{Ct}g(t)≤a+CG(t)≤a+a(eCt−1)=aeCt.

#Uniqueness

Two solutions of the same equation cannot separate, because the Lipschitz condition makes their distance control its own growth.

Theorem3

The equation Equation (1) has at most one strong solution, up to indistinguishability.

Proof

Let XXX and YYY be strong solutions with the same initial value. Their difference is

Xt−Yt=∫0t(b(s,Xs)−b(s,Ys)) ds+∫0t(σ(s,Xs)−σ(s,Ys)) dWs,(3)X_t-Y_t=\int_0^t\big(b(s,X_s)-b(s,Y_s)\big)\,ds+\int_0^t\big(\sigma(s,X_s)-\sigma(s,Y_s)\big)\,dW_s, \tag{3}Xt​−Yt​=∫0t​(b(s,Xs​)−b(s,Ys​))ds+∫0t​(σ(s,Xs​)−σ(s,Ys​))dWs​,(3)

the initial values cancelling. Squaring and using (u+v)2≤2u2+2v2(u+v)^2\le 2u^2+2v^2(u+v)2≤2u2+2v2, then taking expectations, the Ito isometry turns the stochastic integral's second moment into an ordinary integral and the Cauchy-Schwarz inequality bounds the drift integral,

E[(Xt−Yt)2]≤2t∫0tE[(b(s,Xs)−b(s,Ys))2] ds+2∫0tE[(σ(s,Xs)−σ(s,Ys))2] ds.(4)\E\big[(X_t-Y_t)^2\big]\le 2t\int_0^t\E\big[(b(s,X_s)-b(s,Y_s))^2\big]\,ds+2\int_0^t\E\big[(\sigma(s,X_s)- \sigma(s,Y_s))^2\big]\,ds. \tag{4}E[(Xt​−Yt​)2]≤2t∫0t​E[(b(s,Xs​)−b(s,Ys​))2]ds+2∫0t​E[(σ(s,Xs​)−σ(s,Ys​))2]ds.(4)

The Lipschitz condition Equation (2) bounds each integrand by K2E[(Xs−Ys)2]K^2\E[(X_s-Y_s)^2]K2E[(Xs​−Ys​)2], so with C=2(T+1)K2C=2(T+1)K^2C=2(T+1)K2 the function g(t)=E[(Xt−Yt)2]g(t)=\E[(X_t-Y_t)^2]g(t)=E[(Xt​−Yt​)2] satisfies g(t)≤C∫0tg(s) dsg(t)\le C\int_0^t g(s)\,dsg(t)≤C∫0t​g(s)ds, the Gronwall hypothesis with a=0a=0a=0. Here ggg is finite and integrable on [0,T][0,T][0,T], since the strong-solution definition gives E[Xt2]≤3x02+3T∫0tE[b2] ds+3∫0tE[σ2] ds<∞\E[X_t^2]\le 3x_0^2+3T\int_0^t\E[b^2]\,ds+3\int_0^t\E[\sigma^2]\,ds<\inftyE[Xt2​]≤3x02​+3T∫0t​E[b2]ds+3∫0t​E[σ2]ds<∞ by Cauchy-Schwarz and the Ito isometry, and likewise for YYY, so g(t)≤2E[Xt2]+2E[Yt2]g(t)\le 2\E[X_t^2]+2\E[Y_t^2]g(t)≤2E[Xt2​]+2E[Yt2​] is bounded. By Lemma 2, g≡0g\equiv 0g≡0, so Xt=YtX_t=Y_tXt​=Yt​ almost surely for each ttt, and both being continuous they are indistinguishable.

#Existence

The solution is built by iterating the equation, starting from the constant path and feeding each approximation back through the drift and diffusion.

Theorem4

Under Equation (2) the equation Equation (1) has a strong solution on [0,T][0,T][0,T].

Proof

Define the Picard iterates Xt0=x0X^0_t=x_0Xt0​=x0​ and

Xtn+1=x0+∫0tb(s,Xsn) ds+∫0tσ(s,Xsn) dWs.(5)X^{n+1}_t=x_0+\int_0^t b(s,X^n_s)\,ds+\int_0^t\sigma(s,X^n_s)\,dW_s. \tag{5}Xtn+1​=x0​+∫0t​b(s,Xsn​)ds+∫0t​σ(s,Xsn​)dWs​.(5)

Each iterate is adapted and continuous, and a finite second moment propagates along the recursion. Let An=E[sup⁡s≤T(Xsn)2]A_n=\E[\sup_{s\le T}(X^n_s)^2]An​=E[sups≤T​(Xsn​)2]. Then A0=x02<∞A_0=x_0^2<\inftyA0​=x02​<∞, and if An<∞A_n<\inftyAn​<∞ then linear growth gives ∫0TE[σ(s,Xsn)2] ds≤K2T(1+An)<∞\int_0^T\E[\sigma(s,X^n_s)^2]\,ds\le K^2T(1+A_n)<\infty∫0T​E[σ(s,Xsn​)2]ds≤K2T(1+An​)<∞, so the Ito integral in Equation (5) is a genuine L2L^2L2 martingale and both the Ito isometry and Doob's L2L^2L2 inequality apply; bounding Xn+1X^{n+1}Xn+1 by (u+v+w)2≤3u2+3v2+3w2(u+v+w)^2\le 3u^2+3v^2+3w^2(u+v+w)2≤3u2+3v2+3w2, Cauchy-Schwarz on the drift, and Doob on the martingale, An+1≤3x02+3TK2T(1+An)+12K2T(1+An)<∞A_{n+1}\le 3x_0^2+3TK^2T(1+A_n)+12K^2T(1+A_n)<\inftyAn+1​≤3x02​+3TK2T(1+An​)+12K2T(1+An​)<∞. By induction every AnA_nAn​ is finite, so each integral is square-integrable and the recursion is well defined. Write Δtn=Xtn+1−Xtn\Delta^n_t=X^{n+1}_t-X^n_tΔtn​=Xtn+1​−Xtn​ and φn(t)=E[sup⁡s≤t(Δsn)2]\varphi_n(t)=\E\big[\sup_{s\le t}(\Delta ^n_s)^2\big]φn​(t)=E[sups≤t​(Δsn​)2], which the bound An<∞A_n<\inftyAn​<∞ renders finite. For n≥1n\ge 1n≥1, the increment Δn\Delta^nΔn is Equation (3) with Xn,Xn−1X^n,X^{n-1}Xn,Xn−1 in place of X,YX,YX,Y, so by (u+v)2≤2u2+2v2(u+v)^2\le 2u^2+2v^2(u+v)2≤2u2+2v2, the Cauchy-Schwarz bound on the drift, and Doob's L2L^2L2 inequality with the Ito isometry on the martingale part,

φn(t)≤2T∫0tE[(b(s,Xn)−b(s,Xn−1))2] ds+8∫0tE[(σ(s,Xn)−σ(s,Xn−1))2] ds,(6)\varphi_n(t)\le 2T\int_0^t\E\big[(b(s,X^n)-b(s,X^{n-1}))^2\big]\,ds+8\int_0^t\E\big[(\sigma(s,X^n)-\sigma(s, X^{n-1}))^2\big]\,ds, \tag{6}φn​(t)≤2T∫0t​E[(b(s,Xn)−b(s,Xn−1))2]ds+8∫0t​E[(σ(s,Xn)−σ(s,Xn−1))2]ds,(6)

the factor 888 being twice Doob's constant 444. The Lipschitz bound makes both integrands at most K2E[(Δsn−1)2]≤K2φn−1(s)K^2\E[(\Delta^{n-1}_s)^2]\le K^2\varphi_{n-1}(s)K2E[(Δsn−1​)2]≤K2φn−1​(s), so with C=2(T+4)K2C=2(T+4)K^2C=2(T+4)K2,

φn(t)≤C∫0tφn−1(s) ds.(7)\varphi_n(t)\le C\int_0^t\varphi_{n-1}(s)\,ds. \tag{7}φn​(t)≤C∫0t​φn−1​(s)ds.(7)

Since φ0(T)=E[sup⁡s≤T(Xs1−x0)2]\varphi_0(T)=\E[\sup_{s\le T}(X^1_s-x_0)^2]φ0​(T)=E[sups≤T​(Xs1​−x0​)2] is a finite constant MMM by linear growth, iterating Equation (7) gives φn(t)≤M (Ct)n/n!\varphi_n(t)\le M\,(Ct)^n/n!φn​(t)≤M(Ct)n/n!, whose square roots sum to a finite value, ∑nφn(T)1/2≤∑nM(CT)n/n!<∞\sum_n\varphi_n(T)^{1 /2}\le\sum_n\sqrt{M(CT)^n/n!}<\infty∑n​φn​(T)1/2≤∑n​M(CT)n/n!​<∞, making the iterates a Cauchy sequence in the complete space of adapted processes with finite E[sup⁡s≤T(⋅)2]\E[\sup_{s\le T}(\cdot)^2]E[sups≤T​(⋅)2], and converge uniformly in mean square to a limit XXX, which is adapted and has a continuous version, since a subsequence converges uniformly in ttt almost surely.

Passing to the limit in Equation (5) identifies XXX as a solution. The drift integral converges because ∫0t(b(s,Xn)−b(s,X)) ds\int_0^t(b(s,X^n)-b(s,X))\,ds∫0t​(b(s,Xn)−b(s,X))ds has mean square at most TK2∫0tE[(Xsn−Xs)2] ds→0T K^2\int_0^t\E[(X^n_s-X_s)^2]\,ds\to 0TK2∫0t​E[(Xsn​−Xs​)2]ds→0 by Lipschitz continuity, and the stochastic integral converges by the Ito isometry and the same Lipschitz bound. The mean-square convergence gives E[sup⁡s≤TXs2]<∞\E[\sup_{s\le T}X_s^2]<\inftyE[sups≤T​Xs2​]<∞, so linear growth yields ∫0TE[σ(t,Xt)2] dt≤K2T(1+E[sup⁡X2])<∞\int_0^T\E[\sigma(t,X_t)^2]\,dt\le K^2T(1+\E[\sup X^2])<\infty∫0T​E[σ(t,Xt​)2]dt≤K2T(1+E[supX2])<∞ and likewise for bbb, the integrability required by Definition 1. So Xt=x0+∫0tb(s,Xs) ds+∫0tσ(s,Xs) dWsX_t=x_0+\int_0^t b(s,X_s)\,ds+\int_0^t\sigma(s,X_s)\,dW_sXt​=x0​+∫0t​b(s,Xs​)ds+∫0t​σ(s,Xs​)dWs​, a strong solution.

The factorial (Ct)n/n!(Ct)^n/n!(Ct)n/n! is the stochastic echo of the same factor in the Picard proof for ordinary differential equations, the contraction that beats the linear accumulation of error. Together the two theorems make a stochastic differential equation a genuine definition of a process, licensing a model written from its drift and diffusion and treated as a well-defined object. The Ornstein-Uhlenbeck process is the solution of the linear equation with mean reversion, the geometric Brownian motion verified by Ito's formula is the solution of the proportional-growth equation, and the diffusion models of derivatives pricing and optimal execution are instances covered by this theorem.

[1]
B. Øksendal, Stochastic Differential Equations: An Introduction with Applications, 6th ed. Springer, 2003.
[2]
I. Karatzas and S. E. Shreve, Brownian Motion and Stochastic Calculus, 2nd ed. Springer, 1991.

Part 7 of 8 in Stochastic Calculus

← previousChange of Measure and Girsanov's Theoremnext →The Ornstein-Uhlenbeck Process

Explore connections

see in the atlas →

related

  • Ito's Formula
  • Predictable Processes and Stopping Times
  • Martingales

referenced by (1)

  • Ito's Formula
cite
@misc{stochastic-differential-equations,
  author = {Zac Kienzle},
  title  = {Stochastic Differential Equations},
  year   = {2026},
  month  = {06},
  url    = {https://zackienzle.com/blog/stochastic-differential-equations}
}