Derivation of the Black-Scholes-Merton PDE › Worked Example: Verifying the PDE

Michael BrenndoerferNovember 26, 20254 min read
Quantitative FinanceAll chapters
Part 5 of 15Stochastic Calculus and Pricing Theory
Chapter 29Derivation of the Black-Scholes-Merton PDE › Worked Example: Verifying the PDE
Reading section 8 of 16 · Worked Example: Verifying the PDE · 4 minSection 8 of 16 · 4 min

Book contents

Section 8 of 16 · 4 min read

Derivation of the Black-Scholes-Merton PDE › Worked Example: Verifying the PDEOne section · Sections in this chapter8 / 16 ↓
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 verify that a proposed solution satisfies the Black-Scholes PDE. This exercise builds intuition for what it means to "satisfy" a partial differential equation and demonstrates how to check whether a candidate pricing formula is valid. Consider a simple derivative whose price is V(S,t)=SV(S, t) = S, representing a forward contract with zero delivery price and zero interest rates.

We need to check whether V=SV = S satisfies:

∂V∂t+rS∂V∂S+12σ2S2∂2V∂S2=rV\frac{\partial V}{\partial t} + rS\frac{\partial V}{\partial S} + \frac{1}{2} \sigma^2 S^2 \frac{\partial^2 V}{\partial S^2} = rV

where:

  • VV: price of the derivative
  • SS: stock price
  • rr: risk-free interest rate
  • σ\sigma: volatility
  • tt: time

Computing the partial derivatives is straightforward since V=SV = S is a simple function:

∂V∂t=0,∂V∂S=1,∂2V∂S2=0\frac{\partial V}{\partial t} = 0, \quad \frac{\partial V}{\partial S} = 1, \quad \frac{\partial^2 V}{\partial S^2} = 0

where:

  • VV: derivative price (SS)

The time derivative is zero because V=SV = S has no explicit dependence on time. The first derivative with respect to SS is one because VV increases linearly with SS. The second derivative is zero because the relationship is perfectly linear with no curvature.

Substituting into the PDE:

0+rS(1)+12σ2S2(0)=rS0 + rS(1) + \frac{1}{2}\sigma^2 S^2(0) = rS

where:

  • r,S,σr, S, \sigma: model parameters

The left side equals rSrS, and the right side is rV=rSrV = rS. The equation is satisfied, confirming that the stock price itself is a trivial solution to the BSM PDE. This makes economic sense: holding the stock is a valid "derivative strategy" that trivially satisfies no-arbitrage pricing.

Now consider a more interesting example: a derivative that pays S2S^2 at maturity TT. Such a derivative might seem exotic, but it illustrates how the PDE handles nonlinear payoffs. We'll verify that a solution of the following form satisfies the PDE:

V(S,t)=S2e(r+σ2)(T−t)V(S, t) = S^2 e^{(r + \sigma^2)(T-t)}

where:

  • V(S,t)V(S, t): price of the derivative
  • SS: stock price
  • rr: risk-free interest rate
  • σ\sigma: volatility
  • TT: maturity time
  • tt: current time

Notice that this formula has an interesting structure: it equals S2S^2 at maturity (when T−t=0T - t = 0), satisfying the terminal condition, and it includes an exponential factor that accounts for the time value and volatility effects.

Computing the partial derivatives requires careful application of the chain rule:

∂V∂t=−S2(r+σ2)e(r+σ2)(T−t)∂V∂S=2Se(r+σ2)(T−t)∂2V∂S2=2e(r+σ2)(T−t)\begin{aligned} \frac{\partial V}{\partial t} &= -S^2(r + \sigma^2)e^{(r + \sigma^2)(T-t)} \\ \frac{\partial V}{\partial S} &= 2Se^{(r + \sigma^2)(T-t)} \\ \frac{\partial^2 V}{\partial S^2} &= 2e^{(r + \sigma^2)(T-t)} \end{aligned}

where:

  • VV: derivative price (S2e(r+σ2)(T−t)S^2 e^{(r + \sigma^2)(T-t)})
  • S,t,r,σ,TS, t, r, \sigma, T: model parameters

The time derivative is negative because the exponential factor decreases as we approach maturity (when T−tT - t shrinks). The first spatial derivative reflects that the price increases with SS, at a rate proportional to SS itself. The second spatial derivative captures the constant curvature of the quadratic relationship.

Substituting these derivatives into the left-hand side (LHS) of the PDE and simplifying:

LHS=−S2(r+σ2)e(r+σ2)(T−t)+rS⋅2Se(r+σ2)(T−t)+12σ2S2⋅2e(r+σ2)(T−t)=S2e(r+σ2)(T−t)[−(r+σ2)+2r+σ2](factor out common term)=S2e(r+σ2)(T−t)[−r−σ2+2r+σ2](expand terms)=S2e(r+σ2)(T−t)⋅r(simplify)=rV(definition of V)\begin{aligned} \text{LHS} &= -S^2(r + \sigma^2)e^{(r+\sigma^2)(T-t)} + rS \cdot 2Se^{(r+\sigma^2)(T-t)} + \frac{1}{2}\sigma^2 S^2 \cdot 2e^{(r+\sigma^2)(T-t)} \\ &= S^2 e^{(r+\sigma^2)(T-t)} \left[ -(r + \sigma^2) + 2r + \sigma^2 \right] && \text{(factor out common term)} \\ &= S^2 e^{(r+\sigma^2)(T-t)} \left[ -r - \sigma^2 + 2r + \sigma^2 \right] && \text{(expand terms)} \\ &= S^2 e^{(r+\sigma^2)(T-t)} \cdot r && \text{(simplify)} \\ &= rV && \text{(definition of V)} \end{aligned}

where:

  • LHS: value of left-hand side of PDE
  • VV: derivative price

The algebraic simplification is satisfying: the terms involving σ2\sigma^2 cancel (one from the time derivative and one from the gamma term), and the terms involving rr combine to give exactly what we need.

The right side of the PDE is:

rV=r⋅S2e(r+σ2)(T−t)rV = r \cdot S^2 e^{(r+\sigma^2)(T-t)}

where:

  • rr: risk-free rate
  • VV: derivative price

Both sides are equal, confirming our solution. This example illustrates that derivatives with nonlinear payoffs require appropriate adjustments to their prices, captured by factors that depend on volatility as well as interest rates.