AN EFFICIENT AND ROBUST NUMERICAL METHOD FOR OPTION PRICES IN A TWO-ASSET JUMP-DIFFUSION MODEL

We present an efficient and robust finite difference method for a two-asset jump diffusion model, which is a partial integro-differential equation (PIDE). To speed up a computational time, we compute a matrix so that we can calculate the non-local integral term fast by a simple matrix-vector operation. In addition, we use bilinear interpolation to solve integral term of PIDE. We can obtain more stable value by using the payoff-consistent extrapolation. We provide numerical experiments to demonstrate a performance of the proposed numerical method. The numerical results show the robustness and accuracy of the proposed method.