One of the possible improvements of the classical Black–Scholes option pricing model is to incorporate the stochastic nature of the short rate dynamics in option valuation. In this paper, we present the numerical scheme, based on the discontinuous Galerkin method, for European option pricing when the short rate follows the Merton model. The pricing function satisfies a partial differential equation with two underlying variables — stock price and short rate value. With a localization to a bounded spatial domain, including setting the proper boundary conditions, the governing equation is discretized by the discontinuous Galerkin method over a finite element grid and Crank-Nicolson time integration is applied, consequently. As a result the numerical scheme is represented by a sequence of linear algebraic systems with sparse matrices. Moreover, the numerical simulations reflect the capability of the scheme presented.

1.
F.
Black
and
M.
Scholes
, “
The pricing of options and corporate liabilities
,”
J. Polit. Econ.
81
,
637
659
(
1973
).
2.
R. C.
Merton
, “
Theory of rational option pricing
,”
Bell J. Econ.
4
,
141
183
(
1973
).
3.
J.
Hozman
and
T.
Tichý
, “A discontinuous Galerkin method for numerical pricing of European options under Heston stochastic volatility,” in
AMEE ’16
, AIP Conf. Proc No.
1789
, edited by
V.
Pasheva
,
N.
Popivanov
, and
G.
Venkov
(
American Institute of Physics
,
Melville, NY
,
2016
) p.
030003
.
4.
J.
Hozman
and
T.
Tichý
, “A DG approach to the numerical solution of the Stein–Stein stochastic volatility option pricing model,” in
AMEE ’17, AIP Conf. Proc No
.
1910
, edited by
V.
Pasheva
,
N.
Popivanov
, and
G.
Venkov
(
American Institute of Physics
,
Melville, NY
,
2017
) p.
030006
.
5.
J.
Hozman
and
T.
Tichý
, “Numerical pricing of options under the exponential Ornstein–Uhlenbeck stochastic volatility model based on a DG technique,” in
AMEE ’18, AIP Conf. Proc No
.
2048
, edited by
V.
Pasheva
,
N.
Popivanov
, and
G.
Venkov
(
American Institute of Physics
,
Melville, NY
,
2018
) p.
030012
.
6.
B.
Rivière
, Discontinuous Galerkin Methods for Solving Elliptic and Parabolic Equations: Theory and Implementation (
SIAM
,
Philadelphia
,
2008
).
7.
V.
Dolejší
and
M.
Feistauer
,
Discontinuous Galerkin Method. Analysis and Applications to Compressible Flow
(
Springer
,
Cham
,
2015
).
8.
R. C.
Merton
, “
A dynamic general equilibrium model of the asset market and its application to the pricing of the capital structure of the firm
,” (
1970
), (working paper no. 497).
9.
J. J.
Kung
and
L.-S.
Lee
, “
Option pricing under the merton model of the short rate. math comput simulation
,”
Math. Comput. Simulat.
80
,
378
386
(
2009
).
10.
Z.
Guo
, “
Option pricing under the Merton model of the short rate in subdiffusive Brownian motion regime
,”
J. Stat. Comput. Sim.
87
,
519
529
(
2017
).
11.
O. A.
Oleinik
and
E. V.
Radkevič
,
Second Order Equations with Nonnegative Characteristic Form
(
Springer–Verlag
,
Boston
,
1973
).
12.
H. R.
Stoll
, “
The relationship between put and call option prices
,”
J. Financ.
24
,
802
824
(
1969
).
13.
J.
Hozman
and
T.
Tichý
, “
DG framework for pricing European options under one-factor stochastic volatility models
,”
J. Comput. Appl. Math.
344
,
585
600
(
2018
).
14.
Z.
Cui
and
D.
McLeish
, “
Comment on “Option pricing under the Merton model of the short rate
" by Kung and Lee [Math. Comput. Simul. 80 (2009) 378–386],”
Math. Comput. Simulat.
81
,
1
4
(
2010
).
15.
F.
Hecht
, “
New development in freefem+
,”
J. Numer. Math.
20
,
251
265
(
2012
).
This content is only available via PDF.
You do not currently have access to this content.