Journal of Mathematics
Volume 2021 (2021), Article ID 6679484, 12 pages
https://doi.org/10.1155/2021/6679484

Research Article

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

Mohammad Mehdizadeh Khalsaraei,¹ Ali Shokri,¹ Zahra Mohammadnia,¹ and Hamid Mohammad Sedighi²

¹Department of Mathematics, Faculty of Science, University of Maragheh, Maragheh, Iran
²Mechanical Engineering Department, Faculty of Engineering, Shahid Chamran University of Ahvaz, Ahvaz, P.O. Box 61357-43337, Iran

Correspondence should be addressed to Mohammad Mehdizadeh Khalsaraei; mehdizadeh@maragheh.ac.ir

Received 19 December 2020; Accepted 26 April 2021; Published 11 May 2021

Academic Editor: Efthymios G. Tsionas

Copyright © 2021 Mohammad Mehdizadeh Khalsaraei et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

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.

1. Introduction

Financial mathematics is a branch of applied mathematics that deals with financial markets. Option contracts are the most traditionally financial instruments. Hence, with the increasing demand for this financial instrument, the pricing of contracts has very important role to play in the economy. An option on an underlying asset traded on financial markets is an agreement that gives its holder the option to buy or sell the asset mentioned in the contract on a specified date at a certain and predetermined price. The date specified in the contract also calls to the maturity date or expiration date, and the price specified in the contract is called the strike price or exercise price. The biggest advantage of an option contract is limited loss and unlimited profits. Option contracts are divided into two types such as European and American types according to the time of exercise. A European option is only applicable at the time specified in the contract, but the American option is applicable at any time until maturity date . One of the models in the discussion of financial derivatives pricing is the Black–Scholes model. This model was introduced by Fisher Black and Miran Scholes and developed by Robert Merton [1]. It is represented by a partial differential equation in the following form:

where is the value of the option, is the current underlying asset price, t = t is the time, is the expiration date of the option, is the risk-free interest rate, and is the volatility (also called implicit volatility). and are fixed. Partial differential equation paid attention in the field of economics with the Black–Scholes model. Therefore, in order to obtain the price of these contracts, we need a simple and accurate method to solve this differential equation. The standard Black–Scholes equation is considered in a complete market under certain simple assumptions such as the cash market, without transaction cost. Also, the most important feature in the classical Black–Scholes model, the volatility of an underlying asset is assumed to be constant. In 1987, all over the world, stock the value decreased observably in the short course of time. After this decrease, analysts have concluded that the options with profits and losses are traded with higher implicit volatility versus the option that do not have profit and loss. So, this fact is not compatible with the assumption of constant volatility in the Black–Scholes model. In addition to despite the transaction cost, the Black–Scholes equation changes to an equation with a nonlinear volatility function, which depends on the time, current underlying asset price, and the price derivatives of the option relative to the underlying asset price. Then, the classic Black–Scholes equation changes to the following nonlinear Black–Scholes equation:

in which is the modified volatility. Various nonlinear volatility models have been presented in [2–4]. We consider the nonlinear Barles and Soner model which is obtained as follows:

where , in which is the risk-adjusted factor and is the number of options that are sold. The function is the solution of the following nonlinear initial value problem [5]:

For function , an exact implicit solution is presented as follows [6]:

Let us consider equation (2) with the nonlinear modified volatility of Barles and Soner (3) and (4). By applying the variable , equation (1) is converted to the following equation:

for and ,

as in [5] we consider the initial condition, the first boundary condition, and the second boundary condition for and as follows:

where , and denote the strike prices of the options, is a positive constant, and is the Heaviside function. Conditions and are established.

So far, no classical solution has been given for equation (2), but the viscosity solution is used as the solution of equation (2).

The numerical solution of the nonlinear Black–Scholes equation has been considered in some works. For example, in [7], an explicit method has been presented that preserves positivity property and is unconditional stable but not consistent under any conditions. A variational principle is established for the generalized KdV-Burgers equation by the semiinverse method in [8–11]. In [12–14], the homotopy perturbation method is designed to obtain a solution to the Black–Scholes equation with boundary conditions for a European option pricing problem and nonlinear problems. Methods based on the variational iteration method (VIM) in [15–17] are proposed for the Black–Scholes model. In [18], an explicit nonstandard finite difference method has been proposed that it is stable and conditionally acceptable. Also, a finite difference method has been proposed which is the stability condition for it and is established by limiting the length of the step time and the spatial step size. In [5], a completely implicit time step method has been presented for (6) and represented convergence approximate solution to the viscosity solution and used Newton’s algorithm to solve a nonlinear system. They caused numerical instability in the solution that has been resolved by an iterative method. This process leads to very expensive computing costs. Also in [6], the MOL method has been used to solve (6). In [19, 20], equation (6) is solved by the variable which in addition to having very expensive computational cost, convergence, and other results of this work is under relatively unrealistic conditions or limits the step-size time and space.

In [21], considering a class of nonlinear option pricing models, Newton and Picard iterative procedures are proposed for solving the nonlinear systems of algebraic equations. Also, to improve the computational fast, two-grid algorithms are developed. In [22], combination of a sixth-order finite difference scheme in space and a third-order strong stability preserving Runge–Kutta in time has been used to obtain numerical solutions of the linear and nonlinear European put option models. In [23], a new fourth-order finite difference approximation is contributed and equipped with the fourth-order Runge–Kutta scheme for the nonlinear Black–Scholes. They proved that under several criteria, the procedure is time stable. Also, the proposed technique reduces the computational cost when more accurate results are requested. In [24], nonlinear generalization of the Black– Scholes equation for pricing and American-style call options is investigated. They proposed a numerical method that involves transforming the free boundary problem for a nonlinear Black–Scholes equation into the so-called Gamma variational inequality with a new variable depending on the Gamma of the option. Also, a modified projected successive overrelaxation method in order to construct a numerical scheme for discretization of the Gamma variational inequality is presented.

In this paper, we solve equation (6) by using an explicit method that has very low computational cost, and the method is free of any spurious oscillations. It is unconditionally positive, stable, consistent, and monotone. We show that the solution of the proposed method converges to the solution of the viscosity of the equation.

The rest of the paper is organized as follows. In section 2, we give a brief summary of the NSFD methods. In section 3, we examine the application of nonstandard methods through spatial nonlocalized discretization for the nonlinear Black–Scholes under transaction costs. In section 4, we analyze the new method. The numerical results are given in section 5.

2. Nonstandard Finite Difference Scheme

In this section, we give some preliminaries including nonstandard finite difference methods. The initial foundation of NSFD schemes came from the exact finite difference schemes. Numerical methods based on the standard finite difference approach are consistent with the original differential equation and guarantee convergence of the discrete solution to the exact one, but they impose severe restrictions on the time step, and essential qualitative properties of the solution are not transferred to the numerical solution. One response to this situation was the initiation by Mickens [25] of a research program for the investigation of new methods for constructing finite difference schemes that are convergent for any step size. Nonstandard finite difference methods (NSFDs) in addition to the usual properties of consistency, stability, and hence convergence produce numerical solutions that also exhibit essential properties of solutions [3–5, 13, 15, 16, 18, 21, 26–38].

This class of schemes and their formulations centers on two issues: first, how should discrete representations for derivatives be determined and second what are the proper forms to be used for nonlinear and reaction terms. The forward finite difference method is one of the simplest discretization schemes. In this method, the derivative is replaced by , where is step size. However, in the Mickens schemes, this term is replaced by , where is an increasing continuous function of , and the function satisfies the following conditions:

Note that in taking to obtain the derivative, the use of any of this will lead to the usual result for the first derivative:

where and are continuous functions of the step size . A scheme is called a nonstandard finite difference method if at least one of the following conditions is met:

(i)In the discrete derivatives, the traditional denominator is replaced by a nonnegative function such that

 or

 Examples of functions that satisfy these conditions are as follows [25]:

(ii)Nonlinear terms are approximated in a nonlocal way, i.e., by a suitable function of several points of the mesh. If there are reaction terms of the differential equation [3–5, 13, 15, 16, 18, 21, 26–28], these are replaced by

One can say that there is no appropriate general method to choose the function or to choose which nonlinear terms are to be replaced, and some special techniques may be found in [3–5, 13, 15, 16, 18, 21, 26–28].

3. Construction of the New NSFD

In this section, by using the mentioned rules for NSFD schemes, we propose our new NSFD scheme. Let us consider the computational domain and discretize it in the following form. We introduce a grid of mesh points where , , and the spatial step-size is given by , and the time step size is . We denote the approximation of by . The new NSFD is of the following form:

The explicit form of (17) is as follows:

The finite difference approximation provides a difference equation as follows:

with (8)–(10), and we define the initial and boundary conditions as follows:

The matrix form (18) for all and is given in the following form:

in which is a tridiagonal matrix with positive components and is known as boundary values. The new scheme can be written in the following form:

where

4. Analysis of the New Method

In this section, we investigate the monotonicity, positivity preserving, stability, consistency, and convergency properties of the new method.

For and , we define the function as follows:

where

Lemma 1 (monotonicity). The new NSFD method is monotone.

Proof. As in [5] which is enough to show that for each and ,

Since , , and , the first three sentences of the right of the relationship (26), respectively, are nonincreasing in , increasing in , and decreasing in .
Suppose is a column vector . We have

Assume that we take into account the nonlinear sentence is to the right of the relation (26). From (3), we have

where

for each and ; independent of , by derivative phrase relative to the variable and using (5) and (24) and , we have

So, is an increasing function of . Using the above symbols and by monotone combination of and properties of the linear sentences of the relationship (26) and (30), we have

This is the same as (28). Similarly from monotonicity and (31), it can be shown that (29) is confirmed, and this completes the proof.
Now, we investigate the positivity property of the constructed methods in the previous section. With positivity, we mean that the component-wise non-negativity of the initial vector is preserved in time for the approximated solution. Several NSFD schemes have been constructed in the literature (see, for example, [1, 25, 29, 39–42]).

Proposition 1 (positivity). Suppose , , and are the real numbers that are nonnegative, in this case, method (18) offers a nonnegative approximation of for the solution of equation (2).

Proof. Since the parameter

also

then

Therefore, this completes the proof.

Lemma 2 (stability). The new NSFD method is unconditionally stable.

Proof. As in [5], for each , we know , that is the solution of the NSFD method, and then for the initial and boundary conditions, , , and are defined, and just shows

represents infinite norm. For each and , we have

From (23)–(25), it is clear that , , , so we will have

Now, for , let us consider the following two modes:

(i)Assume , then since ,

 Therefore,

(ii)Suppose that , then if or , it is seen from the boundary conditions that

By combining two states (43) and (44), stability is achieved and this completes the proof.

Lemma 3 (consistency). The proposed method (18) is consistent, and its local truncation error is .

Proof. The local truncation error for (18) is as follows:

We have Taylor expansions

replacing these expressions in the term as follows:

However, is the solution of the Black–Scholes equation, so the proposed method is consistent:

Therefore, the main part of the local truncation error is as follows:

so . Therefore, this completes the proof.

Theorem 1. The solution of method (18) convergences to the viscosity solution of equation (2), when

Proof. If the discretized method for nonlinear partial differential equations be stable, consistent, and monotone, then the approximate solution converges to the viscosity solution of the equation [43].

5. Numerical Results

To illustrate the effectiveness of the new NSFD, the problem in the independent variables is solved for various values of parameter transaction cost on a number of uniform meshes by the numerical method presented in the previous sections. The numerical solutions are then transformed back in . We have programmed these methods in MATLAB.

Example 1. The first test problem is which we consider European call option contracts for solving equation (6) with initial (8) and boundary conditions (9) and (10) by the following parameters:

Option value with uniform meshes and is drawn up in Figure 1(a). It can be seen that the solution of the proposed method is positive and numerically stable. Also, the comparison of price options with different transaction costs at is shown in Figure 1(b). It can be seen that with the increase in the transaction cost, the price of the option will also increase. This is also true in practice. The value of the parameters used in this section has been taken from [5].

(a)

(b)

Figure 1: (a) European call option price charts of the explicit nonstandard method for and (b) European option price charts of the explicit nonstandard method for different transaction costs in .

Example 2. As our the second numerical experiment, consider put option contracts for solving equation (6) with initial (8) and boundary conditions (9) and (10) by the following parameters:

Option value with uniform meshes and is drawn up in Figure 2(a). As can be seen from Figure 2, the function of the value of the option is an increasing function of the transaction cost parameter. Also, the price of the put option with different transaction costs at is shown in Figure 2(b).

(a)

(b)

Figure 2: (a) European put option price charts of the explicit nonstandard method for and (b) European option price charts of the explicit nonstandard method for different transaction costs in .

Example 3. Let us consider European butterfly spread option contracts for solving equation (6) with initial (8) and boundary conditions (9) and (10) by the following parameters:

Option value with uniform meshes and is drawn up in Figure 3(a). As can be seen from Figure 3, the function of the value of the option is an increasing function and is the transaction cost parameter. Also, the price of the European put option with different transaction costs at is shown in Figure 3(b).

(a)

(b)

Figure 3: (a) Butterfly spread option price charts of the explicit nonstandard method for and (b) butterfly spread option price charts of the explicit nonstandard method for different transaction costs in .

Example 4. In our last example, we consider cash-or-nothing option contracts for solving equation (6) with initial (8) and boundary conditions (9) and (10) by the following parameters:

Option value with uniform meshes and is drawn up in Figure 4(a). As can be seen from Figure 4, the function of the value of the option is an increasing function of the transaction cost parameter. Also, the price of a cash-or-nothing option with different transaction costs at is shown in Figure 4(b).

(a)

(b)

Figure 4: (a) European cash-or-nothing option price charts of the explicit nonstandard method for and (b) cash-or-nothing charts of the explicit nonstandard method for different transaction costs in .

6. Discussion and Conclusion

In this paper, we proposed a new nonstandard finite difference method for the Black–Scholes nonlinear equation under transaction costs, which, in addition to being free from spurious oscillation and preserving positivity property for the desired step size, is unconditionally stable. Also, we showed the convergence of the solution of the new method to the viscosity solution of the equation by proving stability, consistency, and monotonicity and observed that the price of a European option contract is an increasing function of the parameter transaction cost. Our results indicate that a properly implemented version of our scheme is useful for the numerical integration of the considered Black–Scholes equation. In our future research, we intend to solve time fractional Black–Scholes European option pricing equation using the two-scale fractal derivatives.

Data Availability

The data used to support this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

References

R. C. Merton, “Theory of rational option pricing,” The Bell Journal of Economics and Management Science, vol. 4, no. 1, pp. 141–183, 1973.

T. Hoggard, A. W. Whalley, and P. Wilmott, “Hedging option portfolios in the presence of transaction costs,” Adv. Futures Options Res, vol. 7, pp. 21–35, 1994.

M. Jandačka, “And Ševčovič, on the risk-adjusted pricing-methodology-based valuation of vanilla options and explanation of the volatility smile,” Journal of Applied Mathematics, vol. 3, pp. 235–258, 2005.

H. E. Leland, “Option pricing and replication with transactions costs,” The Journal of Finance, vol. 40, no. 5, pp. 1283–1301, 1985.

D. C. Lesmana and S. Wang, “An upwind finite difference method for a nonlinear Black-Scholes equation governing European option valuation under transaction costs,” Applied Mathematics and Computation, vol. 219, no. 16, pp. 8811–8828, 2013.

R. Company, L. Jódar, and J.-R. Pintos, “A numerical method for European Option Pricing with transaction costs nonlinear equation,” Mathematical and Computer Modelling, vol. 50, no. 5-6, pp. 910–920, 2009.

S. Zhou, W. Li, Y. Wei, and C. Wen, “A positivity-preserving numerical scheme for nonlinear option pricing models,” Journal of Applied Mathematics, vol. 2012, Article ID 205686, 20 pages, 2012.

S. Abbas, M. Benchohra, N. Hamidi, and J. J. Nieto, “Hilfer and Hadamard fractional differential equations in Fréchet spaces,” TWMS Journal of Pure and Applied Mathematics, vol. 10, no. 1, pp. 102–116, 2019.

S. Harikrishnan, K. Kanagarajan, and E. M. Elsayed, “Existence and stability results for differential equations with complex order involving Hilfer fractional derivative,” TWMS Journal of Pure and Applied Mathematics, vol. 10, no. 1, pp. 94–101, 2019.

J.-H. He, F.-Y. Ji, and H. Mohammad-Sedighi, “Difference equation vs differential equation on different scales,” International Journal of Numerical Methods for Heat & Fluid Flow, vol. 31, no. 1, pp. 391–401, 2020.

N. Moshtaghi and A. Saadatmandi, “Numerical solution for diffusion equations with distributed-order in time based on sinc-legendre collocation method,” Applied and Computational Mathematics, vol. 19, no. 3, pp. 317–355, 2020.

V. Gülkaç, “The homotopy perturbation method for the Black Scholes equation,” Journal of Statistical Computation and Simulation, vol. 80, no. 12, pp. 1349–1354, 2010.

J.-H. He, “Homotopy perturbation method: a new nonlinear analytical technique,” Applied Mathematics and Computation, vol. 135, no. 1, pp. 73–79, 2003.

K. Trachoo, W. Sawangtong, and P. Sawangtong, “Laplace transform homotopy perturbation method for the two dimensional Black Scholes model with European call Option,” Mathematical and Computational Applications, vol. 22, pp. 1–11, 2017.

J.-H. He and X.-H. Wu, “Variational iteration method: new development and applications,” Computers & Mathematics with Applications, vol. 54, no. 7-8, pp. 881–894, 2007.

P. Kashyap and A. Kumar Sharma, “Numerical solution of Black-Scholes equation using variational iteration method,” International Journal of Innovative Research in Science, Engineering and Technology, vol. 9, no. 1, pp. 12–19, 2020.

H. Rouhparvary and M. Salamatbakhshz, “Analytical solution of the Black-Scholes equation by using variational iteration method,” Applied Mathematics E-Notes, vol. 13, pp. 243–248, 2013.

S. Mashayekhi and J. hugger, “Finite difference schemes for a nonlinear black-scholes model with transaction cost and volatility risk,” Acta Mathematica Universitatis Comenianae, vol. 2, pp. 255–266, 2015.

B. Düring, M. Fournié, and A. Jüngel, “High order compact finite difference schemes for a nonlinear Black-Scholes equation,” International Journal of Theoretical and Applied Finance, vol. 6, no. 7, pp. 767–789, 2003.

P. Heider, “Numerical methods for non-linear black-scholes equations,” Applied Mathematical Finance, vol. 17, no. 1, pp. 59–81, 2010.

M. N. Koleva and L. G. Vulkov, “Fast computational approach to the Delta Greek of non-linear Black-Scholes equations,” Journal of Computational and Applied Mathematics, vol. 340, pp. 508–522, 2018.

S. Gulen, C. Popescu, and M. Sari, “A new approach for the black-scholes model with linear and nonlinear volatilities,” Mathematics, vol. 7, no. 8, p. 760, 2019.

A. Ghanadian and T. Lotfi, “Approximate solution of nonlinear Black-Scholes equation via a fully discretized fourth-order method,” AIMS Mathematics, vol. 5, no. 2, pp. 879–893, 2020.

M. R. Grossinho, Y. F. Kord, and D. Sevcovic, “Pricing American call options using the Black-Scholes equation with a nonlinear volatility function,” Journal of Computational Finance, vol. 23, no. 4, pp. 93–113, 2020.

R. E. Mickens, Nonstandard Finite Difference Models of Differential Equations, World Scientific, Singapore, 1994.

M. Mehdizadeh Khalsaraei and R. Shokri Jahandizi, “A modified explicit method for the Black-Scholes equation with positivity preserving property,” Journal of Mathematical and Computational Science, vol. 15, pp. 299–305, 2015.

M. Mehdizadeh Khalsaraei and R. Shokri Jahandizi, “A positivity preserving schemes for black-scholes equation,” Journal of Finance and Accounting, vol. 6, pp. 101–105, 2015.

M. Mehdizadeh Khalsaraei and R. Shokri Jahandizi, “An efficient nonstandard numerical method with positivity preserving property,” Journal of Mathematical Modeling, vol. 4, no. 2, pp. 161–169, 2016.

M. M. Khalsaraei, A. Shokri, H. Ramos, and S. Heydari, “A positive and elementary stable nonstandard explicit scheme for a mathematical model of the influenza disease,” Mathematics and Computers in Simulation, vol. 182, pp. 397–410, 2021.

M. M. Khalsaraei, “Positivity of an explicit Runge-Kutta method,” Ain Shams Engineering Journal, vol. 6, no. 4, pp. 1217–1223, 2015.

M. M. Khalsaraei and R. S. Jahandizi, “A family of positivity preserving schemes for numerical solution of Black-Scholes equation,” International Journal of Financial Engineering, vol. 3, no. 4, p. 1650025, 2016.

M. Mehdizadeh Khalsaraei and R. Shokri Jahandizi, “Efficient explicit nonstandard finite difference scheme with positivity-preserving property,” Journal of Science, vol. 30, no. 1, pp. 259–268, 2017.

R. Company, E. Navarro, J. Ramón Pintos, and E. Ponsoda, “Numerical solution of linear and nonlinear Black-Scholes option pricing equations,” Computers & Mathematics with Applications, vol. 56, no. 3, pp. 813–821, 2008.

M. Mehdizadeh Khalsaraei, “An improvement on the positivity results for 2-stage explicit Runge-Kutta methods,” Journal of Computational and Applied Mathematics, vol. 235, no. 1, pp. 137–143, 2010.

M. Mehdizadeh Khalsaraei and N. Usmani, “A family of positive nonstandard numerical methods with application to Black-Scholes equation,” Sahand Communications in Mathematical Analysis, vol. 5, no. 1, pp. 31–40, 2017.

M. Milev and A. Tagliani, “Low volatility options and numerical diffusion of finite difference schemes,” Serdica Mathematical Journal, vol. 36, no. n.3, pp. 223–236, 2010.

M. Milev and A. Tagliani, “Nonstandard finite difference schemes with application to finance: option pricing,” Serdica Mathematical Journal, vol. 36, no. 1, pp. 75–88, 2010.

M. Milev and A. Tagliani, “Numerical valuation of discrete double barrier options,” Journal of Computational and Applied Mathematics, vol. 233, no. 10, pp. 2468–2480, 2010.

M. M. Khalsaraei and F. Khodadosti, “Qualitatively stability of nonstandard 2-stage explicit Runge-Kutta methods of order two,” Computational Mathematics and Mathematical Physics, vol. 56, no. 2, pp. 235–242, 2016.

R. E. Mickens and P. M. Jordan, “A new positivity-preserving nonstandard finite difference scheme for the DWE,” Numerical Methods for Partial Differential Equations, vol. 21, no. 5, pp. 976–985, 2005.

R. E. Mickens and P. M. Jordan, “A positivity-preserving nonstandard finite difference scheme for the Damped Wave Equation,” Numerical Methods for Partial Differential Equations, vol. 20, no. 5, pp. 639–649, 2004.

M. Milev and A. Tagliani, “Efficient implicit scheme with positivity preserving and smoothing properties,” Journal of Computational and Applied Mathematics, vol. 243, pp. 1–9, 2013.

G. Barles, “Convergence of numerical schemes for degenerate parabolic equations arising in finance,” in Numerical Methods in Finance, L. C. G. Rogers and D. Talay, Eds., Cambridge University Press, Cambridge, UK, 1997.