A simple recursive numerical method for Bermudan option pricing under Lévy processes

A numerical method is developed that can price options, including exotic options that can be priced recursively such as Bermudan options, when the underlying process is an exponential Lévy process with closed form conditional characteristic function. The numerical method is an
extension of a recent quadrature option pricing method so that it can be applied with the use of fast Fourier transforms. Thus the method possesses desirable features of both transform and quadrature option pricing techniques since it can be applied for a very general set of underlying
Lévy processes and can handle certain exotic features. To illustrate the method it is applied to European and Bermudan options for a log normal process, a jump diffusion process, a variance gamma process and a normal inverse Gaussian process.