Maths BSM Derivation
Maths BSM Derivation
Ser-Huang Poon
1 Black-Scholes Model
The Black-Scholes formula below is for pricing European call and put op-
tions:
where c (p) is the price of the European call (put), S0 is the current price of
the underlying assset, K is the strike or exercise price, r is the continuously
compounded risk free interest rate, T is the time to option maturity. N (d1 )
is the cumulativ probability distribution of a standard Normal distribution
for area below d1 , and N (−d1 ) = 1 − N (d1 ).
As T → 0,
d1 and d2 → ∞
N (d1 ) and N (d2 ) → 1
N (−d1 ) and N (−d2 ) → 0
1
which means
As σ → 0, again
The boundary conditions (3), (4), (5) and (6) are the boundary conditions
for checking option prices before using them for empirical tests. These
conditions are not specific to Black-Scholes. Option with market prices
(transaction or quote) violating these boundary conditions should be dis-
carded.
• The stock price follows a Geometric Brownian Motion with drift μ and
volatility σ
dS
= μ dt + σ dW (7)
S
where W is a Wiener process
• The underlying stock does not pay any dividend (cf. this assumption
can easily be extended to include dividend payments)
2
• Continuous trading (so that rebalancing of portfolio is done instanta-
neously)
which means that stock price has a lognormal distribution or the logarithm
of stock price has a normal distribution. In discrete time
µ ¶
1 2
d ln S = μ − σ dt + σdz
2
µ ¶
1 √
∆ ln S = μ − σ 2 ∆t + σ ε ∆t
2
∙µ ¶ ¸
1 2 √
ln ST − ln S0 ∼ N μ − σ T, σ T
2
∙ µ ¶ ¸
1 2 √
ln ST ∼ N ln S0 + μ − σ T, σ T . (11)
2
3
2.2 Black-Scholes PDE
The derivation of the Black-Scholes partial differential equation is based
on the fundamental fact that the option price and the stock price depend
on the same underlying source of uncertainty. A portfolio can then be
created consisting of the stock and the option which eliminates this source
of uncertainty. Given that this portfolio is riskless and must therefore earn
the risk-free rate of return. Here is how the logic works:
rΠ∆t = ∆Π
∂f
rΠ∆t = −∆f + ∆S
µ ¶ ∙ ∂S ¸
∂f ∂f ∂f 1 ∂2f 2 2 ∂f
r −f + S ∆t = − μS + + 2
σ S ∆t − σS∆z
∂S ∂S ∂t 2 ∂S ∂S
∂f
+ [μS∆t + σS∆z]
∂S
∂f ∂f ∂f 1 ∂2f 2 2
r (−f ) ∆t = −rS ∆t − μS∆t − ∆t − σ S ∆t
∂S ∂S ∂t 2 ∂S 2
∂f ∂f ∂f
− σS∆z + μS∆t + σS∆z
∂S ∂S ∂S
4
and finally we get the well known Black-Scholes PDE
∂f ∂f 1 ∂2f 2 2
rf = rS + + σ S (14)
∂S ∂t 2 ∂S 2
Now we are ready to solve the BS formular. First, the terminal value
of a call is
cT = E [max (S − K, 0)]
Z ∞
= (S − K) f (S) dS
K
Z ∞ Z ∞
ln S
= e f (ln S) d ln S − K f (ln S) d ln S.
ln K ln K
5
Substituting facts (ii) and (iii) and using information from (i) to set
1
μy = ln S0 + r − σ 2
2
σy = σ
a = ln K
we get
à ! à !
ln S0 + r + 12 σ 2 − ln K ln S0 + r − 12 σ 2 − ln K
cT = S0 er N − KN
ln K ln K
= S0 er N (d1 ) − KN (d2 ) (15)
where
ln SK0 + r − 12 σ 2
d1 =
σ
d2 = d1 − σ.
The present value of the call option is derived by applying e−r to both
sides. The put option price can be derived using put-call parity or use the
same argument above. The σ in the above formular is volatility over the
option maturity. If we use σ as the annualised volatility then we replace σ
√
with σ T in the formula.
There are important insights from (15), all valid only in a “risk neutral”
world:
(ii) Alternatively, N (d2 ) is the probability that call finishes in the money.
(iv) S0 erT N (d1 ) is the expected value E [ST − X]+ , where E [·]+ is expec-
tation computed for positive values only.