DG Method for Pricing European Options under Merton Jump-Diffusion Model

Under real market conditions, there exist many cases when it is inevitable to adopt numerical approximations of option prices due to non-existence of analytical formulae. Obviously, any numerical technique should be tested for the cases when the analytical solution is well known. The paper is devoted to the discontinuous Galerkin method applied to European option pricing under the Merton jump-diffusion model, when the evolution of the asset prices is driven by a Lévy process with finite activity. The valuation of options under such a model with lognormally distributed jumps requires solving a parabolic partial integro-differential equation which involves both the integrals and the derivatives of the unknown pricing function. The integral term related to jumps leads to new theoretical and numerical issues regarding the solving of the pricing equation in comparison with the standard approach for the Black-Scholes equation. Here we adopt the idea of the relatively modern technique that the integral terms in Merton-type models can be viewed as solutions of proper differential equations, which can be accurately solved in a simple way. For practical purposes of numerical pricing of options in such models we propose a two-stage implicit-explicit scheme arising from the discontinuous piecewise polynomial approximation, i.e., the discontinuous Galerkin method. This solution procedure is accompanied with theoretical results and discussed within the numerical results on reference benchmarks.