Solving partial integro-differential option pricing problems for a wide class of infinite activity Lévy processes

In this paper, numerical analysis of finite difference schemes for partial integro-differential models related to European and American option pricing problems under a wide class of Lévy models is studied. Apart from computational and accuracy issues, qualitative properties such as positivity are treated. Consistency of the proposed numerical scheme and stability in the von Neumann sense are included. Gauss–Laguerre quadrature formula is used for the discretization of the integral part. Numerical examples illustrating the potential advantages of the presented results are included.