This paper considers the problem of numerically evaluating American option prices when the dynamics of the underlying are driven by both stochastic volatility following the square root process of Heston , and by a Poisson jump process of the type originally introduced by Merton . We develop a method of lines algorithm to evaluate the price as well as the delta and gamma of the option, thereby extending the method developed by Meyer  for the case of jump-diffusion dynamics. The accuracy of the method is tested against two numerical methods that directly solve the integro-partial differential pricing equation. The first is an extension to the jump-diffusion situation of the componentwise splitting method of Ikonen and Toivanen . The second method is a Crank-Nicolson scheme that is solved using projected successive over relaxation and which is taken as the benchmark for the price. The relative efficiency of these methods for computing the American call option price, delta, gamma and free boundary is analysed. If one seeks an algorithm that gives not only the price but also the delta and gamma to the same level of accuracy for a given computational effort then the method of lines seems to perform best amongst the methods considered.
|Number of pages
|International Journal of Theoretical and Applied Finance
|Published - 1 May 2009