Soleymani, FazlollahAkgul, AliAkgul, Esra Karatas2024-12-242024-12-2420191224-17841844-0835https://doi.org/10.2478/auom-2019-0042https://hdl.handle.net/20.500.12604/7686The aim of this work is to tackle the three-dimensional (3D) Heston-Cox-Ingersoll-Ross (HCIR) time-dependent partial differential equation (PDE) computationally by employing a non-uniform discretization and gathering the finite difference (FD) weighting coefficients into differentiation matrices. In fact, a non-uniform discretization of the 3D computational domain is employed to achieve the second-order of accuracy for all the spatial variables. It is contributed that under what conditions the proposed procedure is stable. This stability bound is novel in literature for solving this model. Several financial experiments are worked out along with computation of the hedging quantities Delta and Gamma.eninfo:eu-repo/semantics/openAccessFinancial option pricingstochastic interest rateHeston modelnon-uniform finite difference methodquadratically convergentArticle273207230Q2WOS:000498863500013Q32-s2.0-8507815714710.2478/auom-2019-0042