Efficient numerical methods based on integral transforms to solve option pricing problems

Abstract

In this thesis, we design and implement a class of numerical methods (based on integral transforms) to solve PDEs for pricing a variety of financial derivatives. Our approach is based on spectral discretization of the spatial (asset) derivatives and the use of inverse Laplace transforms to solve the resulting problem in time. The conventional spectral methods are further modified by using piecewise high order rational interpolants on the Chebyshev mesh within each sub-domain with the boundary domain placed at the strike price where the discontinuity is located. The resulting system is then solved by applying Laplace transform method through deformation of a contour integral. Firstly, we use this approach to price plain vanilla options and then extend it to price options described by a jump-diffusion model, barrier options and the Heston’s volatility model. To approximate the integral part in the jump-diffusion model, we use the Gauss-Legendre quadrature method. Finally, we carry out extensive numerical simulations to value these options and associated Greeks (the measures of sensitivity). The results presented in this thesis demonstrate the spectral accuracy and efficiency of our approach, which can therefore be considered as an alternative approach to price these class of options.