Chapter V: Mathematical Finance in Discrete Time: Finite
Chapter V: Mathematical Finance in Discrete Time: Finite
We follow [BK], Ch. 4 (or see the other sources cited in Ch0).
Here Hni denotes the number of shares of asset i held in the portfolio at
time n – to be determined on the basis of information available before time
n; the vector Hn = (Hn0 , Hn1 , · · · , Hnd ) is the portfolio at time n. Writing
Sn = (Sn0 , Sn1 , · · · , Snd ) for the vector of prices at time n, the value of the
portfolio at time n is the scalar product
1
The discounted value is
where βn := 1/Sn0 and S̃n = (1, βn Sn1 , · · · , βn Snd ) is the vector of discounted
prices.
Note. The previsibility of H reflects that there is no insider trading.
Interpretation. When new prices Sn are quoted at time n, the investor adjusts
his portfolio from Hn to Hn+1 , without bringing in or consuming any wealth.
For a SF strategy, the second term on the right is zero. Then the LHS, the
net increase in the value of the portfolio, is shown as due only to the price
changes Sn+1 − Sn . So for H ∈ SF ,
2
Proof. If H is SF, then as above
Hn0 = V0 + Σn1 (Hj1 ∆S̃j1 + · · · + Hjd ∆S̃jd ) − (Hn1 S̃n1 + · · · + Hnd S̃nd ),
which defines Hn0 uniquely. The terms in S̃ni are Hni ∆S̃ni − Hni S̃ni = −Hni S̃n−1
i
,
which is Fn−1 -measurable. So
Hn0 = V0 + Σn−1
1 (Hj1 ∆S̃j1 + · · · + Hjd ∆S̃jd ) − (Hn1 Sn−1
1
+ · · · + Hnd S̃n−1
d
),
where as H 1 , · · · , H d are predictable, all terms on the RHS are Fn−1 -measurable,
so H 0 is predictable. //
Numéraire. What units do we reckon value in? All that is really necessary is
that our chosen unit of account should always be positive (as we then reckon
our holdings by dividing by it, and one cannot divide by zero). Common
choices are pounds sterling (UK), dollars (US), euros etc. Gold is also pos-
sible (now priced in sterling etc. – but the pound sterling represented an
amount of gold, till the UK ‘went off the gold standard’). By contrast, risky
stocks can have value 0 (if the company goes bankrupt). We call such an
always-positive asset, used to reckon values in, a numéraire.
Of course, one has to be able to change numéraire – e.g. when going
from UK to the US or eurozone. As one would expect, this changes nothing
important. In particular, we quote (numéraire invariance theorem – see e.g.
[BK] Prop. 4.1.1) that the set SF of self-financing strategies is invariant un-
der change of numéraire.
Note. 1. This alerts us to what is meant by ‘risky’. To the owner of a gold-
mine, sterling is risky. The danger is not that the UK government might go
bankrupt, but that sterling might depreciate against the dollar, or euro, etc.
2. With this understood, we shall feel free to refer to our numéraire as ‘bank
account’. The point is that we don’t trade in it (why would a goldmine owner
trade in gold?); it is the other – ‘risky’ – assets that we trade in.
3
§2. Viability (NA): Existence of Equivalent Martingale Measures.
Although we are allowed to borrow (from the bank), and sell (stocks)
short, we are – naturally – required to stay solvent (recall that trading while
insolvent is an offence under the Companies Act!).
If the strategy is admissible and its initial value – the RHS above – is zero,
the LHS E ∗ [ṼN (H)] is zero, but ṼN (H) ≥ 0 (by admissibility). Since each
P ({ω}) > 0 (by assumption), each P ∗ ({ω}) > 0 (by equivalence). This and
ṼN (H) ≥ 0 force ṼN (H) = 0 (sum of non-negatives can only be 0 if each
4
term is 0). So no arbitrage is possible. //
The converse is true, but harder, and needs a preparatory result – which
is interesting and important in its own right.
Separating Hyperplane Theorem (SHT).
In a vector space V , a hyperplane is a translate of a (vector) subspace
U of codimension 1 – that is, U and some one-dimensional subspace, say R,
together span V : V is the direct sum V = U ⊕ R (e.g., R3 = R2 ⊕ R). Then
H = [f, α] := {x : f (x) = α}
for some α and linear functional f . In the finite-dimensional case, of dimen-
sion n, say, one can think of f (x) as an inner product,
f (x) = f.x = f1 x1 + . . . + fn xn .
The hyperplane H = [f, α] separates sets A, B ⊂ V if
f (x) ≥ α ∀ x ∈ A, f (x) ≤ α ∀x∈B
(or the same inequalities with A, B, or ≥, ≤, interchanged).
Call a set A in a vector space V convex if
x, y ∈ A, 0≤λ≤1 ⇒ λx + (1 − λ)y ∈ A
– that is, A contains the line-segment joining any pair of its points.
We can now state (without proof) the SHT (see e,g, [BK] App. C).
SHT. Any two non-empty disjoint convex sets in a vector space can be sep-
arated by a hyperplane.
A cone is a subset of a vector space closed under vector addition and
multiplication by positive constants (so: like a vector subspace, but with a
sign-restriction in scalar multiplication).
We turn now to the proof of the converse.
Proof of the converse (not examinable). ⇒: Write Γ for the cone of strictly
positive random variables. Viability (NA) says that for any admissible strat-
egy H,
V0 (H) = 0 ⇒ ṼN (H) ∈
/ Γ. (∗)
To any admissible process (Hn1 , · · · , Hnd ), we associate its discounted cu-
mulative gain process
G̃n (H) := Σn1 (Hj1 ∆S̃j1 + · · · + Hjd ∆S̃jd ).
5
By the Proposition, we can extend (H1 , · · · , Hd ) to a unique predictable pro-
cess (Hn0 ) such that the strategy H = ((Hn0 , Hn1 , · · · , Hnd )) is self-financing
with initial value zero. By NA, G̃N (H) = 0 – that is, G̃N (H) ∈ / Γ.
We now form the set V of random variables G̃N (H), with H = (H 1 , · · · , H d )
a previsible process. This is a vector subspace of the vector space RΩ of ran-
dom variables on Ω, by linearity of the gain process G(H) in H. By (∗), this
subspace V does not meet Γ. So V does not meet the subset
K := {X ∈ Γ : Σω X(ω) = 1}.
Now K is a convex set not meeting the origin. By the Separating Hyper-
plane Theorem, there is a vector λ = (λ(ω) : ω ∈ Ω) such that for all X ∈ K
E ∗ [ΣN i i
1 Hj ∆S̃j ] = 0 (i = 1, · · · , d).
By the Martingale Transform Lemma, this says that the discounted price
processes (S̃ni ) are P ∗ -martingales. //
6
(so that we know how to evaluate h at the terminal time N ).
7
actually worth something. These sound similar; the Lemma shows they are
the same here. So we only need one term; we use SF as it is shorter, but
w.l.o.g. this means admissible also.]
Lemma. In a viable market, any attainable h (i.e., any h that can be repli-
cated by a SF strategy H) can also be replicated by an admissible strategy H.
U0 + ΣN
1 Hn .∆S̃n
8
Write kXk∞ := max{|X(ω)| : ω ∈ Ω}, and define P ∗∗ by
X(ω) ∗
P ∗∗ ({ω}) = 1 + P ({ω}).
2kXk∞
By construction, P ∗∗ is equivalent to P ∗ (same null-sets - actually, as P ∗ ∼ P
and P has no non-empty null-sets, neither do P ∗ , P ∗∗ ). As X is non-zero,
P ∗∗ and P ∗ are different. Now
E ∗∗ [ΣN
1 Hn .∆S̃ n ] = Σω P ∗∗
(ω) ΣN
1 Hn .∆S̃n (ω)
X(ω) ∗ N
= Σω 1 + P (ω) Σ1 Hn .∆S̃n (ω).
2kXk∞
The ‘1’ term on the right gives E ∗ [ΣN1 Hn .∆S̃n ], which is zero since this is a
∗
martingale transform of the E -martingale S̃n . The ‘X’ term gives a multiple
of the inner product
(X, ΣN 1 Hn .∆S̃n ),
9
a P ∗ equivalent to the true probability measure P is called an equivalent
martingale measure. Then
1 (No-Arbitrage Theorem: §2). If the market is viable (arbitrage-free),
equivalent martingale measures P ∗ exist.
2 (Completeness Theorem: §3). If the market is complete (all contingent
claims can be replicated), equivalent martingale measures are unique. Com-
bining:
V0 (H) = E ∗ [h/SN
0
].
More generally, the same argument gives Ṽn (H) = Vn (H)/Sn0 = E ∗ [(h/SN
0
)|Fn ]:
h
Vn (H) = Sn0 E ∗ [ 0
|Fn ] (n = 0, 1, · · · , N ).
SN
10
the empty set is null).
Now option pricing is our central task, and for pricing purposes P ∗ is
vital and P itself irrelevant. We thus may – and shall – focus attention on
P ∗ , which is called the risk-neutral probability measure. Risk-neutrality is
the central concept of the subject. The concept of risk-neutrality is due in
its modern form to Harrison and Pliska [HP] in 1981 – though the idea can
be traced back to actuarial practice much earlier. Harrison and Pliska call
P ∗ the reference measure; other names are risk-adjusted or martingale mea-
sure. The term ‘risk-neutral’ reflects the P ∗ -martingale property of the risky
assets, since martingales model fair games.
To summarise, we have the Risk-Neutral Valuation (or Pricing) Formula:
11
since Sn = S̃n (1 + r)n , Tn+1 = Sn+1 /Sn = (S̃n+1 /S̃n )(1 + r). But
p∗ = (b − r)/(b − a).
here j, N − n − j are the numbers of times Ti takes the two possible values
1 + a, 1 + b. This is the discrete Black-Scholes formula of Cox, Ross &
Rubinstein (1979) for pricing a European call option in the binomial model.
The European put is similar – or use put-call parity (II.3).
To find the (perfect-hedge) strategy for replicating this explicitly: write
N −n N − n
c(n, x) := Σj=0 p∗ j (1 − p∗ )N −n−j (x(1 + a)j (1 + b)N −n−j − K)+ .
j
12
Then c(n, x) is the undiscounted P ∗ -expectation of the call at time n given
that Sn = x. This must be the value of the portfolio at time n if the strategy
H = (Hn ) replicates the claim:
Subtract:
Notice that the numerator is the difference of two values of c(n, x) with
the larger value of x in the first term (recall b > a). When the payoff function
c(n, x) is an increasing function of x, as for the European call option consid-
ered here, this is non-negative. In this case, the Proposition gives Hn ≥ 0:
the replicating strategy does not involve short-selling. We record this as:
13
(temporarily) ρ ≥ 0 for the instantaneous interest rate in continuous time,
and define (again temporarily) r by
ρT N
r := ρT /N : eρT = limN →∞ (1 + ) = limN →∞ (1 + r)N .
N
Here r, which tends to zero as N → ∞, represents the interest rate in discrete
time for the approximating binomial model.
For σ > 0 fixed (σ 2 plays the role of a variance, corresponding in continuous
time to the volatility of the stock – below), define a, b (→ 0 as N → ∞) by
√ √
log((1 + a)/(1 + r)) = −σ/ N , log((1 + b)/(1 + r)) = σ/ N .
We now have a sequence of binomial models, for each of which we can price
options as in §5. We shall show that the pricing formula converges as N → ∞
to a limit. This is the famous Black-Scholes formula, the central result of
the course. We shall meet it later, and re-derive it, in continuous time, its
natural setting, in Ch. VII; see also e.g. [BK], 4.6.2. Fortunately, the con-
tinuous Black-Scholes formula is much neater than its discrete counterpart,
which involves the unwieldy binomial sum above.
N µN → µ (N → ∞)
YN → Y = N (µ, σ) (N → ∞).
Proof. Use characteristic functions (CFs), I.4: since YN has mean and vari-
ance as given, it also has second moment σ 2 (1 + o(1))/N , so has CF
φN (u) := E exp{iuYN } = ΠN N N N
1 E exp{iuXj } = [E exp{iuX1 }]
iuµ 1 σ 2 u2 1 1
= (1 + − + o( ))N → exp{iuµ − σ 2 u2 } (N → ∞),
N 2 N N 2
the CF of the normal law N (µ, σ). Convergence of CFs implies convergence
in distribution by Lévy’s continuity theorem for CFs ([W], §18.1). //
14
We can apply this to pricing the call option above (the details of the
calculation below, which are messy, are not hard, and not crucial; the point
is that this can be done):
(N ) ρT −N ∗
C0 = (1 + ) E [(S0 ΠN
1 Tn − K)+ ]
N
ρT −N
= E ∗ [(S0 exp{YN } − (1 + ) K)+ ], (1)
N
where XN
YN := log(Tn /(1 + r)).
1
so a, b, r → 0 as N → ∞, and
√
(b − r) [(1 + r)eσ/ N − 1 − r]
1 − 2p∗ = 1 − 2 =1−2 √ √
(b − a) [(1 + r)(eσ/ N − e−σ/ N )]
√
[eσ/ N
− 1]
=1−2 √ √ .
[eσ/ N − e−σ/ N ]
Now expand the two [· · · ] terms above by Taylor’s theorem: they give
σ 1 σ 2σ σ2
√ (1 + √ + · · · ), √ (1 + + · · · ).
N 2 N N 6N
√
So, cancelling σ/ N ,
1 √σ
∗
2(1 + 2 N
+ ···) 1 σ
1 − 2p = 1 − σ 2 = − √ + O(1/N ) :
2(1 + 6N + ···) 2 N
15
σ 1 σ 1
N µN = N. √ .(− √ + O(1/N )) → µ := − σ 2 (N → ∞).
N 2 N 2
We now need to change notation:
(i) We replace the variance σ 2 above by σ 2 T . So σ 2 is the variance per unit
time (which is more suited to the work of Ch. VI, VII in continuous time);
the standard deviation (SD) σ is called the volatility. It measures the vari-
ability of the stock, so its riskiness, or its sensitivity to new information.
(ii) We replace ρ in the above by r. This is the standard notation for the
riskless interest rate in continuous time, to which we are now moving.
As usual, we write the standard normal density function as φ and distri-
bution function as Φ:
1 2 Z x Z x − 1 u2
e− 2 x e 2
φ(x) := √ , Φ(x) := φ(u)du = √ du.
2π −∞ −∞ 2π
Note that as φ is even, the left and right tails of Φ are equal:
Z −x Z ∞
φ(x) = φ(−x), so φ(u)du = φ(u)du : Φ(−x) = 1−Φ(x).
−∞ x
where St is the stock price at time t ∈ [0, T ], K is the strike price, r is the
riskless interest rate, σ is the volatility and
1 √ √
d± := [log(S/K) + (r ± σ 2 )(T − t)]/σ T − t : d+ = d− + σ T − t.
2
For completeness, we state the corresponding Black-Scholes formula for
puts. The proofs of the two results are closely analogous, and one can derive
either from the other by put-call parity.
16
useful, as a benchmark and first approximation.
17
√ √
as −c + σ T = d− + σ T = d+ when t = 0. So the option price is given
in terms of the initial price S0 , strike price K, expiry T , interest rate r and
volatility σ by
1 √
C0 = S0 Φ(d+ ) − Ke−rT Φ(d− ), d± := [log(S/K) + (r ± σ 2 )T ]/σ T . //
2
Note. 1. Normal approximation to binomial. The proof above starts from a
binomial distribution and ends with a normal distribution. The binomial dis-
tribution is that of a sum of independent Bernoulli random variables. That
sums (or averages) of independent random variables with finite means and
variances gives a normal limit is the content of the Central Limit Theorem
or CLT (the Law of Errors, as physicists would say). This form of the CLT
is the de Moivre-Laplace limit theorem.
The picture for this is familiar. The Binomial distribution B(n, p) has a
histogram with n + 1 bars, whose heights peak at the mode and decrease to
either side. For large n, one can draw a smooth curve through the histogram.
This curve is the relevant approximating normal density.
2. The Cox-Ross-Rubinstein binomial model above goes over in the passage
to the limit to the geometric Brownian motion model of VII.1. We will later
re-derive the continuous Black-Scholes formula in Ch. VII, using continuous-
time methods (Itô calculus), rather than, as above, deriving the discrete BS
formula and going to the limit on the formula, rather than the model.
3. For similar derivations of the discrete Black-Scholes formula and the pas-
sage to the limit to the continuous Black-Scholes formula, see e.g. [CR], §5.6.
4. One of the most striking features of the Black-Scholes formula is that it
does not involve the mean rate of return µ of the stock – only the riskless
interest-rate r and the volatility of the stock σ. Mathematically, this reflects
the fact that the change of measure involved in the passage to the risk-neutral
measure involves a change of drift. This eliminates the µ term; see VII.
5. Volatility. The volatility σ can be estimated in two ways:
a. Historic volatility. Directly from the movement of a stock price in time,
using Time Series methods in discrete time [see Ch. VI for continuous time].
b. Implied volatility. From the observed market prices of options: if we know
everything in the Black-Scholes formula (including the price at which the
option is traded) except the volatility σ, we can solve for σ. Since σ appears
inside the argument of the normal distribution function Φ as well as outside
it, this is a transcendental equation for σ and has to be solved numerically by
iteration (Newton-Raphson method). We quote (see ‘The Greeks’ below, and
18
Problems 7) that the Black-Scholes price is a monotone (increasing) function
of the volatility (more volatility doesn’t make us ‘more likely to win’, but
when we do win, we ‘win bigger’), so there is a unique root of the equation.
In practice, one sees discrepancies between historic and implied volatility,
which show limitations to the accuracy of the Black-Scholes model. But it is
the standard ‘benchmark model’, and useful as a first approximation.
The classical view of volatility is that it is caused by future uncertainty,
and shows the market’s reaction to the stream of new information. How-
ever, studies taking into account periods when the markets are open and
closed [there are only about 250 trading days in the year] have shown that
the volatility is less when markets are closed than when they are open. This
suggests that trading itself is one of the main causes of volatility.
Note. This observation has deep implications for the macro-prudential and
regulatory issues discussed in Ch. II. The real economy cannot afford too
much volatility. Volatility is (at least partly) caused by trading. Conclusion:
there is too much trading. Policy question: how can we reduce the volume
of trading (much of it speculative, designed to enrich traders, and not serv-
ing a more widely useful economic purpose)? One answer is the so-called
Tobin tax (also known as the ”Robin Hood tax”) (James Tobin (1918-2002),
American economist; Nobel Prize for Economics, 1981). This would levy
a small charge (e.g. 0.01%) on all financial transactions. This would both
provide a major and useful source of tax revenue, and – more importantly –
would discourage a lot of speculative trading, thereby (shrinking the size of
the financial services industry, but) diminishing volatility, to the benefit of
the general economy (Problems 7 again).
If the Black-Scholes model were perfect, historic and implied volatility
estimates would agree (to within sampling error). But discrepancies can be
observed, which shows the imperfections of our model.
Volatility estimation is a major topic, both theoretically and in practice.
We return to this in V.7.3-4 below and VII.7.5-8. But we note here:
(i) trading is itself one of the major causes of volatility, as above;
(ii) options like volatility [i.e., option prices go up with volatility].
Recalling Ch. II, this shows that volatility is a ‘bad thing’ from the point
of view of the real economy (uncertainty about, e.g., future material costs
is nothing but a nuisance to manufacturers), but a ‘good thing’ for financial
markets (trading increases volatility, which increases option prices, which
generates more trade ...) – at the cost of increased instability.
19
§7. More on European Options
1. Bounds. We use the notation above. We also write c, p for the values of
European calls and puts, C, P for the values of the American counterparts.
Obvious upper bounds are c ≤ S, C ≤ S, where S is the stock price (we
can buy for S on the market without worrying about options, so would not
pay more than this for the right to buy). For puts, one has correspondingly
the obvious upper bounds p ≤ K, P ≤ K, where K is the strike price: one
would not pay more than K for the right to sell at price K, as one would not
spend more than one’s maximum return. For lower bounds:
c0 ≥ max(S0 − Ke−rT , 0).
Proof. Consider the following two portfolios:
I: one European call plus Ke−rT in cash; II: one share. Show ”I ≥ II”.
p0 ≥ max(Ke−rT − S0 , 0) (proof: by above and put-call parity).
2. Dependence of the Black-Scholes price on the parameters.
Recall the Black-Scholes formulae for the values ct , pt for the European
call and put: with
1 p
d± := [log(St /K) + (r ± σ 2 )(T − t)]/σ (T − t),
2
ct = St Φ(d+ ) − Ke−r(T −t) Φ(d− ), pt = Ke−r(T −t) Φ(−d− ) − St Φ(−d+ ),
(a). S. As the stock price S increases, the call option becomes more likely to
be exercised. As S → ∞, d± → ∞, Φ(d± ) → 1, so ct ∼ St − Ke−r(T −t) . This
has a natural economic interpretation: as the value of a forward contract
with delivery price K (Hull [H1] Ch. 3, [H2] Ch. 3).
(b). σ. When the volatility σ → 0, the stock becomes riskless, and behaves
like money in the bank. Again, d± → ∞, as above, with the same economic
interpretation.
3. The Greeks.
These are the partial derivatives of the option price with respect to the
input parameters. They have the interpretation of sensitivities.
(i) For a call, say, ∂c/∂S is called the delta, ∆. Adjusting our holdings of
stock to eliminate our portfolio’s dependence on S is called delta-hedging.
(ii) Second-order effects involve gamma := ∂(∆)/∂S.
(iii) Time-dependence is given by Theta is ∂c/∂t.
(iv) Volatility dependence is given by vega := ∂c/∂σ.1
1
Of course, vega is not a letter of the Greek alphabet! (it is the Spanish word for
‘meadow’, as in Las Vegas) – presumably so named for ”v for volatility, variance and
vega”, and because vega sounds quite like beta, etc.
20
Vega.
From the Black-Scholes formula (which gives the price explicitly as a
function of σ), one can check by calculus (Problems 7) that
∂c/∂σ > 0,
and similarly for puts (or, use the result for calls and put-call parity). In sum:
options like volatility. This fits our intuition. The more uncertain things are
(the higher the volatility), the more valuable protection against adversity –
or insurance against it – becomes (the higher the option price).
(v) rho is ∂c/∂r, the sensitivity to interest rates.
C = c.
St − K + Ke−r(T −t) .
Since r > 0 and t < T , this is < St , the value of Portfolio II at t. So Portfolio
I is always worth less than Portfolio II if exercised early.
If however the option is exercised instead at expiry, T , the American call
option is then the same as a European call option. Then at time T , Portfolio
I is worth max(ST , K) and Portfolio II is worth ST . So:
21
This direct comparison with the underlying [the share in Portfolio II] shows
that early exercise is never optimal. Since an American option at expiry is
the same as a European one, this completes the proof. //
Second Proof. One can instead use the bounds of §7.1. For details, see e.g.
[BK, Th. 4.7.1].
Financial Interpretation.
There are two reasons why an American call should not be exercised early:
1. Insurance. Consider an investor choosing to hold a call option instead of
the underlying stock. He does not care if the share price falls below the strike
price (as he can then just discard his option) – but if he held the stock, he
would. Thus the option insures the investor against such a fall in stock price,
and if he exercises early, he loses this insurance.
2. Interest on the strike price. When the holder exercises the option, he
buys the stock and pays the strike price, K. Early exercise at t < T loses the
interest on K between times t and T : the later he pays out K, the better.
Economic Note. Despite Merton’s theorem, and the interpretation above,
there are plenty of real-life situations where early exercise of an American
call might be sensible, and indeed it is done routinely. Consider, for exam-
ple, a manufacturer of electrical goods, in bulk. He needs a regular supply
of large amounts of copper. The danger is future price increases; the obvious
precaution is to hedge against this by buying call options. If the expiry is a
year but copper stocks are running low after six months, he would exercise
his American call early, to keep an adequate inventory of copper, his crucial
raw material. This ensures that his main business activity – manufacturing
– can continue unobstructed. Neither of the reasons above applies here:
Insurance. He doesn’t care if the price of copper falls: he isn’t going to sell
his copper stocks, but use them.
Interest. He doesn’t care about losing interest on cash over the remaining six
months. He is in manufacturing to use his money to make things, and then
sell them, rather than put it in the bank.
This neatly illustrates the contrast between finance (money, options etc.)
and economics (the real economy – goods and services).
Put-Call Symmetry.
The BS formulae for puts and calls resemble each other, with stock price
S and discounted strike K interchanged. Results of this type are called put-
call symmetry.
22
American Puts.
Recall the put-call parity of Ch. II (valid only for European options):
c − p = S − Ke−rT . A partial analogue for American options is given by the
inequalities below:
For proof (as above) and background, see e.g. Ch. 8 (p. 216) of [H1].
We now consider how to evaluate an American put option, European
and American call options having been treated already. First, we will need
to work in discrete time. We do this by dividing the time-interval [0, T ]
into N equal subintervals of length ∆t say. Next, we take the values of the
underlying stock to be discrete: we use the binomial model of §5, with a
slight change of notation: we write u, d (‘up’, ‘down’) for (1 + b), (1 + a):
thus stock with initial value S is worth Sui dj after i steps up and j steps
down. Consequently, after N steps, there are N + 1 possible prices, Sui dN −i
(i = 0, · · · , N ). It is convenient to display the possible paths followed by the
stock price as a binomial tree, with time going left to right and two paths, up
and down, leaving each node in the tree, until we reach the N + 1 terminal
nodes at expiry. There are 2N possible paths through the tree. It is common
to take N of the order of 30, for two reasons:
(i) typical lengths of time to expiry are measured in months (9 months, say);
this gives a time-step around the corresponding number of days,
(ii) 230 paths is about the order of magnitude that can be easily handled by
computers (recall that 210 = 1, 024, so 230 is somewhat over a billion).
We now return to the binomial model in §§5,6, with a slight change of
notation. Recall that in §5 (discrete time) we used 1 + r for the discount
factor. Now call this 1 + ρ instead, freeing r for its usual use as the short
rate of interest in continuous time. Thus 1 + ρ = er∆t , and the risk-neutrality
condition p∗ = (b − r)/(b − a) of §5 becomes
23
√
Since ∆t is small, its square ∆t is a second-order term. So to first order,
ud = 1, which simplifies filling in the terminal values in the binary tree.
We begin again: define our up and down factors u, d so that
ud = 1; (∗)
p∗ = (u − er∆t )/(u − d)
(to get the mean return from the risky stock the same as that from the
riskless bank account), and the volatility σ to get the variance of the stock
price S 0 after one time-step when it is worth S initially as S 2 σ 2 ∆t:
S 2 σ 2 ∆t = p∗ S 2 u2 + (1 − p∗ )S 2 d2 − S 2 [p∗ u + (1 − p∗ )d]2
√
(using varS 0 = E(S 0 2 ) − [ES 0 ]2 ). Then to first order in ∆t (which is all the
accuracy we shall need), one can check that we have as before
√ √
u = exp(σ ∆t), d = exp(−σ ∆t). (∗∗)
We can now calculate both the value of an American put option and
the optimal exercise strategy by working backwards through the tree (this
method of backward recursion in time is a form of the Dynamic Programming
[DP] technique (Richard Bellman (1920-84) in 1953, book, 1957), which is
important in many areas of optimization and Operational Research (OR)).
1. Draw a binary tree showing the initial stock value and having the right
number, N , of time-intervals.
2. Fill in the stock prices: after one time interval, these are Su (upper) and
Sd (lower); after two time-intervals, Su2 , S and Sd2 = S/u2 ; after i time-
intervals, these are Suj di−j = Su2j−i at the node with j ‘up’ steps and i − j
‘down’ steps (the ‘(i, j)’ node).
3. Using the strike price K and the prices at the terminal nodes, fill in the
payoffs (fN,j = max[K − Suj dN −j , 0]) from the option at the terminal nodes
(where, at expiry, the values of the European and American options coincide)
underneath the terminal prices.
4. Work back down the tree one time-step. Fill in the ‘European’ value at
the penultimate nodes as the discounted values of the upper and lower right
(terminal node) values, under the risk-neutral measure – ‘p∗ times lower right
24
plus 1−p∗ times upper right’ [notation of V.6]. Fill in the ‘intrinsic’ (or early-
exercise) value - the value if the option is exercised. Fill in the American put
value as the higher of these.
5. Treat these values as ‘terminal node values’, and fill in the values one
time-step earlier by repeating Step 4 for this ‘reduced tree’.
6. Iterate. The value of the American put at time 0 is the value at the root -
the last node to be filled in. The ‘early-exercise region’ is the node set where
the early-exercise value is the higher; the rest is the ‘continuation region’.
Note. The above procedure is simple to describe and understand, and simple
to programme. It is laborious to implement numerically by hand, on exam-
ples big enough to be non-trivial. Numerical examples are worked through
in detail in [H1], 359-360 and [CR], 241-242.
Mathematically, the task remains of describing the continuation region -
the part of the tree where early exercise is not optimal. This is a classical
optimal stopping problem. No explicit solution is known (and presumably
there isn’t one). We will, however, connect the work above with that of IV.7
on the Snell envelope. Consider the pricing of an American put, strike price
K, expiry N , in discrete time, with discount factor 1 + r per unit time as
earlier. Let Z = (Zn )N n=0 be the payoff on exercising at time n. We want
to price Zn , by Un say (to conform to our earlier notation), so as to avoid
arbitrage; again, we work backwards in time. The recursive step is
1
Un−1 = max(Zn−1 , E ∗ [Un |Fn−1 ]),
1+r
the first alternative on the right corresponding to early exercise, the second
to the discounted expectation under P ∗ , as usual. Let Ũn = Un /(1 + r)n be
the discounted price of the American option. Then
(Ũn ) is the Snell envelope (III.7) of the discounted payoff process (Z̃n ). So:
(i) a P ∗ -supermartingale,
(ii) the smallest supermartingale dominating (Z̃n ),
(iii) the solution of the optimal stopping problem for Z̃.
Note. One can use the Snell envelope to prove Merton’s theorem (equiva-
lence of American and European calls) without using arbitrage arguments.
For details see e.g. [BK, Th. 4.7.1 and Cor. 4.7.1].
25
P -measure and P ∗ − (or Q−) measure.
We use P and P ∗ in the above, as E and E ∗ are convenient, but P and
Q when the emphasis is on Q, for brevity.
The measure P , the real (or real-world) probability measure, models the
uncertainty driving prices, which are indeed uncertain, thus allowing us to
bring mathematics to bear on financial problems. But P is difficult to get
at directly. By contrast, Q is more accessible: the market tells us about Q,
or more specifically, trading does. In addition, trading also tells us about
the volatility σ, via implied volatility, which we can infer from observing the
prices at which options are traded. So Q is certainly more accessible than P .
There is thus a sense in which it is Q, rather than P , which is the more real.
It is as well to bear all this in mind when looking at specific problems, par-
ticularly numerical ones. Now that we know the CRR binomial-tree model,
which gives us the Black-Scholes formula in discrete time (and hence also, by
the limiting argument above, the Black- Scholes formula in continuous time,
the main result of the course), we can recognise the ‘one-period, up or down’
model ($/SFr in II.8, price of gold in Problems 5), though clearly artificial
and stylised, as a workable ‘building block’ of the whole theory. Because P
itself does not occur in the Black-Scholes formula(e), from a purely financial
point of view there is little need to try to construct more realistic, and so
more complicated, models of P . Instead, one can exploit what one can infer
about Q, which does occur in Black-Scholes, from seeing the prices at which
options trade.
From the economic point of view, it is the real world, the real economy,
and so the real probability measure P , that matters. The ‘Q-measure-eye
view of the world’ has a degree of artificiality, in so far as options do. One
can eat food, and needs to. One can’t eat options.
A fuller discussion of Q-measure involves Arrow-Debreu prices, equilibria
etc., but we omit this for lack of time, and because it would take us too far
into Economics.
26