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. )e 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.