Abstract
The Heston-Hull-White three-dimensional time-dependent partial differential equation (PDE) is one of the important models in mathematical finance, at which not only the volatility is modeled based on a stochastic process but also the rate of interest is assumed to follow a stochastic dynamic. Hence, an efficient method is derived in this paper based on the methodology of the localized radial basis function generated finite difference (RBF-FD) scheme. The proposed solver uses the RBF-FD approximations on graded meshes along all three spatial variables and a high order time-stepping scheme. Stability is also studied in detail to show under what conditions the proposed method is stable. Computational simulations are given to support the theoretical discussions.