A Radial Basis Function Scheme for Option Pricing in Exponential Lévy Models

We use Radial Basis Function (RBF) interpolation to price options in exponential Lévy models by numerically solving the fundamental pricing PIDE (Partial integro-differential equations). Our RBF scheme can handle arbitrary singularities of the Lévy measure in 0 without introducing further approximations, making it simpler to implement than competing methods. In numerical experiments using processes from the CGMY-KoBoL class, the scheme is found to be second order convergent in the number of interpolation points, including for processes of unbounded variation.