“…The celebrated Black-Scholes model 1,2 is based on assumption that the price of the underlying asset behaves like a geometric Brownian motion with a drift and a constant volatility, which cannot explain the market prices of options with various strike prices and maturities. To explain this behavior, a number of alternative models has appeared in the financial literatures, for example, nonlinear models, [3][4][5][6][7][8] stochastic volatility models, [9][10][11][12] jump-diffusion models, [13][14][15][16] regime-switching models, 17,18 and regime-switching jump-diffusion models, 19,20 which are given by coupled partial integro-differential equations (PIDEs). However, these models are more difficult to handle numerically in contrast to the celebrated Black-Scholes model.…”