An extremely efficient numerical method for pricing options in the Black–Scholes model with jumps

We propose a new numerical method for pricing options in the Black–Scholesmodel with jumps. Specifically, we consider the partial integro-differential problem that yields the option price, and we solve it by means of a finitedifference scheme that combines a fixed-point iteration technique and arepeated space-time Richardson extrapolation procedure. Such an approachturns out to be not only extremely accurate and fast but also very simple toimplement, since the use of fast convolution techniques for handling the jumpintegral is not required. Numerical experiments are presented in which vanilla,barrier, and American options are considered.