The Heston-Hull-White (HHW) model is a generalization of the classical Heston approach that incorporates stochastic interest rates, making it a more accurate representation of financial markets. In this work, we investigate a computational procedure via a three-dimensional partial differential equation (PDE) to solve option pricing problems under the HHW framework. We propose a local radial basis function–finite difference (RBF–FD) framework under the integration of a new variant of the multiquadric function for efficiently resolving the model. Our study highlights the error analysis of the proposed weights for the first and second derivatives of a suitable function and demonstrates the effectiveness of the RBF–FD approach for this high-dimensional financial model.

The matrix sign function has a key role in several applications in numerical linear algebra. This paper presents a novel iterative approach with a sixth order of convergence to efficiently compute this function. The scheme is constructed via the employment of a nonlinear equations solver for simple roots. Then, the convergence of the extended matrix procedure is investigated to demonstrate the sixth rate of convergence. Basins of attractions for the proposed solver are given to show its global convergence behavior as well. Finally, the numerical experiments demonstrate the effectiveness of our approach compared to classical methods.

This article introduces a multi-step solver for sets of nonlinear equations. To achieve this, we consider and develop a multi-step Steffensen-type method without memory, which does not require evaluations of the Fréchet derivatives, and subsequently extend it to a method with memory. The resulting order is $\sqrt{5}+2$, utilizing the identical number of functional evaluations as the solver without memory, thereby demonstrating a higher computational index of efficiency. Finally, we illustrate the advantages of the proposed scheme with memory through various test problems.

We develop a high-order compact numerical scheme for solving a nonlinear Black–Scholes equation arising in option pricing under transaction costs. By leveraging a Hermite-enhanced Radial Basis Function-Finite Difference (RBF-HFD) method with three-point stencils, we achieve fourth-order spatial accuracy. The fully nonlinear PDE, driven by Gamma-dependent volatility models, is discretized via RBF-HFD in space and integrated using an explicit sixth-order Runge–Kutta scheme. Numerical results confirm the proposed method’s accuracy, stability, and its capability to capture sharp gradient behavior near strike prices.

In the article, we propose a combination method based on the multigrid method and constraint data to solve the inverse problem in the context of the nonlinear convection–diffusion equation in the multiphase porous media flow. The inverse problem consists of a data-fitting term involving the discretization of a direct problem, a constraint term concerning the incorporation of constraint data, and a regularization term dealing with the improvement of stability. A multigrid method, which is specialized for large-scale problems and works by keeping the consistence of objective functionals between different grids, is applied in the process of inversion. Based on the numerical results, the proposed combination method has the advantages of fast calculation, high precision, good stability, and strong anti-noise ability in computation. It obtains good performance under various noise levels, as well as outperforming any one method used alone.

We develop a novel Steffensen-type iterative solver to solve nonlinear scalar equations without requiring derivatives. A two-parameter one-step scheme without memory is first introduced and analyzed. Its optimal quadratic convergence is then established. To enhance the convergence rate without additional functional evaluations, we extend the scheme by incorporating memory through adaptively updated accelerator parameters. These parameters are approximated by Newton interpolation polynomials constructed from previously computed values, yielding a derivative-free method with R-rate of convergence of approximately 3.56155. A dynamical system analysis based on attraction basins demonstrates enlarged convergence regions compared to Steffensen-type methods without memory. Numerical experiments further confirm the accuracy of the proposed scheme for solving nonlinear equations.