Qualitatively Stable Nonstandard Finite Difference Scheme for Numerical Solution of the Nonlinear Black–Scholes Equation

In this paper, we use a numerical method for solving the nonlinear Black–Scholes partial differential equation of the European option under transaction costs, which is based on the nonstandard discretization of the spatial derivatives. The proposed scheme, in addition to the unconditional positivity, is stable, consistent, and monotone. In order to illustrate the efficiency of the new method, numerical results have been performed by four models.