Numerical investigation of high-dimensional option pricing PDEs by utilizing a hybrid radial basis function- finite difference procedure

The target of this research is to resolve high -dimensional partial differential equations (PDEs) for multi -asset options, modeled as parabolic time -dependent PDEs. We present a hybrid radial basis function - finite difference (RBF-FD) solver, which combines the advantages of Gaussian and multiquadric functions. Additionally, we employ the Krylov subspace method on the resulting system of ordinary differential equations, reducing the computation load for finding the numerical solution. Computational tests, up to seven dimensions, and comparisons support the superiority of our hybrid solver.