Lectures On Financial Mathematics: Harald Lang
Lectures On Financial Mathematics: Harald Lang
Harald Lang
c Harald Lang, KTH Mathematics 2012
Preface
Preface
My main goal with this text is to present the mathematical modelling
of financial markets in a mathematically rigorous way, yet avoiding math-
ematical technicalities that tends to deter people from trying to access
it.
Trade takes place in discrete time; the continuous case is considered
as the limiting case when the length of the time intervals tend to zero.
However, the dynamics of asset values are modelled in continuous time as
in the usual Black-Scholes model. This avoids some mathematical techni-
calities that seem irrelevant to the reality we are modelling.
The text focuses on the price dynamics of forward (or futures) prices
rather than spot prices, which is more traditional. The rationale for this is
that forward and futures prices for any good—also consumption goods—
exhibit a Martingale property on an arbitrage free market, whereas this is
not true in general for spot prices (other than for pure investment assets.)
It also simplifies computations when derivatives on investment assets that
pay dividends are studied.
Another departure from more traditional texts is that I avoid the no-
tion of “objective” probabilities or probability distributions. I think they
are suspect constructs in this context. We can in a meaningful way assign
probabilities to outcomes of experiments that can be repeated under simi-
lar circumstances, or where there are strong symmetries between possible
outcomes. But it is unclear to me what the “objective” probability distri-
bution for the price of crude oil, say, at some future point in time would
be. In fact, I don’t think this is a well defined concept.
The text presents the mathematical modelling of financial markets. In
order to get familiar with the workings of these markets in practice, the
reader is encouraged to supplement this text with some text on financial
economics. A good such text book is John C. Hull’s: Options, Futures, &
Other Derivatives (Prentice Hall,) which I will refer to in some places.
13/9–2012
Harald Lang
Contents
Contents
I: Introduction to Present-, Forward- and Futures Prices . . . 1
Zero Coupon Bonds . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
Money Market Account . . . . . . . . . . . . . . . . . . . . . . . . . . 2
Relations between Present-, Forward- and Futures Prices . . . . . 3
Comparison of Forward- and Futures Prices . . . . . . . . . . . . . 4
Spot Prices, Storage Cost and Dividends . . . . . . . . . . . . . . . 6
Comments . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
II: Forwards, FRA:s and Swaps . . . . . . . . . . . . . . . . . . . 8
Forward Prices . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
Forward Rate Agreements . . . . . . . . . . . . . . . . . . . . . . . . 10
Plain Vanilla Interest Rate Swap . . . . . . . . . . . . . . . . . . . . 11
Exercises and Examples . . . . . . . . . . . . . . . . . . . . . . . . . 12
Solutions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
III: Optimal Hedge Ratio . . . . . . . . . . . . . . . . . . . . . . 17
Exercises and Examples . . . . . . . . . . . . . . . . . . . . . . . . . . 19
Solutions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
IV: Conditions for No Arbitrage . . . . . . . . . . . . . . . . . 21
Theorem (The No Arbitrage Theorem) . . . . . . . . . . . . . . . . 22
The No Arbitrage Assumption . . . . . . . . . . . . . . . . . . . . . . 23
V: Pricing European Derivatives . . . . . . . . . . . . . . . . . 25
Black’s Model . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
The Black-Scholes Pricing Formula . . . . . . . . . . . . . . . . . . . 26
Put and Call Options . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
The Interpretation of σ and the Market Price of Risk . . . . . . . 26
Comments . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
Exercises and Examples . . . . . . . . . . . . . . . . . . . . . . . . . 28
Solutions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
VI: Yield and Duration . . . . . . . . . . . . . . . . . . . . . . . 32
Forward Yield and Forward Duration . . . . . . . . . . . . . . . . . 34
Black’s Model for Bond Options . . . . . . . . . . . . . . . . . . . . 36
Portfolio Immunising . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
Exercises and Examples . . . . . . . . . . . . . . . . . . . . . . . . . . 37
Solutions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
VII: Risk Adjusted Probability Distributions . . . . . . . . . 42
An Example . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 43
Forward Distributions for Different Maturities . . . . . . . . . . . . 44
Exercises and Examples . . . . . . . . . . . . . . . . . . . . . . . . . . 46
Solutions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
VIII: Conditional Expectations and Martingales . . . . . . . 48
Martingale Prices . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
Contents
Forward contract
day 0 1 2 3 ··· T−1 T
cash flow 0 0 0 0 ··· 0 X −G
1
I: Present-, Forward and Futures Prices
Futures contract
day 0 1 2 3 ··· T−1 T
cash flow 0 F1 −F0 F2 −F1 F3 −F2 ··· FT−1 −FT−2 X−FT−1
The simplest of these three contracts is the one when we pay in advance,
at least if the good that is delivered is non-pecuniary, since in that case
the interest does not play a part. For futures contracts, the interest rate
clearly plays a part, since the return of the contract is spread out over
time.
We will derive some book-keeping relations between the present prices,
forward prices and futures prices, but first we need some interest rate
securities.
We need a notation for the currency, and we use a dollar sign $, even
though the currency may be Euro or any other currency.
2
I: Present-, Forward and Futures Prices
Theorem
The following relations hold:
(T ) (T ) (T )
a) P0 , G0 and F0 are linear functions, i.e., if X and Y are random
payments made at time T , then for any constants a and b
(T ) (T ) (T )
P0 [aX + bY ] = aP0 [X] + bP0 [Y ],
(T ) (T )
and similarly for G0 and F0 .
(T )
b) For any deterministic (i.e., known today) value V, G0 [V ] = V ,
(T ) (T )
F0 [V ] = V and P0 [V ] = ZT V
(T ) (T )
c) P0 [X] = ZT G0 [X]
(T ) (T )
d) P0 [XeR(0,T ) ] = F0 [X]
(T ) (T )
e) P0 [X] = F0 [Xe−R(0,T ) ]
Proof
The proof relies on an assumption of the model: the law of one price.
It means that there can not be two contracts that both yield the same
payoff X at time T , but have different prices today. Indeed, if there were
two such contracts, we would buy the cheaper and sell the more expensive,
and make a profit today, and have no further cash flows in the future. But
so would everyone else, and this is inconsistent with market equilibrium.
In Ch. IV we will extend this model assumption somewhat.
To prove c), note that if we take a long position on a forward contract
on X and at the same time a long position on a zero coupon bond with
(T ) (T )
face value G0 [X], then we have a portfolio which costs ZT G0 [X] today,
and yields the income X at time T . By the law of one price, it hence must
be that c) holds.
To prove d), consider the following strategy: Deposit F0 on the money
market account, and take er1 long positions on the futures contract on X
for delivery at time T .
The next day the total balance is then F0 er1 + er1 (F1 − F0 ) = F1 er1 .
Deposit this on the money market account, and increase the futures posi-
tion to er1 +r2 contracts.
3
I: Present-, Forward and Futures Prices
The next day, day 2, the total balance is then F1 er1 +r2 + er1 +r2 (F2 −
F1 ) = F2 er1 +r2 . Deposit this on the money market account, and increase
the futures position to er1 +r2 +r3 contracts.
And so on, up to day T when the total balance is FT er1 +···+rT =
Xer1 +···+rT . In this way we have a strategy which is equivalent to a
contract where we pay F0 day zero, and receive the value FT er1 +···+rT =
Xer1 +···+rT day T . This proves d).
Since relation d) is true for any random variable X whose outcome
is known day T , we may replace X by Xe−R(0,T ) in that relation. This
proves e).
(T )
It is now easy to prove b). The fact that P0 [V ] = ZT V is simply
(T )
the definition of ZT . The relation G0 [V ] = V now follows from c) with
(T )
X = V . In order to prove that F0 [V ] = V , note that by the definition
of money market account, the price needed to be paid day zero in order
(T )
to receive V eR(0,T ) day T is V ; hence V = P0 [V eR(0,T ) ]. The relation
(T )
F0 [V ] = V now follows from d) with X = V .
(T )
Finally, to prove a), note that if we buy a contracts which cost P0 [X]
day zero and gives the payoff X day T , and b contracts that gives payoff Y ,
(T ) (T )
then we have a portfolio that costs aP0 [X]+bP0 [Y ] day zero and yields
(T ) (T ) (T )
the payoff aX + bY day T ; hence P0 [aX + bY ] = aP0 [X] + bP0 [Y ].
The other two relations now follow immediately employing c) and d). This
completes the proof.
(T ) (T )
$1 = P0 [$eR(0,T ) ] = P0 [$1]eR(0,T ) = $ZT eR(0,T ) ,
hence
ZT = e−R(0,T ) .
Therefore
(T ) (T ) (T )
ZT G0 [X] = P0 [X] = F0 [Xe−R(0,T ) ]
(T ) (T )
= F0 [X]e−R(0,T ) = F0 [X]ZT ,
and hence
(T ) (T )
G0 [X] = F0 [X].
4
I: Present-, Forward and Futures Prices
Corollary
If interest rates are deterministic, the forward price and the futures
(T ) (T )
price coincide: G0 [X] = F0 [X]
The equality of forward- and futures prices does not in general hold if
interest rates are random, though. To see this, we show as an example
(T ) (T )
that if eR(0,T ) is random, then F0 [$eR(0,T ) ] > G0 [$eR(0,T ) ].
Indeed, note that the function y = x1 is convex for x > 0. This implies
1
that its graph lies over its tangent. Let y = m + k(x − m) be the tangent
1 1 1
line through the point (m, m ). Then x ≥ m + k(x − m) with equality
only for x = m (we consider only positive values of x.) Now use this with
x = e−R(0,T ) and m = ZT . We then have
where the equality holds only for one particular value of R(0, T ). In the
absence of arbitrage (we will come back to this in Ch. IV,) the futures
price of the value of the left hand side is greater than the futures price of
the value of the right hand side, i.e.,
(T ) (T )
F0 [$eR(0,T ) ] > $ZT−1 + k(F0 [$e−R(0,T ) ] − $ZT )
(T )
But F0 [$e−R(0,T ) ] = P (T ) [$1] = $ZT , so the parenthesis following k is
equal to zero, hence
(T )
F0 [$eR(0,T ) ] > $ZT−1
(T ) (T )
ZT G0 [$eR(0,T ) ] = P0 [$eR(0,T ) ] = $1
(T )
so G0 [$eR(0,T ) ] = $ZT−1 , and it follows that
(T ) (T )
F0 [$eR(0,T ) ] > G0 [$eR(0,T ) ].
5
I: Present-, Forward and Futures Prices
Comments
You can read about forward and futures contracts John Hull’s book
“Options, Futures, & other Derivatives 1 ”. He describes in detail how fu-
tures contracts work, and why they are specified in the somewhat peculiar
way they are.
The mathematical modelling of a futures contract is a slight simplifi-
cation of the real contract. We disregard the maintenance account, thus
avoiding any problems with interest on the balance. Furthermore, we as-
sume that the delivery date is defined as a certain day, not a whole month.
We also disregard the issues on accounting and tax.
We use “continuous compounding” of interest rates (use of the expo-
nential function.) If you feel uncomfortable with this, you may want to
read relevant chapters in Hull’s book. We will use continuous compound-
ing unless otherwise explicitly stated, since it is the most convenient way
to handle interest. It is of course easy to convert between continuous
compounded interest and any other compounding.
1
Prentice Hall
6
I: Present-, Forward and Futures Prices
It is important not to confuse the present price with the spot price
of the same type of good. The spot price is the price of the good for
immediate delivery.
It is also extremely important to distinguish between constants (val-
ues that are currently known) and random variables. For instance, as-
sume that ZT = e−rT (where r hence is a number.) Then it is true that
(T ) (T ) (T )
P0 [X] = e−rT G0 [X] (Theorem c,) however, the relation P0 [X] =
(T )
e−R(0,T ) G0 [X] is invalid and nonsense! Indeed, R(0, T ) is a random
(T ) (T )
variable; its outcome is not known until time T−1, whereas G0 and P0
are known prices today.
7
II: Forwards, FRA:s and Swaps
Forward Prices
In many cases the theorem of Ch. I can be used to calculate forward
prices. As we will see later, in order to calculate option prices, it is essential
to first calculate the forward price of the underlying asset.
Example 1.
Consider a share of a stock which costs S0 today, and which gives
a known dividend amount d in t years, and whose (random) spot price
at time T > t is ST . Assume that there are no other dividends or other
convenience yield during the time up to T . What is the forward price G
on this stock for delivery at time T ?
Assume that we buy the stock today, and sell it at time T . The cash
flow is
day 0 t T
cash flow −S0 d ST
The present value of the dividend is Zt d and the present value of the
income ST at time T is Zt G. Hence we have the relation
S0 = Zt d + ZT G
Example 2.
Consider a share of a stock which costs S0 now, and which gives a
known dividend yield d St i t years, where St is the spot price immediately
before the dividend is paid out. Let the (random) spot price at time T > t
be ST . Assume that there are no other dividends or other convenience
yield during the time up to T . What is the forward price G on this stock
for delivery at time T ?
Consider the strategy of buying the stock now, and sell it at time t
immediately before the dividend is paid out.
day 0 t
cash flow −S0 St
With the notation of Ch. I, we have the relation
(t)
P0 [St ] = S0 (1)
Consider now the strategy of buying the stock now, cash the dividend at
time t, and eventually sell the stock at time T .
8
II: Forwards, FRA:s and Swaps
day 0 t T
cash flow −S0 dSt ST
The present value of the dividend is d P (t) [St ] and the present value of the
income ST at time T is ZT G. Hence we have the relation
(t)
S0 = d P0 [St ] + ZT G
(1 − d)S0 = ZT G
Example 3.
With the same setting as in example 2, assume that there are dividend
yield payments at several points in time t1 < . . . < tn , where tn < T, each
time with the amount d Stj . As in the above example, we can buy the
stock today and sell it just before the first dividend is paid out, so the
relation
(t1 )
P0 [St1 ] = S0 (2)
holds. For any k = 2, 3, . . . , n we can buy the stock at time tk−1 imme-
diately before the payment of the dividend, collect the dividend, and sell
the stock immediately before the dividend is paid out at time tk .
day 0 tk−1 tk
cash flow 0 −Stk−1 + dStk−1 Stk
The price of this strategy today is zero, so we have
(tk−1 ) (tk−1 ) (tk )
0 = −P0 [Stk−1 ] + d P0 [Stk−1 ] + P0 [Stk ], i.e.,
(t ) (t )
P0 k [Stk ] = (1 − d)P0 k−1 [Stk−1 ]
(T ) (tn )
and repeated use of this relation and P0 [ST ] = (1 − d)P0 [Stn ] gives
(T ) (t1 )
P0 [ST ] = (1 − d)n P0 [St1 ] = (1 − d)n S0
where we have used (2) to obtain the last equality. Hence, by the theorem
of Ch. I, we have the relation
ZT G = (1 − d)n S0 (3)
9
II: Forwards, FRA:s and Swaps
Example 4.
We now consider the setting in example 3, but with a continuous
dividend yield ρ, i.e., for any small interval in time (t, t + δt) we get the
dividend ρSt δt. If we divide the time interval (0, T ) into a large number
n of intervals of length δt = T /n, we see from (3) that
ZT G = (1 − ρ δt)n S0
ZT G = e−ρT S0
Derivation
of the limit: ln (1 −ρ δt)n = n ln(1 − ρ δt) = n − ρ δt +
O(δt2 ) = n − ρT /n + O(1/n2 ) → −ρT. Taking exponential gives the
limit.
Example 5.
Assume we want to buy a foreign currency in t years time at an
exchange rate, the forward rate, agreed upon today. Assume that the
interest on the foreign currency is ρ and the domestic rate is r per year.
Let X0 be the exchange rate now (one unit of foreign currency costs X0 in
domestic currency,) and let Xt be the (random) exchange rate as of time
t. Let G be the forward exchange rate.
Consider now the strategy: buy one unit of foreign currency today,
buy foreign zero coupon bonds for the amount, so that we have eρt worth
of bonds in foreign currency at time t when we sell the bonds.
day 0 t
ρt
cash flow −X0 e Xt
Since the exchange rate at that time is Xt , we get Xt eρt in domestic
currency. Since the price we have paid today is X0 we have the relation
(t)
X0 = P 0 Xt eρt = e−rt G eρt = G e(ρ−r)t
i.e.,
G = X0 e(r−ρ)t
10
II: Forwards, FRA:s and Swaps
0 = Zt L − ZT Lef (T −t) .
From this we can solve for f . The interest rate f is the forward rate from
t to T .
The floating rate between tj and tk is the zero coupon rate that prevails
between these two points in time. The amount that A pays at time tk is
thus Lk · (1/Z(tk−1 , tk ) − 1), where Z(tk−1 , tk ) of course is the price of the
zero coupon bond at time tk−1 that matures at tk . Note that this price
is unknown today but known at time tk−1 . The total amount that B will
receive, and A will pay is thus random.
On the other hand, party B pays A a fixed amount c at each of the
times t1 , . . . , tn . The question is what c ought to be in order to make this
deal “fair”.
The notation with continuous compounding is here bit inconvenient.
Let us introduce the one period floating rate r̂j : the interest from time
tj−1 to time tj , i.e., if I deposit an amount a on a bank account at time
tj−1 the balance at time tj is a + ar̂j . This is the same as to say that a
zero coupon bond issued at time tj−1 with maturity at tj is 1/(1 + r̂j ).
Note that r̂j is random whose value becomes known at time tj−1 . The
cash flow that A pays to B is then
day t0 t1 t2 t3 ··· tn
cash flow 0 r̂1 L1 r̂2 L2 r̂3 L3 ··· r̂n Ln
In order to calculate the present value of this cash flow, we first determine
(t )
the present value P0 k (r̂k Lk ).
Consider the strategy: buy Lk worth of zero coupon bonds at time
tk−1 with maturity at tk and face value (1 + r̂k )Lk . This costs nothing
today, so the cash flow is
11
II: Forwards, FRA:s and Swaps
(t ) (tk ) (tk )
hence 0 = −P0 k−1 [Lk ] + P0 [Lk ] + P0 [r̂k Lk ] = −Ztk−1 Lk + Ztk Lk +
(t )
P0 k [r̂k Lk ], so we get
(tk )
P0 [r̂k Lk ] = (Ztk−1 − Ztk )Lk .
It is now easy to calculate the present value of the cash flow from A to B:
it is
n
PAB = (Ztk−1 − Ztk )Lk
1
4. The exchange rate of US dollars is today 8.50 SEK per dollar. The
forward price of a dollar to be delivered in six months is 8.40 SEK.
If the Swedish six month interest rate is 4% a year, determine the
American six month interest rate. (6.367%)
12
II: Forwards, FRA:s and Swaps
(94.05 EUR)
b) If the share is worth 45 dollars six months later, what is the value
of the original forward contract at this time? If another forward
contract is to be written with the same date of maturity, what
should the forward price be? (2.949 dollars, 47.307 dollars)
13
II: Forwards, FRA:s and Swaps
9. Let ri be the random daily rate of interest (per day) from day i − 1
to day i, and R(0, t) = r1 + · · · + rt . The random variable Xt is the
stock exchange index day t (today is day 0.) The random variables
Xt and R(0, t) are not independent.
The forward price of a contract for delivery of the payment
Xt eR(0,t) EUR day t is 115 EUR, a zero coupon bond which pays
1 EUR day t costs 0.96 EUR.
Determine the futures price of a contract for delivery of Xt EUR
day t. The stock exchange index is today X0 = 100.
10. A share of a stock currently costs 80 EUR. One year from now, it will
pay a dividend of 5% of its price at the time of the payment, and the
same happens two years from now. Determine the forward price of
the asset for delivery in 2.5 years. The interest rate for all maturities
is 6% per year. (83.88 EUR)
11. The 6 month zero rate is 5% per year and the one year rate is 5.2%
per year. What is the forward rate from 6 to 12 months? (5.4% per
year)
12. Show that the present value of a cash flow at time T that equals the
floating interest from t to T (t < T ) on a principal $L is the same
as if the floating rate is replaced by the current forward rate. (Note
that this is an easy way to value interest rate swaps: just replace any
floating rate by the corresponding forward rate.)
Solutions
3. (The problem is admittedly somewhat artificial, but serves as an ex-
ercise.)
If we buy the stock today and sell it after one year, the cash flow can
be represented:
month 0 4 10 12
cash flow −S0 2 0.02S10 S12
i.e.,
(10) (12)
S0 = 2e−0.06·4/12 + 0.02 P0 S10 + P0 S12 . (1)
On the other hand, if we sell the stock after 10 months, before the
dividend, then the cash flow is
month 0 4 10
cash flow −S0 2 S10
i.e.,
14
II: Forwards, FRA:s and Swaps
(10)
S0 = 2e−0.06·4/12 + P0 S10
(10)
If we here solve for P0 S10 and substitute into (1) we get
(12)
0.98S0 = 1.96e−0.06·4/12 + P0 S12
(12) (12)
But P0 S12 = e−0.06 G0 S12 = $e−0.06 110, hence
6. The present value P0 of the cash flow of the underlying bond after 2
years (that is, after the bond has been delivered) is
Note that we don’t need the interest rates for shorter duration than
two years!
7a. Since there are no dividends or other convenience yield, the present
(1) (1) (1)
price P0 is the same as the spot price: P0 = $40. Hence G0 =
(1)
Z1−1 P0 = e0.1 · 40 ≈ $44.207.
15
II: Forwards, FRA:s and Swaps
(12)
b. By the same argument, six months later the present price P6 S12
= $45; S12 is the (random) spot price of the stock at 12 months.
The cash flow of the forward contract is S12 − $44.207 at 12 months,
(12)
so the present value of this cash flow at 6 months is P6 S12 −
e−0.1·0.5 44.207 ≈ $2.949.
If a forward contract were drawn up at six months, the forward
(12) −1
price would be G6 = Z6,12 S6 = e0.1·0.5 45 ≈ $47.303, given that
the interest rate is still 10% per year.
P0 = $(Zt − ZT ) L
The
deterministic cash flow at T when the forward rate is applied is
$ ef (T −t) − 1 L where f is the forward rate, and its present value is
of course
P0 = $ZT ef (T −t) − 1 L = $ ZT ef (T −t) − ZT L
P0 = $(Zt − ZT ) L = P0
16
III: Optimal Hedge Ratio
A company that knows that it is due to buy an asset in the future can
hedge by taking a long futures position on the asset. Similarly, a firm that
is going to sell may take a short hedge. But it may happen that there is
no futures contract on the market for the exact product or delivery date.
For instance, the firm might want to buy petrol or diesel, whereas the
closest futures contract is on crude oil, or the delivery date of the futures
contract is a month later than the date of the hedge. In this case one
might want to use several futures on different assets, or delivery times,
to hedge. Let S be the price of the asset to be hedged at the date t of
delivery, and Ft1 , . . . , Ftn the futures prices at date t of the contracts that
are being used to hedge. All these prices are of course random as seen
from today, but we assume that there are enough price data so that it is
possible to estimate their variances and covariances.
Assume that for each unit of volume of S we use futures contracts
corresponding to βi units of volume for contract F i . The difference be-
tween the spot price S and n the total futures price at the date of the end
of the hedge is e = S − i=1 βi Fti , or
n
S= βi Fti + e
i=1
n
The total price
paid for the asset including the hedge is S − i=1 βi (Fti −
n
F0i ) = e + i=1 βi F0i , where only e is random, i.e., unknown to us to-
day. The task is to choose βi such that the variance of the residual e is
minimised.
Lemma
If we choose the coefficients βi such that Cov(Fti , e) = 0 for i =
1, . . . , n, then the variance Var(e) is minimised.
Proof
Assume that we have chosen βi such that Cov(Fti , e) = 0 for i =
1, . . . , n, and consider any other choice of coefficients:
n
S= γi Fti + f
i=1
17
III: Optimal Hedge Ratio
n
n
n
f =S− γi Fti = βi Fti + e − γi Fti
i=1 i=1 i=1
n
= δi Fti + e
i=1
n
Note that since Cov(Fti , e) = 0 it holds that Cov( i=1 δi Fti , e) = 0, hence
n
Theorem
The set of coefficients βi that minimises the variance Var(e) of the
residual is the solution to the system
n
Cov(Ftj , Fti )βi = Cov(Ftj , S) j = 1, . . . , n.
i=1
Proof
The condition that Cov(Ftj , e) = 0 means that
n
0 = Cov(Ftj , S − βi Fti )
i=1
n
= Cov(Ftj , S) − βi Cov(Ftj , Fti )
i=1
Cov2 (S, Ft )
Var(S) = β 2 Var(Ft ) + Var(e) = + Var(e)
Var(Ft )
= ρ2 (S, Ft ) Var(S) + Var(e)
18
III: Optimal Hedge Ratio
where ρ(S, Ft ) = Cov(S, Ft )/ Var(S) Var(Ft ) is the correlation coeffi-
cient between S and Ft .
Solving
for Var(e) we get the pleasant relation Var(e) = Var(S) 1 −
ρ2 (S, Ft ) which means that the standard deviation of the hedged position
is 1 − ρ2 (S, Ft ) times that of the unhedged position.
Solutions
1. The optimal hedge ratio β is the solution to
Cov(F, F ) β = Cov(F, S)
where S = spot price and F = futures price at the time of the hedge.
Since Cov(F, F ) = Var(F ) we get (SD = standard deviation; ρ =
correlation coefficient)
19
III: Optimal Hedge Ratio
2. Let X be the value of one share some months from now, X0 = $30
is the current value. Likewise, let I be the value of the index some
months from now, I0 = 1’500 is the current value. The beta value β
is by definition such that
X I
=β +e
X0 I0
where the covariance Cov(I, e) = 0. In other words,
X0
X=β I + X0 e
I0
The idea of the strategy is to hedge the term containing I; this is the
impact on the value of the share due to general market movements.
We want to capitalise on the idiosyncratic term containing e.
Assume now that we own one share and short k futures. Each
future is for delivery of n = $50 times the index. The value of the
hedged position after some months is then
X0
X − knI = β − kn I + X0 e
I0
X0
so, obviously, in order to get rid if the I-term we should have β I0 −
kn = 0, i.e.,
1 X0
k= β = 0.00052
n I0
If we own 50’000 shares, we should thus short 50’000 · 0.00052 = 26
futures.
20
IV: Conditions for No Arbitrage
21
IV: Conditions for No Arbitrage
n may be huge!). Then the payoff Y of any derivative can also take on
only n values y1 , . . . , yn (some of which may coincide.)
We assume that the market is arbitrage free in the sense that the law
of one price prevails (Ch. I,) but also in the sense that
Proof
We associate with each contract (G, Y ) the vector (−G, y1 , . . . , yn ) in
Rn+1 . The set of such vectors constitutes a linear subspace V of Rn+1 .
We now prove that there is a vector q̄ which is orthogonal to V, all of
whose coordinates are positive.
Let K be the subset of Rn+1 : K = {ū = (u0 , . . . , un ) ∈ Rn+1 |
u0 + · · · + un = 1 and ui ≥ 0 for all i}. Obviously K and V have no
vector in common; indeed, it is easy to see that any such common vector
would represent an arbitrage. Now let q̄ be the vector of shortest Euclidean
length such that q̄ = ū−v̄ for some vectors ū and v̄ in K and V respectively.
The fact that such a vector exists needs to be proved, however, we will
skip that proof. We write q̄ = ū∗ − v̄ ∗ where ū∗ ∈ K and v̄ ∗ ∈ V.
Now note that for any t ∈ [0, 1] and any ū ∈ K, v̄ ∈ V , the vector tū+
(1−t)ū∗ ∈ K and tv̄+(1−t)v̄ ∗ ∈ V, hence |(tū+(1−t)ū∗ )−(tv̄+(1−t)v̄ ∗ )|
as a function of t on [0, 1] has a minimum at t = 0, by definition of ū∗ and
v̄ ∗ , i.e., |t(ū − v̄) + (1 − t)q̄|2 = t2 |ū − v̄|2 + 2t(1 − t)(ū − v̄) · q̄ + (1 − t)2 |q̄|2
has minimum at t = 0 which implies that the derivative w.r.t. t at t = 0
is ≥ 0. This gives (ū − v̄) · q̄ − |q̄|2 ≥ 0 or, equivalently, ū · q̄ − |q̄|2 ≥ v̄ · q̄
for all v̄ ∈ V and ū ∈ K. But since V is a linear space, it follows that we
must have v̄ · q̄ = 0 for all v̄ ∈ V (the inequality must hold if we replace
v̄ by λv̄ for any scalar λ.) It remains to prove that q̄ has strictly positive
entries. But we have ū · q̄ − |q̄|2 ≥ 0 for all ū ∈ K, in particular we can
take ū = (1, 0, . . . , 0) which shows that the first entry of q̄ is > 0 and so
on.
22
IV: Conditions for No Arbitrage
−q0 G + q1 y1 + · · · + qn yn = 0
for all contracts, and by scaling the qk :s we can arrange that q0 is equal
to one, hence
n
G= qk yk .
1
n
Since G[1] = 1 we have 1 qk = 1, so pX (xk ) = qk defines a probability
distribution for X such that for any payoff function f (X) it holds that
A Generalisation
In the discussion above we have considered only one underlying asset
X. But the arguments apply equally well to the situation when X repre-
sents a vector of underlying values X = (X1 , . . . , Xj ). Hence, if we have
a finite number of underlying assets X1 , . . . , Xj , all of which can take on
only a finite number of values, then, on an arbitrage free market, there
is a joint probability distribution for X = (X1 , . . . , Xj ) such that for any
derivative on X with payoff f (X) it holds that
Even though the assumption that the value of the underlying assets
can take on only a finite (but possibly huge) number of values is not
entirely unrealistic, it is inconvenient for the theoretical modelling. So we
do away with this assumption and henceforth assume the following No
Arbitrage Condition:
23
IV: Conditions for No Arbitrage
24
V: Pricing European Derivatives
Black’s Model
A European derivative on some underlying value X is a contract which
at a certain date, the date of maturity, pays the long holder of the contract
a certain function f (X) of the value X. A typical example is an option to
buy X for the strike price K. The strike price K is a deterministic value
written in the contract, and the payoff is f (X) = max(X − K, 0).
In order to compute the price of a European derivative on some un-
derlying asset value X, we assume that there is already a forward contract
for X on the market, or that we can compute the forward price of X, and
(t)
denote the current forward price for delivery at t by G0 [X]. Black’s
model assumes that the value X has a log-Normal forward probability
distribution (see comments below:)
√
X = Aeσ tz
where z ∈ N (0, 1). (1)
The current price of any derivative f (X) of X is thus, by (1) of Ch. IV,
(t) (t) 1 2
√
P0 f (X) = Zt E(t) [f (X)] = Zt E(t) f G0 [X] e− 2 σ t+σ t z
∞
Zt (t) 1 2
√ 1 2
=√ f G0 [X]e− 2 σ t+σ t x e− 2 x dx. (3)
2π −∞
This is Black’s pricing formula.
25
V: Pricing European Derivatives
is called the market price of risk. It measures the expected rate of return
on a forward contract on X per unit of σ.
26
V: Pricing European Derivatives
Comments
The specification (1) may need some clarification. The idea is this.
Let X be the value at some future time t of some asset, for example a share
of a stock, or a commodity like a barrel of oil. There are several different
(t)
random events which can influence this value. If G0 is the forward price
today, then tomorrow the forward price may change by a factor (1 + w1 ),
(t) (t)
i.e. tomorrow the forward price is G1 = G0 (1 + w1 ). Here w1 is a
random variable with rather small standard deviation.
The day after tomorrow, the forward price has changed by a new
(t) (t)
factor (1 + w2 ) such that the forward price then is G2 = G1 (1 + w2 ) =
(t)
G0 (1 + w1 )(1 + w2 ) and so on. Eventually, at time t when the forward
expires, the “forward” price is the same as the spot price, and we have
(t)
X = G0 (1 + w1 )(1 + w2 ) · · · (1 + wt ), i.e.,
(t)
ln X = ln G0 + ln(1 + w1 ) + ln(1 + w2 ) + · · · + ln(1 + wt )
27
V: Pricing European Derivatives
There is one problem with Black’s model applied to futures options. The
problem is that it is the forward price that should go into the formula, not
the futures price. However, the forward price equals the futures price if
the value of the asset is essentially independent of the interest rate. This
is proven in Ch. VII.
3. The same question as above, but now we assume the share pays a
dividend of 0.50 SEK in 3 months, all other assumptions are the
same. (4.115 SEK)
1
John C. Hull: Options, Futures, & Other Derivatives; Prentice Hall.
28
V: Pricing European Derivatives
6. Let S(t) be the spot price of a share at time t (year) which does not
pay dividends the following year. Determine the price of a contract
S(1)2
which after one year gives the owner . The risk-free rate of
S(0)
interest is 6% a year, and the share’s volatility is assumed to be 30%
for one year. (S(0) e0.15 )
7. Let S(t) be the spot price of a share at time t (year) which does not
pay dividends the following year. Determine the price of a contract
which after one year gives the owner $100 if S(1) > $50 and nothing if
S(1) ≤ $50. The current spot price is S0 = $45, the share’s volatility
is assumed to be 30% for one year and the risk-free rate of interest is
6% a year. ($35.94)
8. Show the put-call parity: If c and p are the prices of a call and a put
option with the same strike price K, then
c − p = Zt (G0 − K)
29
V: Pricing European Derivatives
Solutions
I suggest that you make a spread-
sheet in Excel (or some such)
where you can calculate call- and
put options with Black’s formula.
You input the parameters K, G0 ,
σ, time to maturity T and dis-
count factor ZT . You find the
relevant formulas in this chapter.
You can then use this spreadsheet
to solve most of the problems be-
low; the major problem is to find
the correct value for G0 to go into
the formula.
2. Since the share doesn’t pay any dividend and there is no convenience
(4)
yield, we get the forward price G0 from the relation (see example 1
(4) (4)
of Ch. II with d = 0) 45 SEK = e−0.09/3 G0 , i.e., G0 ≈ 46.37 SEK.
30
V: Pricing European Derivatives
The payoff of the derivative is then $100 · IA (S(1)), hence the value
of the contract is, with the usual notation
(1) 1 2
√
p = Z1 E 100 · IA (G0 e− 2 σ 1+σ 1 w )
= e−0.06 E 100 · IA (45 · e0.06−0.045+0.3w )
= e−0.06 100 · Pr 45 · e0.06−0.045+0.3w > 50
1 50
= e−0.06 100 · Pr w > ln − 0.05 ≈ $35.94
0.3 45
8. Let x+ denote max(x, 0). Then the price of a call option with maturity
(t)
t and strike price K on the underlying asset value X is P0 [(X −K)+ ],
(t)
and similarly, the price of the corresponding put option is P0 [(K −
+
X) ]. Hence
(t) (t)
c − p = P0 [(X − K)+ ] − P0 [(K − X)+ ]
(t)
= P0 [(X − K)+ − (K − X)+ ]
(t) (t) (t)
= P0 [X − K] = P0 [X] − P0 [K]
(t)
= Zt G0 [X] − K
31
VI: Yield and Duration
P0 (y)
D=−
P0 (y)
For coupon bearing bonds, and portfolios of bonds, the situation is more
complicated, since the value depends on more than one point on the zero
rate curve. The consequence is that the sensitivity of the value to changes
in the zero rate curve can not be exactly measured by a single number, and
furthermore that it is hard to come up with a model for pricing options
on such instruments.
The remedy is approximation. We will show that a coupon bearing
bond, or portfolio of bonds, can be approximated by a zero coupon bond
with a certain maturity, the duration of the bond or portfolio.
Note that for a zero coupon bond with time to maturity D and cor-
responding zero rate y0 , it holds that
1
D=− P (y0 ) (1)
P0 0
P0 (y) = PD e−yD (2)
−y(D−t)
Pt (y) = PD e (3)
where Ps (y) is the value at time s, if the zero coupon interest rate at that
time equals y, which may differ from y0 .
The idea of the concepts yield (y) and duration (D) is to define a
number D for a coupon paying bond, or a portfolio of bonds, such that
the relations (2) and (3) still hold approximately. Hence, let P0 be the
present value of a bond, or a portfolio of bonds:
n
P0 = Zti ci
i=1
where Zti is the price today of a zero coupon bond maturing at time ti
and ci is the payment received at time ti . Now define the yield y0 by the
relation
32
VI: Yield and Duration
n
P0 = e−y0 ti ci = P0 (y0 ) (4)
i=1
Here (2 ) follows from the definition of P̂D and (3 ) holds if the yield is y
at time t, as is seen as follows:
n
n
Pt (y) = e−y(ti −t) ci = eyt e−yti ci = eyt P0 (y) = P̂D e−y(D−t)
i=1 i=1
Summary
A bond, or a portfolio of bonds, has a yield y0 , defined by (4), and
a duration D, defined by (5). At any time t before any coupon or other
payments have been paid out, the value Pt (y) of the asset is approximately
equal to
33
VI: Yield and Duration
a portfolio of bonds,) but it’s growth rate goes down, and at the time
t = duration the two effects nearly net out. Similarly, of course, if the
rate of interest goes up. All this is captured in (3 ). Note, however, that
in the meantime coupon dividends have been paid out, and when these
are re-invested, the duration goes up, so in order to keep a certain date
as a target for the value of the portfolio, it has to be re-balanced.
Example
Consider a portfolio of bonds that gives the payment 1’000 after one
year, 1’000 after two years and 2’000 after three years. Assume that
Z1 = 0.945, Z2 = 0.890 and Z3 = 0.830. The present value of the portfolio
is then
so the duration is D = 7’716.5 years = 2.208 years. The value of the port-
3’495
folio at time t, if the yield at that time is y and t is less than one year,
can thus be approximated by
34
VI: Yield and Duration
n
(t)
Zt G0 = Zti ci
i=1
Pt = PF e−DF y
Summary
An interest rate security that after time t gives the deterministic
payments c1 , c2 , . . . , cn at times t1 < t2 < · · · < tn has a forward yield yF ,
defined by (6), and a forward duration DF defined by (7). The value of
the asset as of time t, 0 < t < t1 is then, to a first order approximation
(t)
Pt = PF e−DF y where PF = G0 eDF yF
35
VI: Yield and Duration
Example
In the previous example, assume that Z1.5 = 0.914. The forward price
G0 (1.5) of the portfolio for delivery in one and a half years is obtained
from the relation (note that the first payment of 1’000 has already been
paid out, and hence does not contribute to the forward value:)
(1.5) (1.5)
0.914 G0 = 0.890 · 1’000 + 0.830 · 2’000 ⇒ G0 = 2’790
(1.5)
2’790 = G0 = e−0.5yF 1’000 + e−1.5yF 2’000 ⇒ yF = 0.06258
Now, by (7)
(1.5)
DF G0 = 0.5 e−0.5yF 1’000 + 1.5 e−1.5yF 2’000 = 3’215.8
The minus sign is there for notational convenience later on; note that z
and −z have the same distribution. Hence, the present value Pt at time t
of the asset is
√ √
Pt = PF e−DF (α−σ t z)
= A eDF σ tz
where A = PF e−DF α .
Note that this is exactly the same specification as (1) in Ch. V with σ
replaced by DF σ.√We can thus use Black’s formula for European options
of Ch. V, where σ t in that formula stands for the product of the standard
deviation of the yield and the forward duration of the underlying asset.
36
VI: Yield and Duration
Portfolio Immunisation
Assume that we have portfolio P of bonds whose duration is DP . The
duration reflects its sensitivity to changes in the yield, so we might well
want to change that without making any further investments. One way
to do this is to add a bond futures to the portfolio (or possibly short such
a futures.)
Consider again the security discussed under “Forward Yield and For-
ward Duration” above. We will look at a futures contract on the value
Pt . Using the notation of Ch. I, the futures price is F0 = F0 [Pt ] =
F0 [PF e−DF y ], where the yield y at time t is random as seen from to-
day. Assume now that the present yield y0 of the portfolio P changes to
y0 + δy. The question arises what impact this has on y. If y = (current
yield) + (random variable independent of current yield), then y should
simply be replaced by y + δy. The new futures price should then be
F0
Dtot = DP + N DF
P
37
VI: Yield and Duration
1. Determine the
a) forward price (1’993.29 SEK)
b) forward yield (9.939%)
c) forward duration (1.861 years)
two years into the future for a bond that pays out 100 SEK in 2.5, 3
and 3.5 years and 2’100 SEK in 4 years. Zero-coupon interest rates
are at present 6%, 6.5%, 7%, 7.5% and 8% for the duration of 2, 2.5,
3, 3.5 and 4 years.
4. Calculate the value of a one-year put option on a ten year bond issued
today. Assume that the present price of the bond is 1’250 SEK, the
strike price of the option is 1’200 SEK, the one-year interest rate is
10% per year, the yield has a standard deviation of 0.013 in one year,
the forward duration of the bond as of the time of maturity for the
option is 6.00 years and the present value of the coupon payments
which will be paid out during the lifetime of the option is 133 SEK.
(20.90 SEK)
38
VI: Yield and Duration
Solutions
1. First we calculate the present price P0 of the bond:
We get the duration of the bond DP in the portfolio from the relation
39
VI: Yield and Duration
Next we calculate the present value of the cash flow from the
bond underlying the futures:
F0 50’577.84
Dtot = DP − DF = 1.915 − 2.791 years ≈ 0.514 years.
P 0 100’791.43
40
VI: Yield and Duration
4. Again we need the forward price of the underlying asset and the
appropriate value of σ. The value of G is easy; the present price
of the cash flow from the bond after the maturity of the option is
1’250 − 133 = 1’117 SEK, hence G = 1117e0,1 ≈ 1’234.48 SEK. The
relevant value of σ is σ = 0.013 · 6 = 0.0780. We get the option’s price
20.90 SEK from Black’s formula.
41
VII: Risk Adjusted Probability Distributions
(t) (t)
P0 [X] = P0 [N ] EN [ X
N] (1)
(t)
Here P0 [·] is the present price in any currency. In this way we have
defined a new probability distribution Pr N such that the probability of
any event A, whose outcome is known at time t, is defined by
(t)
P0 [IA · N ]
Pr N (A) = EN [IA ] = (t)
P0 [N ]
42
VII: Risk Adjusted Probability Distributions
The probability of any event is the ratio of the value of receiving one
unit of N conditional upon the happening of the event, and the value
of receiving one unit of N unconditionally.1
This seems to be a reasonable definition of the market’s assessment of
the probability of an event. It is called the risk adjusted probability with
respect to the numeraire N . However, this probability may depend on the
numeraire N employed. For example, if the event A is “the price of one
kilo of gold exceeds 10’000 EUR”, then the probability of A, defined by
employing the numeraire “one kilo of gold” is reasonably higher than if
the numeraire “one EUR” is employed.
If the numeraire N = one unit of some currency, i.e., the numeraire
is a zero coupon bond, then the probability distribution is the “usual”
forward probability distribution in that currency.
An Example
Assume that we want to price a call option on 1’000 barrels of crude
oil when the “strike price” is 3 kilos of gold and the date of maturity is
t. We model the value of one barrel of crude oil in terms of gold to be
log-Normal with respect to the forward-gold distribution:
X √
= Aeσ t w w ∈ N (0, 1)
N
where X is the value at time t of one barrel of crude oil, and N is the
value at time t of one kilo of gold. The forward-gold distribution of w is
N (0, 1). The parameter σ has to be somehow assessed by the writer of the
contract, e.g. by studying historical data (this parameter is assumed to
be the same for the “true” probability distribution as for the forward-oil
distribution.) The value of A can be computed using (1):
(t)
P0 [X] 1 2
(t) = EN [ X
N ] = Ae
2 σ t , i.e.,
P0 [N ]
(t)
P0 [X] − 1 σ 2 t
A= (t) e 2 .
P0 [N ]
Hence
X (t) √
P0 [X] − 1 σ 2 t+σ t w
= (t) e 2 w ∈ N (0, 1)
N P0 [N ]
The value-flow at time t for the long holder of the call option is (x+ denotes
max(x, 0))
1
Thomas Bayes’ (1702–1761) definition was: “The probability of any event is the
ratio of the value of which an expectation depending on the event ought to be computed,
and the value of the thing expected upon its happening.”
43
VII: Risk Adjusted Probability Distributions
(t) √
+
P [X] 1 2
(1’000X − 3N )+ = 1’000 0(t) e− 2 σ t+σ t w − 3 N
P0 [N ]
(t) (t)
It may be preferable to replace P0 [·] by Zt G0 [·] and factor out Zt to get
∞
Zt (t) 1 2
√
(t) + 1 2
c= √ 1’000G0 [X] e− 2 σ t+σ t x − 3 G0 [N ] e− 2 x dx
2π −∞
(t) (t)
The values G0 [X] and G0 [N ] can be observed on the futures market
(we assume that forward-price = futures price for these commodities.)
The expression for the option value is thus the same as Black’s “usual”
(t)
formula, when K is replaced by KG0 [N ], but the volatility σ is now the
volatility of the relative price of X compared to N .
Theorem 1
Assume that t < T and that the zero coupon bond price Z(t, T ) is
known at the present time 0. Then, for any random variable X whose
outcome is known at time t,
44
VII: Risk Adjusted Probability Distributions
Proof:
Consider the following strategy: Enter a forward contract on X ma-
turing at t (so G(t) = E(t) [X],) invest the net payment X − G(t) in zero
coupon bonds maturing at T . In this way we have constructed a contract
which pays Z(t, T )−1 (X − G(t) ) at time T and costs 0 today. Hence
Theorem 2
Assume that t < T and that the random interest rate R(t, T ) is inde-
pendent of the random variable X whose outcome is known at time t.
Then
Proof:
(t) (t) (T ) (1)
Zt E(t) [X] = Zt G0 [X] = P0 [X] = P0 [eR(t,T ) X]
(T )
= ZT G0 [eR(t,T ) X] = ZT E(T ) [eR(t,T ) X]
(2)
= ZT E(T ) [eR(t,T ) ] E(T ) [X]
(T )
= $−1 ZT G0 [$eR(t,T ) ] E(T ) [X]
(T ) (t) (3)
= $−1 P0 [$eR(t,T ) ] E(T ) [X] = $−1 P0 [$1] E(T ) [X]
= Zt E(T ) [X]
When E(t) [X] is independent of the date of maturity t (as long as the
outcome of X precedes time t,) we denote this expectation by a star:
E∗ [X].
45
VII: Risk Adjusted Probability Distributions
Theorem 3
Assume the random interest rate R(0, t) is independent of the random
variable X whose outcome is known at time t. Then the forward price
(t) (t)
and the futures price coincide: G0 [X] = F0 [X].
1. Show that
(t) (t) X
G0 [X] = G0 [N ] EN N
4. Prove theorem 3.
Solutions
1. This follows immediately from (1) if we divide both sides of the equal-
ity with Zt (see Ch. I.)
2. Let g[X] be the forward price of X in terms of gold. For the long
holder of the contract, the cash flow at maturity is then X − g[X]N .
The present value of this cash flow is zero (= value of contract.) Hence
(t) (t) (t)
0 = P0 X − g[X]N = P0 [X] − g[X]P0 [N ].
Hence
(t) (t)
P [X] G0 [X]
g[X] = 0(t) = (t)
P0 [N ] G0 [N ]
(the last equality follows as in exercise 1.)
46
VII: Risk Adjusted Probability Distributions
X √
= Aeσ t w (2)
Y
where w ∈ N (0, 1) with respect to√the forward-Y probability distri-
bution. The standard deviation σ t of ln X − ln Y is assumed to
be the same as under the “true” probability distribution (if there is
any such thing) and has to be somehow assessed or estimated. Taking
EY of (2) gives, as in the example above,
(t) (t)
1 2 P0 [X] G0 [X]
A e2σ t
= (t) = (t)
P0 [Y ] G0 [Y ]
so
X (t) √
G0 [X] − 1 σ 2 t+σ t w
= (t) e 2
Y G0 [Y ]
4. The fact that the interest rate R(0, t) is independent of X (under the
t-forward distribution) implies that
Hence
(t) (t) (t)
F0 [X] = P0 [XeR(0,t) ] = Zt G0 [XeR(0,t) ]
= Zt E(t) [XeR(0,t) ] = Zt E(t) [X] E(t) [eR(0,t) ]
(t) (t) (t)
= G0 [X] P0 [$eR(0,t) ]$−1 = G0 [X]
47
VIII: Conditional Expectations and Martingales
Properties
– If X is stochastically independent of any random variable observed
at time t, then the conditional expectation equals the unconditional
expectation: Et [X] = E[X].
The first property is the special case of the second when Y ≡ 1. In order
to see the second property, recall that if two random variables are inde-
pendent, then the expected value of their product is the product of their
expected values. Hence, if A is any event whose outcome is observed at
time t, then
48
VIII: Conditional Expectations and Martingales
E[IA Y X] = E[IA Y ]E[X] = E IA Y E[X] .
Hence Et [Y X] = Y E[X].
Assume now that the time interval from 0 to T is discretised by the points
0 = t0 < t1 , . . . < tn = T , and that we have a sequence of random variables
Xj , j = 0, . . . , n, where the outcome of Xj is observed at time tj . Such a
sequence is a stochastic process. A stochastic process is a martingale if it
is true that
Martingale Prices
Consider a contract that gives the random payoff X at time T . The
(T )
forward price Gt , 0 ≤ t < T , of this contract at time t is a random
variable whose outcome is determined at time t.
Theorem
(T )
For any t, 0 ≤ t ≤ T , the forward price Gt equals the conditional
expectation
(T ) (T )
Gt = Et [X]
Hence, the forward process has the martingale property (3) w.r.t. the
forward distribution, i.e.,
(T ) (T ) (T )
Gtj = Etj [Gtj+1 ]
49
VIII: Conditional Expectations and Martingales
Proof
Consider the following strategy: Let A be any event whose outcome
occurs at time t. Wait until time t and then, if A has occurred, enter a
(T )
forward contract on X maturing at T with forward price Gt , but if A
has not occurred, do nothing. We have thus constructed a contract which
(T )
gives payoff IA (X − Gt ) at time T which costs 0 today. Hence
(T ) (T )
ZT E(T ) [IA (X − Gt )] = 0 i.e., E(T ) [IA X] = E(T ) [IA Gt ]
Since this is true for any event A whose outcome is known at time t, we
have by the definition of conditional expectation
(T ) (T )
Gt = Et [X]
Q.E.D.
Remark 1
From this we deduce (c.f. Ch. I) that the present price pt at time t
of the contract yielding X at time T is
(T )
pt = Z(t, T )Et [X]
Remark 2
The Martingale property holds for distributions referring to any nu-
meraire. By the same token as above, it follows that
(T )
Pt [X] X
(T ) = EN
t [N ]
Pt [N ]
50
IX: Asset Price Dynamics and Binomial Trees
In Black’s model, there are only two relevant points in time: the time
when the contract is written, and the time when it matures. In many
cases we have to consider also what happens in between these two points
in time.
First some notation. Time is discrete, t = t0 < . . . < tn = T , and
we let for notational simplicity t0 = 0 and all time spells tk − tk−1 = Δt
be equally long; hence tk = kΔt. Let X be some random value whose
outcome becomes known at time T ; we are interested in evaluating the
prices of derivatives of X.
Black-Scholes Dynamics
The Black-Scholes assumption about the random behaviour of X is
that
n
X = Aeσ Σ1 wj (1)
Binomial Approximation
In almost all cases when the whole price path has to be taken into
account, it is necessary to use some numerical procedure to calculate the
price of the derivative. A useful numerical procedure is to approximate
√ In (1), we replace wj by a binary
the specification (1) by a binomial tree.
variable bj which takes the value − Δt with probability
√ (under the for-
ward probability distribution) 0.5 and the value + Δt with probability
0.5. For large values of n, and small Δt, this is a good approximation—in
fact, as n → ∞ and Δt → 0 the price of a derivative calculated from the
binomial tree will converge to the theoretical value it would have from the
specification (1).
51
IX: Asset Price Dynamics and Binomial Trees
Combining this relation with the same for Gk+1 , we get the relation
eσ bk+1 √
Gk+1 = Gk √ = Gk 1 ± ε) where ε = tanh(σ Δt)
cosh(σ Δt)
and the plus and minus sign occur with probability (under the forward
probability distribution) 0.5 each.
G0 G0 u G 0 u 2 G0 u 3 G0 u 4 ··· G0 un
G0 d G0 du G0 du G0 du3
2
··· G0 dun−1
G 0 d2 G0 d 2 u G 0 d2 u 2 ··· G0 d2 un−2
G0 d 3 G0 d3 u ··· G0 d3 un−3
.. ..
. ··· .
G 0 dn
52
IX: Asset Price Dynamics and Binomial Trees
where Δt is the time spell between columns, and σ is the volatility of the
forward price.
Here the prices in the rightmost column are known, since they represent
the value of the option at maturity when the underlying futures price
is that in the corresponding position of the tree of futures prices. Then
option prices are calculated backwards:
p̂j,k = e−rΔt 0.5 pj+1,k + 0.5 pj+1,k+1 (1)
pj,k = max(p̂j,k , value if exercised) (2)
Indeed, p̂j,k = e−rΔt E∗ [pj+1,• ], which is just another way of writing (1),
is the value if the option is not exercised (see (4) of Ch. VIII,) and the
highest value of this and the value if exercised is the present value. In
order to calculate the value if exercised, we employ the original tree of
futures prices, of course. The end result p is the price of the option today.
53
IX: Asset Price Dynamics and Binomial Trees
The price of the option is then calculated in the same way as above.
It is now easy to build the tree of stock prices: First compute G0 from (4),
then construct the tree for Gj , and from these prices, compute Sj from
(3).
54
IX: Asset Price Dynamics and Binomial Trees
For times j before i we can compute the stock price in accordance with
example 2 of Ch. II:
Sj (1 − d) = e−r(n−j)Δt Gj
It is now easy to build the tree of stock prices: First compute G0 from (7),
then construct the tree for Gj , and from these prices, compute Sj from
(5) and (6).
Comments
There is more than one way to specify the binomial model. The one
given here is slightly different from the original model by Cox, Ross and
Rubenstein. I find the specification given here more appealing, and also
easier to generalise to different situations (see for example Ch. XII, where
we discuss an interest rate model.) It is important not to mix up different
specifications!
55
IX: Asset Price Dynamics and Binomial Trees
One can construct binomial trees for various underlying assets: in-
dices, currencies, futures, and shares paying various forms of dividends.
In order to make these appear less ad-hoc, I do this in two steps: First
we construct a tree over the forward prices of the underlying asset. These
trees are constructed the same way independent of underlying asset. Then
we construct a tree over the relevant spot prices (unless the underlying is
a futures contract,) where we use the techniques from Ch. II. Finally, we
construct the tree over the derivative price.
Many texts construct the binomial tree directly for the spot prices un-
derlying the derivative, without first specifying one for the corresponding
forward prices. This is a somewhat risky approach: the tree of spot prices
must be consistent with the fact that the corresponding forward prices
have the Martingale property. In order to ensure that this happens, I pre-
fer to first make a tree over the forward prices such that the Martingale
property is satisfied. Now I know that there are no arbitrage possibilities
lurking in the model; thereafter I derive the relevant spot prices, and I
can feel safe that the model is sound.
One can also use trinomial trees, where the N (0, Δt) variable is approxi-
mated by the “trinary” variable
⎧√
⎪
⎨ 3Δt with probability 16
u= 0 with probability 23
⎩ √
⎪
− 3Δt with probability 16
This variable shares the same five first moments with the N (0, Δt) vari-
able, so the approximation is better in a corresponding trinomial tree
compared to a binomial tree.
Trinomial trees are also used in many interest rate models (Hull and
White.)
56
IX: Asset Price Dynamics and Binomial Trees
57
IX: Asset Price Dynamics and Binomial Trees
6. A two month American put option on an index of shares has the strike
price 480. The present value of the index is 484, the risk-free rate of
interest is 10% per year, the dividend of the index is 3% per year and
the volatility is 25% over a year. Determine the value of the option
by using a binomial tree with time interval of half a month. (15.2336)
7. The spot price of copper is $0.60 per pound. Assume that, at present,
the futures price of copper is
maturity futures price
3 months 0.59
6 months 0.57
9 months 0.54
12 months 0.50
The volatility of the price of copper is 40% during one year and the
risk-free interest rate is 6% a year. Use a suitable binomial tree with
time interval three months to estimate the price of an American call
option on copper with strike price $0.60 and maturity in one year.
($0.063167.)
Solutions
1.
Here are three binomial trees (the display shows decimal commas,
which is Swedish standard) cut from an Excel spreadsheet. Note
58
IX: Asset Price Dynamics and Binomial Trees
that time goes from bottom to top; I think this is the easiest way to
arrange binomial trees in a spreadsheet.
The first tree displays the futures (=forward) prices of coffee.
Week 0 (today) it is $70, week 1 it is either $71.3998 or $68.6002, each
with probability 0.5 under the forward probability measure. Week
two, the futures price is $72.8276 or $69.9720 if the futures price were
$71.3998 week 1, and $69.9720 or $67.2284 if the futures price were
$68.6002 week 1; and so on. This tree is generated in accordance with
the first tree (forward prices) appearing in Ch. IX.
The second tree displays the value of the European put option.
It is generated from top to bottom: week 4 the option is worth 0 if
the futures price is more than or equal to $72, otherwise it is worth
$72 − futures price (where the futures price equals the spot price of
coffee, since the futures contract matures week 4.) Hence we get
the first row (week 4) easily from the previous tree’s week 4. Next we
compute the week 3 option prices, employing (4) of Ch. VIII: the price
week 3 is the discounted value of the expected value week 4 under the
forward probability measure. Hence, each price week 3 is e−0.001 times
(0.5·p1 +0.5·p2 ) where p1 and p2 are the prices above and right-above
respectively, For example, 3.4238 = e−0.001 (0.5 · 2.0560 + 0.5 · 4.7985).
Going down through the tree, we finally find the price week 0 to
be $2.4255.
The third tree displays the values of the corresponding American
put option. It is obtained as in the tree on the European option,
except that the price is replaced by the exercise value if that is higher.
These numbers are in boldface. For instance, the value $3.4272 is the
value if the option is exercised: 3.4272 = 72 − 68.5728 which is more
than e−0.001 (0.5 · 2.0560 + 0.5 · 4.7985).
Going down through this tree, we end up with the price $2.4287
for the American option.
59
IX: Asset Price Dynamics and Binomial Trees
2a.
First we compute the forward price of the share according to Ch. II:
G0 = S0 e0.004 = $13.3332. Starting out with this value, we construct
the first tree of forward prices. Now we know that the forward prices
of the underlying share indeed have the Martingale property under
the forward probability measure.
Next we construct the tree of spot prices of the share: for exam-
ple, week 2 the spot price S2 = G2 e−2·0.001 , so we get the spot prices
by just discounting the values in the forward price tree. Of course,
week 0 we get the initial spot price of the share.
Now we construct the tree of option prices. We start with the
prices week 4. Next week three: these prices are computed as the
mean of the two prices above and right-above, discounted by the one
period interest rate. However, we must check that it isn’t profitable
to exercise the option. This happens in a few places; these are shown
in boldface. Finally, the option price is found to be $0.464143.
60
IX: Asset Price Dynamics and Binomial Trees
2b.
61
IX: Asset Price Dynamics and Binomial Trees
3.
62
IX: Asset Price Dynamics and Binomial Trees
4.
63
IX: Asset Price Dynamics and Binomial Trees
6.
7.
64
IX: Asset Price Dynamics and Binomial Trees
65
X: Random Interest Rates: The Futures Distribution
Theorem 1
The futures distribution as defined above is independent of the date of
maturity T in the following sense: if the outcome of X is known at
time tn then E (n) [X] = E
(m) [X] if m > n, where we denote by E (k)
the futures distribution for contracts maturing at time tk .
Theorem 2
The futures prices {Fj } have the martingale property w.r.t. the futures
distribution: Fj = E t [Fk ] for j < k. In particular, the futures price
j
F0 is the expected value w.r.t. the futures distribution of X, the spot
price at delivery.
66
X: Random Interest Rates: The Futures Distribution
Proof of Theorem 1
Let P0 be the present price of the value XeR(0,n) to be delivered at
time tn . We can convert this contract to one where instead the value
XeR(0,m) is delivered at time tm by simply depositing the payoff XeR(0,n)
in the money market account up to time tm . The present price is of course
the same. Using the theorem in Ch. I, we have the following equalities
from the two contracts:
(tn ) (tn )
P0 = P 0 [XeR(0,n) ] = F0 [X] and
(t ) (t )
P0 = P0 m [XeR(0,m) ] = F0 m [X]
(n) [X] = E
E (m) [X]
Q.E.D.
Proof of theorem 2
Consider the following strategy: Let A be any event whose outcome is
known at time tj < tk . At time tk−1 enter a long position of eR(0,k) futures
contracts if A has occurred, otherwise do nothing. At time tk collect
IA (Fk − Fk−1 )eR(0,k) (which may be negative) and close the contract.
This gives the payment IA (Fk − Fk−1 )eR(0,k) at time tk and no other cash
flow. The present price is zero, hence
0 = P0
(tk ) A (Fk − Fk−1 )]
[IA (Fk − Fk−1 )eR(0,k) ] = F0 [IA (Fk − Fk−1 )] = E[I
Hence,
A Fk ] = E[I
E[I A Fk−1 ]
A Fj ] = E[I
E[I A Fk ]
Since A can be any event whose outcome is known at time tj , this means
that (Ch. VIII)
t [Fk ]
Fj = E for j < k.
j
67
X: Random Interest Rates: The Futures Distribution
Theorem 3
The following formula for the present price pj as of time tj for the
asset Y to be delivered at time tk > tj holds:
t [Y e−R(j,k) ].
pj = E j
Proof
With the notation of Ch. I,
(k) t [Y e−R(j,k) ]
pj = Fj [Y e−R(j,k) ] = E j
Q.E.D.
68
XI: A Model of the Short Interest Rate: Ho-Lee
where θk are some numbers, the factor σ is the volatility of the short in-
terest rate. The zj :s are independent N(0,1)-variables and the outcome of
each zj occurs at time tj . This is under the futures probability distribu-
tion.
It remains to compute the exact value of θk , but first a notational
simplification: multiplying by Δt yields
3
rk Δt = θk Δt + σΔt 2 (z1 + · · · + zk−1 )
Let $Ztk be the price of a zero coupon bond maturing at tk with face value
$1. Then, according to Theorem 3, Ch. X,
e−r1 −···−rk ]
Ztk = E[1
In order for the model (1) to correctly represent the current term structure,
we must thus have
69
XI: A Model of the Short Interest Rate: Ho-Lee
Z
σ2
t−1
rt = ln + (t − 1)2 + σ(z1 + · · · + zt−1 ) (2)
Zt 2
t [e−rt+1 −···−rT ]
Z(t, T ) = E
ZT − σ2 (t2 +···+(T −1)2 )−σ(T −t)(z1 +···+zt )
= Et [e 2
Zt
× e−σ((T −t−1)zt+1 +(T −t−2)zt+2 +···+1 zT −1 ) ]
ZT − σ2 (t2 +···+(T −1)2 ) −σ(T −t)(z1 +···+zt )
= e 2 e
Zt
×E t [e−σ((T −t−1)zt+1 +(T −t−2)zt+2 +···+1 zT −1 ) ]
ZT −σ(T −t)(z1 +···+zt ) − σ2 (t2 +···+(T −1)2 ) σ2 ((T −t−1)2 +···+12 )
= e e 2 e2
Zt
σ2 2 σ2
− (t + · · · + (T − 1)2 ) + ((T − t − 1)2 + · · · + 12 )
2 2
T −t−1
σ2
=− (t + j)2 − j 2
2 j=0
T −t−1
σ2
=− t (2j + t)
2 j=0
σ2
=− (T − 1)(T − t)t
2
Hence,
ZT − σ2 (T −1)(T −t)t −σ(T −t)(z1 +···+zt )
Z(t, T ) = e 2 e (3)
Zt
70
XI: A Model of the Short Interest Rate: Ho-Lee
(t) ZT
G0 = $
Zt
regardless of the interest rate model. We leave the proof to the reader.
The futures price F0 is
ZT e− σ2
2
F0 = E[$Z(t, T )] = E[$ (T −1)(T −t)t −σ(T −t)(z1 +···+zt )
e ]
Zt
2
(t)
= G0 e − σ2 (T −1)(T −t)t −σ(T −t)(z1 +···+zt ) ]
E[e
σ2 σ2 σ2
(t) (T −t)2 t (t)
= G0 e− 2 (T −1)(T −t)t
e 2 = G0 e− 2 (T −t)(t−1)t
We see that in contrast to the situation when interest rates are determinis-
tic, the forward price G0 and the futures price F0 differ. Since the price of
the underlying asset (the bond) is negatively correlated with the interest
σ2
(T −t)t2
rate, the futures price is lower; in the Ho-Lee model by factor e− 2 .
71
XI: A Model of the Short Interest Rate: Ho-Lee
1 2
= $C3 . . . e− 2 Σ(xj +σ(t−j)) dx1 . . . dxt
D
Here the equalities (1) follow from the theorem in Ch. I. Hence, under
random variables zj (j ≤ t) are
the t-forward probability distribution, the
independent and zj ∈ N − σ(t − j), 1 . In particular, we have from (3)
that under the E(t) -distribution, with wj independent ∈ N (0, 1) :
But $Z T
Zt = G0 , the forward price of the underlying bond, hence, (also
with no normalisation of time,)
p = Zt E(t) f $Z(t, T )
∞
Zt σ2 2
√ 1 2
=√ f G0 e− 2 (T −t) t eσ(T −t) t x e− 2 x dx
2π −∞
Note that the same formula may be obtained by Black’s model for bond
options (Ch. VI.)
72
XII: Ho-Lee’s Binomial Interest Rate Model
where {bk } are binomial random variables which takes the values ±1, each
with probability 0.5 under the futures distribution; they are thus assumed
to be statistically independent. The outcome of the variable bk occurs at
3
time tk−1 . The volatility parameter σ is proportional to Δt 2 just as in
Ch. XI.
We will choose the terms θj :s such that the model is consistent with
the current term structure—it will be close to, but not exactly, the same
as in the Normal distribution case.
Let us compute the price Z4 at time t0 of a zero coupon bond maturing
at time t4 with face value 1:
· e−r1 −r2 −r3 −r4 ] = E[e
Z4 = E[1 −θ1 −θ2 −θ3 −θ4 −σ(3b2 +2b3 +b4 ) ]
= e−θ1 −θ2 −θ3 −θ4 E[e −2σb3 ] E[e
−3σb2 ] E[e −σb4 ]
= e−θ1 −θ2 −θ3 −θ4 cosh(3σ) cosh(2σ) cosh(σ)
And, by the same token, in general
Zk = e−θ1 −···−θk cosh (k − 1)σ · . . . · cosh(σ)
Since this must hold for all k, we must have
Z
k−1
θk = ln + ln cosh (k − 1)σ
Zk
and the dynamics under the futures distribution is thus described by
Z
k
k−1
rk = ln + ln cosh (k − 1)σ + σ bj k = 1, . . . , n
Zk 2
1
bj = 1 with probability 2
where 1
bj = −1 with probability 2
Once we have the parameters of the model, we can price interest rate
derivatives in the binomial tree. We show the procedure by an example,
73
XII: Ho-Lee’s Binomial Interest Rate Model
t0 t1 t2 t3
0.06 0.071 0.082 0.093
0.051 0.062 0.073
0.042 0.053
0.033
The interest rate from one period to the next is obtained by going either
one step to the right on the same line, or step to the right to the line
below; each with a (risk adjusted) probability of 0.5. We can compute the
value at t2 of the bond: since its value at t4 is 100, the value at t3 and
t2 is obtained recursively backwards by use of the formula of theorem 3
in Ch. X, i.e., we take the average of the values in the next period and
discount by the interest rate given in the above tree:
t2 t3
84.794 91.119
88.254 92.960
91.856 94.838
96.754
The value of the option can now also be obtained by backward recursion:
t0 t1 t2
2.309 1.050 0
3.853 2.254
5.856
The value of the option is thus 2.309. It is also easy to price other more
exotic derivatives in this binomial tree model.
Comments
The model presented in this and the previous chapter is Ho-Lee’s
interest rate model. It was originally presented as a binomial model. It is
the simplest interest rate model, and it can be calibrated so as to fit the
current term structure perfectly. There are two features with this model
that many people feel unhappy with: (1) the short rate will eventually
far into the future take negative values with a probability that approaches
one half, and (2) there is no “mean reversion” which means that all shocks
to the short rate are permanent; they don’t wear off with time.
74
XII: Ho-Lee’s Binomial Interest Rate Model
where 0 < β ≤ 1 is the mean-reverting factor. As you can see, the shock
z1 , for instance, wear off exponentially with time, since β (k−2)Δt → 0 as
k → ∞ if β < 1. The Ho-Lee model is the special case β = 1 (compare
with the first expression in Ch. XI.) In Hull-White’s model, the standard
σ
deviation of the short rate is always bounded by √ , so if β is small,
−2 ln β
the probability of negative interest rates can stay small.
This model can not be represented by a (recombining) binomial tree
(if β = 1,) but it can be implemented as a trinomial tree; see the book
by Hull1 . The difference between Ho-Lee’s and the Hull-White model is
solely on technicalities concerning the construction of the tree of interest
rates; the conceptual ideas are the same.
period 0 1 2 3
3.0 3.3 3.5 3.8
2.9 3.1 3.4
2.7 3.0
2.6
The interest from period 0 to period 1 is hence 3.0%, and so on. An
interest rate security is in period 3 worth 400 r where r is the interest
(if r = 3.8% the value is $1’520, if r = 2.6 the value is $1’040 etc.)
a) Determine today’s futures price of the security to two decimal
places. ($1’280.00)
b) Determine today’s forward price of the security to two decimal
places. ($1’279.52)
1
John C. Hull: Options, Futures, & Other Derivatives; Prentice Hall.
75
XII: Ho-Lee’s Binomial Interest Rate Model
The numbers are the interest rate in percent per 6 months, i.e., per
time step in the tree. ($9’896.85)
76
XII: Ho-Lee’s Binomial Interest Rate Model
a) Use a binomial tree with time interval of one year based on Ho-
Lee’s model. (74.5689)
Solutions
1. Here is first the tree of interest rates:
77
XII: Ho-Lee’s Binomial Interest Rate Model
The
value of the bond is calculated backwards in time:
0.5(number above + number above-right) × discount where we get
the discounting rate from the above tree:
78
XII: Ho-Lee’s Binomial Interest Rate Model
b. This is an ordinary call option where the forward price of the under-
lying asset is 8’311.04 SEK and the volatility is 0.015 · 2 = 0.03. The
factor 2 is the forward duration (time to maturity of the underlying
bond when the option is exercised.)
3. Here is the interest rate tree, now with time going from bottom to
top:
79