58
CHAPTER
2.
NUMERICAL
INTEGRATION.
BONDS.
2.3.
SOLUTIONS
TO
SUPPLEJYIENTAL
EXERCISES
59
Note
that
, unlike classical
duration
and
convexity, which can only
be
computed for individual bonds,dollar duration
and
dollar convexity can
be
estimated for any bond portfolio, assuming all bond yields change
by
the
same amount. In particular, for a bond
with
value B,
duration
D ,
and
convexity C,
the
dollar
duration
and
the
dollar convexity can
be
computed as
D$
=
BD
and
G$
=
BG.
You invest
$1
million in a
bond
with
duration 3.2
and
convexity
16
and
$2.5 million
in
a bond with
duration
4
and
convexity
24.
(i)
What
are
the
dollar
d
旧
ation
and
dollar convexity of your portfolio?
(ii)
If
the
yield goes up
by
ten
basis
points
,五叫
new
approximate values
for each of
the
bonds.
What
is
the
new value of
the
portfolio?
(iii) You can buy or sell two
other
bo
时
s
,
one
with
d
旧时
ion
1.
6
and
convexity
12
and
another one
with
duration
3.2
and
convexity
20.
What
positions could you take in these bonds
to
immunize your portfolio
(i.
e.
,
to
obtain
a portfolio
with
zero dollar
duration
and
dollar
co
盯
exity)?
2.3
Solutions
to
Supplemental
Exercises
Problem
1: Assume
that
the
continuously compounded instantaneous
rate
curve r(t) is given by
γ
(t)
二
0.05
xp(
一
(1
+
t)2)
.
Use Simpson's Rule
to
compute
the
1-year
and
2-year discount factors
with
six decimal digits accuracy,
and
compute
the
3-year discount factor
with
eight decimal digits accuracy.
(ii)
Find
the
value of a three year yearly coupon
bond
with
coupon
rate
5%
(and face value 100).
Solution: (i) Recall
that
the
discount factor corresponding
to
time
t is
叫-fa'州市
Using Simpson's Rule,
we
obtain
that
the
1-year
, 2-year,
and
3-year
discount
factors are
disc(l)
= 0.956595; disc(2) = 0.910128; disc(3) = 0.86574100.
(ii)
The
value of
the
three year yearly coupon
bond
is
B = 5 disc(l) +5 disc(2) +105 disc(3) =
100.236424.
口
Problem
2: Consider a six months
plai
丑
vanilla
European
put
option with
strike
50
on a lognormally distributed underlying asset paying dividends
co
坠
tinuouslyat
2%.
Assume
that
interest rates are constant
at
4%.
Use
risk-neutral
valuation
to
write
the
value of
the
put
as
an
integral
over
a
且
nite
interval.
Find
the
value of
the
put
option
with
six decimal digits
accuracy using
the
Midpoint Rule
and
using Simpson's Rule. Also,compute
the
Black-Scholes value PBS of
the
put
and
report
the
approximation errors
of
the
numerical integration approximations
at
each step.
Solution:
If
the
underlying asset follows a lognormal distribution,
the
value
S(T)
of
the
underlying asset
at
maturity
is a lognormal variable given by
/σ2\
广-
In(S(T))
=
In(S(O))
+
(γ
-
q - v
2
) T
+
σ
飞
ITZ
,
where
σis
the
volatility of
the
underlying asset.
Then
,
the
probability density
function
h(
ν)
of
S(T)
is
/
(l叼一附
(0)
)一
(r
-
q
一手)
T)2
i
h(y)
=
~呵!一\时
LL
" (2.14)
if
y > 0,
and
h(y)
= 0 if
y
三
O.
Using
risk-neutral
valuation,
we
find
that
the
value of
the
put
is given
bγ
P -
e-rTE
R
川
max(K
-
S(T)
,
0)]
I'
K
=
e-
rT
/
(K
-
y)h(
ν
)
dy, (2.15)
JQ
where
h(y)
is given by (2.14).
The
Black-Scholes value of
the
put
is PBS = 4.863603. To compute a
numerical appr()ximation of
the
integral (2.15),
we
start
with a partition
of
the
interval
[0
,K] into 4 intervals,
and
double
the
numbers of -intervals
up
to
8192 intervals. We report
the
Midpoint Rule and Simpson's Rule
appro
均
nations
to
(2.15) and
the
correspondi
吨
approximation
errors
to
the
Black-Scholes value PBS in
the
table below:
We
丑
rst
note
that
the
approximation error does
not
go below 6 .
10-
6
.
This is due
to
the
fact
that
the
Black-Scholes value of
the
put
,which is given
by
PBS
Ke-r(T-t)
N(
-d
2
) - Se-q(T-t)
N(
-d
1
)
,