Next Article in Journal
Modeling of Cross Work Hardening and Apparent Normality Loss after Biaxial–Shear Loading Path Change
Next Article in Special Issue
Subordination Involving Regular Coulomb Wave Functions
Previous Article in Journal
SCRnet: A Spatial Consistency Guided Network Using Contrastive Learning for Point Cloud Registration
Previous Article in Special Issue
Multivalent Prestarlike Functions with Respect to Symmetric Points
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Nonstandard Finite Difference Method for a Generalized Black–Scholes Equation

by
Mohammad Mehdizadeh Khalsaraei
1,*,
Mohammad Mehdi Rashidi
2,3,
Ali Shokri
1,
Higinio Ramos
4 and
Pari Khakzad
1
1
Department of Mathematics, Faculty of Science, University of Maragheh, Maragheh 83111-55181, Iran
2
Institute of Fundamental and Frontier Sciences, University of Electronic Science and Technology of China, Chengdu 610054, China
3
Faculty of Mechanical and Industrial Engineering, Quchan University of Technology, Quchan 94771-67335, Iran
4
Scientific Computing Group, Universidad de Salamanca, Plaza de la Merced, 37008 Salamanca, Spain
*
Author to whom correspondence should be addressed.
Symmetry 2022, 14(1), 141; https://doi.org/10.3390/sym14010141
Submission received: 17 October 2021 / Revised: 3 November 2021 / Accepted: 16 December 2021 / Published: 12 January 2022
(This article belongs to the Special Issue Complex Analysis, in Particular Analytic and Univalent Functions)

Abstract

:
An implicit finite difference scheme for the numerical solution of a generalized Black–Scholes equation is presented. The method is based on the nonstandard finite difference technique. The positivity property is discussed and it is shown that the proposed method is consistent, stable and also the order of the scheme respect to the space variable is two. As the Black–Scholes model relies on symmetry of distribution and ignores the skewness of the distribution of the asset, the proposed method will be more appropriate for solving such symmetric models. In order to illustrate the efficiency of the new method, we applied it on some test examples. The obtained results confirm the theoretical behavior regarding the order of convergence. Furthermore, the numerical results are in good agreement with the exact solution and are more accurate than other existing results in the literature.

1. Introduction

An option is the right, but not the obligation, to buy or sell an underlying asset at a specific price on or before a certain date at a fixed strike price. European call (put) option and American call (put) option are two types of known options. European option can be exercised only on expiration date T. For American option, exercise is permitted at any time t T . It was shown in [1] that option pricing can be modeled using a partial differential equation of second order. The resulting equation is called the Black–Scholes equation. If the coefficients are constant or independent of the spatial variable, it can be transformed into the standard heat equation and can be solved exactly. However, in the financial market, the coefficients can depend on the time and the asset price. In such situation the analytical solution of the generalized Black–Scholes model is not available. Therefore, numerical simulations are essential to obtain information about the behavior of the solutions. It is important to develop numerical methods that reproduce the qualitative properties of the solution.
Here, we consider the generalized Black–Scholes model as given in [2]:
C t 1 2 x 2 σ 2 ( x , t ) 2 C x 2 x r ( t ) C x + r ( t ) C = 0 , ( x , t ) R + × ( 0 , T ) ,
with the following conditions:
C ( x , T ) = max ( x K , 0 ) , x R + ,
C ( 0 , t ) = 0 , t [ 0 , T ] .
Equation (1) is a partial differential equation used for pricing options. The function C ( x , t ) is the European call option price with asset price x and time t. T is the maturity, K is the price of strike, r ( t ) is the risk-free interest rate, and σ ( x , t ) is the volatility function of the underlying asset and time. When the coefficients are constant functions, we get the classical Black–Scholes model. Following [2] we consider the above model on a truncated domain R = ( 0 , S max ) × ( 0 , T ) , where S max > K is a suitably chosen positive number so that
C * t 1 2 x 2 σ 2 ( x , t ) 2 C * x 2 x r ( t ) C * x + r ( t ) C * = 0 , ( x , t ) R
C * ( x , T ) = max ( x K , 0 ) , x [ 0 , S max ] ,
C * ( 0 , t ) = 0 , t [ 0 , T ] ,
C * ( S max , t ) = S max K e t T r ( s ) d s , t [ 0 , T ] .
The result about existence and uniqueness of a solution of (4)–(7) when σ ( x , t ) and r ( t ) are constant functions can be found in [3,4]. It is proved in [5] that if C and C * are solutions of (1)–(3) and (4)–(7), respectively, then
| C ( x , t ) C * ( x , t ) | K exp ln S max x 2 2 min R ¯ σ 2 ( T t ) , ( x , t ) R .
The terminal condition (5) is not a smooth one. Therefore, to guarantee the convergence of the numerical solution, we substitute the terminal condition by a four order smooth function (see in [2]) and obtain the following model:
u t 1 2 x 2 σ 2 ( x , t ) 2 u x 2 x r ( t ) u x + r ( t ) u = 0 , ( x , t ) R ,
u ( x , T ) = ψ ε ( x K ) , x [ 0 , S max ] ,
u ( 0 , t ) = 0 , t [ 0 , T ] ,
u ( S max , t ) = S max K e t T r ( s ) d s , t [ 0 , T ] ,
where the function ψ ε is defined as follows:
ψ ε ( y ) = y , y > ε , c 0 + c 1 y + + c 9 y 9 , ε y ε , 0 , y < ε ,
for a sufficiently small ε > 0 , ψ ε ( y ) satisfying
ψ ε ( ε ) = ψ ε ( ε ) = ψ ε ( ε ) = ψ ε ( ε ) = ψ ε ( 4 ) ( ε ) = 0 , ψ ε ( ε ) = ε , ψ ε ( ε ) = 1 , ψ ε ( ε ) = ψ ε ( ε ) = ψ ε ( 4 ) ( ε ) = 0 .
By using these conditions one can easily find that
c 0 = 35 256 ε , c 1 = 1 2 , c 2 = 35 64 ε , c 4 = 35 128 ε 3 , c 6 = 7 64 ε 5 , c 8 = 5 256 ε 7 , c 3 = c 5 = c 7 = c 9 = 0 .
The proof of the existence and uniqueness of a solution of (9)–(12) can be found in [3]. A proof of the following result can also be found there:
Proposition 1.
There exists a positive constant P independent of ψ ε ( x K ) such that
| C * ( x , t ) u ( x , t ) | P ψ ε ( x K ) max ( x K , 0 ) L
for ( x , t ) R ¯ .
From (8) and (16) we have that we can take the solution of the modified models (9)–(12) closer to that of the original models (1)–(3) by choosing a sufficiently large value for S max and a sufficiently small value for ε .
Traditionally, important requirements for constructing a numerical method are the investigation of the consistency of the discrete scheme with the original differential equation and linear stability analysis with smooth solutions [6,7,8,9,10]. These requirements are important, because they are necessary to ensure the convergence of the discrete solution to the exact one, but for the explicit numerical methods based on standard finite difference approaches. However, most of the time, the qualitative properties of the true solution are not reproduced by the numerical solution. Therefore, this might result in a calamitous erroneous outcome. One way to overcome this disadvantage is to use nonstandard finite difference methods (NSFDs). The origin of the NSFD methods began with Mickens in 1989 [11]. A summary of the known results up to 1994 are given in [12]. The simplicity in construction and the possiblity of their extensions allowed researchers to apply NSFD schemes to solve several complex differential equation models that arise at different disciplines. A scheme is called nonstandard if at least one of the following conditions is satisfied:
1.
The usual denominator h or h 2 of the discrete first or second order derivatives is replaced by a non-negative function Φ such that
Φ ( h ) = h + o ( h 2 ) or Φ ( h ) = h 2 + o ( h 3 ) as 0 < h 0 .
Some examples of Φ ( h ) satisfying these conditions are
h , sin h , 1 e λ h λ or 4 sin 2 ( h 2 ) .
2.
The nonlinear terms that appear in the differential equation are approximated in a non-local way, for instance, the nonlinear terms u 2 and u 3 can be replaced by (see [13]):
u 2 a u i 2 + b u i u i + 1 , a + b = 1 , a , b R , u 3 a u i 3 + ( 1 a ) u i 2 u i + 1 , a R .
In addition to the usual properties of consistency, stability and convergence, the NSDF methods give rise to numerical solutions that also maintain essential qualitative properties of the solutions [14,15,16,17,18,19,20,21].
In this paper, we introduce a new scheme for the numerical solution of the generalized Black–Scholes equation. We propose a non-standard finite difference method, where we use a non-local approximation of the reaction term of the generalized Black–Scholes equation, combined with an implicit time step technique. The method is conditionally stable, positivity preserving and of second-order of convergence with respect to the space variable. Our proposed method has some advantages: A uniform mesh is used; thus, the analysis is simple and it has better accuracy and convergence rate.
The rest of the paper is organized as follows. In Section 2, we review the works done on the generalized Black–Scholes equation and the real options. In Section 3, we propose a NSFD scheme for such a model and we present the analysis of the method concerning the positivity preserving property, Stability, and truncation error. Furthermore, the numerical results obtained from the new method are presented and compared with the results in [2]. The last section describes the features of the new method, the results obtained, and ends with the presentation of ideas for future works.

2. Literature Review

In this section, we review the works done on the generalized Black–Scholes equation. In addition, at the end of the section, we will mention some of the work done related to real options. There are several methods for solving the generalized Black–Scholes equation in the literature. Cen and Le [2] presented an accurate finite difference method, based on a central difference discretization in the space variable using a piecewise uniform mesh and an implicit method in the temporal variable. Cen et al. [22] proposed an exponential scheme in the temporal variable combined with a central difference approach on a spatial piecewise uniform mesh. Kadalbajoo et al. [23] suggested a cubic B-spline collocation method for the generalized model with second order accuracy in both space and time. A fitted finite volume method was used in [24] for the generalized model, showing a second-order convergence. A Cubic spline method was proposed in [25] for a generalized Black–Scholes equation. In that work, the implicit Euler method was used in the temporal discretization, together with a cubic polynomial spline approximation for the spatial discretization. In [26], Mohammadi achieved third order of accuracy in space considering a quintic B-spline collocation method and first- and second-order of accuracy in time using the backward Euler method and the Crank–Nicolson method, respectively. In [27], a high-order numerical method was obtained. In this method a two-step backward differentiation formula is used in the temporal discretization while a high-order difference approximation was used in the spatial discretization. This gives a second order of accuracy in time and third order of accuracy in space.
In [28], Clevenhaus et al. added a penalty term to the American option problem, by doing this the problem was reduced to a PDE on a fixed domain. Their presented penalty term is bounded by an initial guess of the free boundary value and is thus independent of the solution at the last time step. The novelty and advantage of their approach consist in introducing a time-dependent penalty function which led to the construction of an efficient adaptive numerical scheme. In [29] a new method for pricing American options by Boire et al. is presented. The numerical results show that their method successfully reduces the bias plaguing the standard importance sampling method across a wide range of moneyness and maturities, with negligible change to estimator variance.
In [30], considering several uncertainty factors: a renewable energy production amount, an R and D investment amount, a unit price, a risk-free interest rate, and the R and D investment in the six types of renewable energy sources: biomass energy, waste energy, photovoltaic energy, solar thermal energy, marine energy, and wind power energy are evaluated. They utilize the system dynamics approach using the Black–Scholes model and conducts sensitivity analyses of the uncertainty factors. Larissa et al. [31] investigate the contribution of adjusted net savings to sustainable economic growth for 10 Central and Eastern European and Baltic nations using data corresponding to the period 2005–2016. Their results indicated that adjusted net savings impacted on the GDP across the 10 countries analyzed. A bibliometric analysis to examine the features of Carbon capture and storage (CCS) literature including the research focus and trends as well the real options uncertainty and models, types of options, and valuation techniques in [32] are presented. Their results provide energy and environmental policy-makers and CCS project planners with valuable insights into various aspects of CCS policy and project design. For real options in the field of agriculture, we can refer to the article [33] in which the authors investigated a hydroponic farm considering uncertainty, sustainability, and the system’s utility in the dominant desert geography.

3. Method and Main Results

3.1. Construction of the New Scheme

In this subsection, by using the strategy of nonstandard discretization methods for generalized Black–Scholes models (9)–(12), we consider the forward difference for u t , the centered differences for u x , 2 u x 2 , and ξ ( U i + 1 j + U i 1 j ) + ( 1 2 ξ ) ( U i + 1 j + 1 + U i 1 j + 1 ) to approximate the u in the reaction term. Thus, we propose a new discretization scheme for (9) given by
U i j + 1 U i j Δ t 1 2 ( σ i j ) 2 x i 2 U i 1 j 2 U i j + U i + 1 j h 2 r j x i U i + 1 j U i 1 j 2 h + r j ξ ( U i + 1 j + U i 1 j ) + ( 1 2 ξ ) ( U i + 1 j + 1 + U i 1 j + 1 ) = 0 ,
where U i j , σ i j and r j are approximations of u ( x , t ) , σ ( x , t ) and r ( t ) , respectively at the grid point ( x i , t j ) with x i = i h , t j = j Δ t ( i = 1 , , N 1 ; j = 0 , , M 1 ) and h = S max N , Δ t = T M .
The above discretization may be written for each j in a matrix form as
A j U j + d j = B j U j + 1 + d j + 1 ,
where U j = ( U 1 j , , U N 1 j ) T and U j + 1 = ( U 1 j + 1 , , U N 1 j + 1 ) T are vectors containing approximations to u at time levels t j and t j + 1 , respectively. Furthermore, d j and d j + 1 are column vectors that contain the known boundary-values and zeros, and A j and B j are tridiagonal matrices of the form
A j = 1 2 σ i j x i h 2 + r j x i 2 h + ξ r j ; 1 Δ t + σ i j x i h 2 ; 1 2 σ i j x i h 2 r j x i 2 h + ξ r j ,
B j = ( 1 2 ξ ) r j ; 1 Δ t ; ( 1 2 ξ ) r j .

3.2. Characteristics of the New Scheme

This subsection is concerned with the positivity preservation, stability, consistency and convergence properties of the proposed method. Firstly, we present some definitions and results that will be needed in the proof of Theorem 3.
Definition 1.
A matrix A = ( a i j ) is called essentially positive if A is irreducible and a i j 0 , i j .
Theorem 1.
([34]). A matrix A Z n × n is a nonsingular M-matrix if and only if A 1 exists and A 1 0 .
Definition 2.
A matrix A = ( a i j ) is called an L-matrix if a i i > 0 , i N and a i j 0 , i j .
Theorem 2.
([34]). Let A be an L-matrix which is strongly row or column diagonally dominant, i.e., A e > 0 or e T A > 0 T . Then, A is a nonsingular M-matrix.
Theorem 3.
If for i = 1 , , N 1 ; j = 0 , , M 1 it is
1 2 ξ σ i j x i h 2 r j x i h 2 r j ,
then, the method given in (21) is positivity preserving.
Proof. 
From (21), we have
U j = ( A j ) 1 B j U j + 1 + ( A j ) 1 ( d j + 1 d j ) ,
In matrix form, we have
b 1 j c 1 j a 2 j b 2 j c 2 j a N 2 j b N 2 j c N 2 j a N 1 j b N 1 j A j U 1 j U N 1 j + a 1 j U 0 j c N 1 j U N j d j = e 1 j f 1 j d 2 j e 2 j f 2 j d N 2 j e N 2 j f N 2 j d N 1 j e N 1 j B j U 1 j + 1 U N 1 j + 1 + d 1 j U 0 j + 1 f N 1 j U N j + 1 d j + 1
where
a i j = 1 2 σ i j x i h 2 + r j x i 2 h + ξ r j , b i j = 1 t + σ i j x i h 2 c i j = 1 2 σ i j x i h 2 r j x i 2 h + ξ r j d i j = ( 1 2 ξ ) r j , e i j = 1 t , f i j = ( 1 2 ξ ) r j , i = 1 , , N 1 .
As the final condition is positive, if ( A j ) 1 > 0 , B j 0 , and ( d j + 1 d j ) > 0 , then using backward recursion we get that U j is positive for every j. The condition ( A j ) 1 > 0 holds if A j is an M-matrix.
To be A j an M-matrix it should happen that
1 2 σ i j x i h 2 + r j x i 2 h + ξ r j 0 ,
1 2 σ i j x i h 2 r j x i 2 h + ξ r j 0 .
From (28), we should have
ξ σ i j x i h 2 r j x i h 2 r j ,
and from (29), we should have
ξ σ i j x i h 2 + r j x i h 2 r j .
On the other hand, for the non-negativity of B j , it is sufficient to have
( 1 2 ξ ) r j 0 ξ 1 2 .
Then, from (30)–(32), if we choose ξ so that
1 2 ξ σ i j x i h 2 r j x i h 2 r j , i = 1 , , N 1 ,
then ( A j ) 1 > 0 and B j 0 . Now, we show that ( d j + 1 d j ) is positive:
d j + 1 d j = d 1 j U 0 j + 1 f N 1 j U N j + 1 a 1 j U 0 j c N 1 j U N j
The first boundary is zero, therefore it is to show that f N 1 j > 0 and c N 1 j < 0 . As ξ 1 2 , then f N 1 j 0 . On the other hand, according to (33), c N 1 j is negative, and thus the theorem is proved. □
Theorem 4.
The proposed method is conditionally stable and convergent. Furthermore, the local truncation error is an O ( h 2 , Δ t ) .
Proof. 
Under condition (24), A j = [ A i r j ] is similar to a symmetric tridiagonal matrix, so the eigenvalues of A j , λ i ( A j ) , i = 1 , , N , are real [34]. Furthermore, A j is row diagonally dominant with
δ i = | A i i j | r i | A i r j | = 1 Δ t + 2 ξ r j , i = 1 , , N 1 ,
which yields ( A j ) 1 max 1 δ i (see in [34], Proposition 1.3). Therefore, we have
( A j ) 1 1 1 Δ t + 2 ξ r j ,
and by combining with B j = 1 Δ t + 2 ξ r j r j , we obtain
ρ ( ( A j ) 1 B j ) ( A j ) 1 B j ( A j ) 1 B j 1 Δ t + 2 ξ r j r j 1 Δ t + 2 ξ r j < 1 ,
where ρ ( ( A j ) 1 B j ) is the spectral radius of ( A j ) 1 B j . Therefore, the method is stable (see [35], p. 54) and then, via the Lax-theorem, convergent. The local truncation error is given by
T i j = u ( x i , t j + 1 ) u ( x i , t j ) Δ t ( σ i j ) 2 x i 2 2 u ( x i 1 , t j ) 2 u ( x i , t j ) + u ( x i + 1 , t j ) h 2 r j x i u ( x i + 1 , t j ) u ( x i 1 , t j ) 2 h + r j ξ ( u ( x i + 1 , t j ) + u ( x i 1 , t j ) ) + ( 1 2 ξ ) ( u ( x i + 1 , t j + 1 ) + u ( x i 1 , t j + 1 ) ) .
By Taylor’s expansion, we have
u ( x i , t j + 1 ) = u i j + Δ t u t i j + Δ t 2 2 ! 2 u t 2 i j + , u ( x i + 1 , t j ) = u i j + h u x i j + h 2 2 ! 2 u x 2 i j + h 3 3 ! 3 u x 3 i j + h 4 4 ! 4 u x 4 i j + , u ( x i 1 , t j ) = u i j h u x i j + h 2 2 ! 2 u x 2 i j h 3 3 ! 3 u x 3 i j + h 4 4 ! 4 u x 4 i j + , u ( x i + 1 , t j + 1 ) = u i j + h u x i j + Δ t u t i j + h 2 2 ! 2 u x 2 i j + Δ t 2 2 ! 2 u t 2 i j + h Δ t 2 u x t i j + , u ( x i 1 , t j + 1 ) = u i j h u x i j + Δ t u t i j + h 2 2 ! 2 u x 2 i j + Δ t 2 2 ! 2 u t 2 i j h Δ t 2 u x t i j + ,
and by substituting these expressions in (38), we obtain
T i j = u t 1 2 x 2 σ 2 ( x , t ) 2 u x 2 x r ( t ) u x + r ( t ) u i j + r j Δ t u t i j + r j 2 h 2 2 u x 2 i j 1 2 Δ t 2 u t 2 i j + r j 2 Δ t 2 2 u t 2 i j .
As u is the solution of Equation (9), then
u t 1 2 x 2 σ 2 ( x , t ) 2 u x 2 x r ( t ) u x + r ( t ) u i j = 0 .
Therefore, the principle term of the local truncation error is
r j Δ t u t i j + r j 2 h 2 2 u x 2 i j 1 2 Δ t 2 u t 2 i j + r j 2 Δ t 2 2 u t 2 i j ,
which completes the proof. □

3.3. Numerical Results

In this subsection, we give some numerical examples to confirm the theoretical results. In the last two examples, the results of our method are compared with the results obtained in [2]. The values of volatility and interest rate parameters have been taken from the works in [2,22,23,25,27]. These values are given in Table 1 with references. The other parameters for calculating the option price of the (9)–(12) model to all examples are equal T = 1 , K = 25 , S max = 100 and ε = 10 4 . As in [2], we have used the approximate solution obtained for N = 2 11 and M = 2 10 as a reference solution on the mesh points. We compute the maximum error as follows:
E N , M = max i , j | U i j N , M U ¯ ( x i , t j ) | ,
where U ¯ ( x , t ) is a linear interpolation of the approximate solution U 2 11 , 2 10 to get solutions at other points. Furthermore, we have used the following formula to compute the convergence rate:
R ¯ N , M = log E N , M E 2 N , M log 2 .
The error estimates along with the convergence rates, obtained with the proposed NSFD method, are listed in Table 2, Table 3, Table 4, Table 5, Table 6, Table 7 and Table 8. From Table 2, Table 3, Table 4, Table 5, Table 6, Table 7 and Table 8, we see that the convergence order in space is almost 2, which confirms the estimate given in Theorem 4.
It can be seen from Table 6 and Table 7 the good performance of the proposed method, concerning accuracy and convergence rates. Figure 1 and Figure 2 show some computed solutions for N = M = 64 and ε = 10 4 , showing the positivity preserving behavior.

4. Conclusions

In this paper, the generalized Black–Scholes equation for European option pricing is studied. To numerically solve this equation, a non-standard finite difference method has been used. In the proposed method, a non-local expression to approximate the reaction sentence of the generalized Black–Scholes and an implicit time step technique is used.
The theory results show that the proposed scheme is positivity preserving, stable and the order of the scheme with respect to the space variable is two. To confirm the theory results and good performance of the new method the option price for different values of volatility and interest rates is calculated. The numerical results confirm the theoretical behavior regarding the order of convergence for Theorem 4 (see Table 2, Table 3, Table 4, Table 5, Table 6, Table 7 and Table 8). For the last two examples, the option price values are plotted for all time levels, which shows that the new method is positive and non-oscillation (see Figure 1 and Figure 2). Furthermore, the obtained results are more accurate than other existing results in literature.
The results obtained indicate that non-standard difference schemes may be promising for solving problems which may have an impact on the stock price such as nonlinear Black–Scholes equations. In our subsequent investigation, we intend to address these problems using non-standard methods.

Author Contributions

Methodology, M.M.K., M.M.R., A.S., H.R. and P.K.; funding acquisition, M.M.R. All authors planned the scheme, developed the mathematical modeling, and examined the theory validation. The manuscript was written through the contribution by all authors. All authors discussed the results, reviewed, and approved the final version of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We are grateful for valuable comments from anonymous reviewers and the editor of this journal.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Black, F.; Scholes, M.S. The pricing of options and corporate liabilities. J. Political Econ. 1973, 81, 637–654. [Google Scholar] [CrossRef] [Green Version]
  2. Cen, Z.; Le, A. A robust and accurate finite difference method for a generalized Black-Scholes equation. J. Comput. Appl. Math. 2011, 235, 3728–3733. [Google Scholar] [CrossRef] [Green Version]
  3. Cho, C.; Kim, T.; Kwon, Y. Estimation of local volatilities in a generalized Black-Scholes model. Appl. Math. Comput. 2005, 162, 1135–1149. [Google Scholar] [CrossRef]
  4. Ladyzenskaja, O.A.; Solonnikov, V.A.; Ural’ceva, N.N. Linear and quasilinear equations of parabolic type. Translations of mathematical monographs. Am. Math. Soc. 1968, 23. [Google Scholar] [CrossRef]
  5. Kangro, R.; Nicolaides, R. Far field boundary conditions for Black-Scholes equations. SIAM J. Numer. Anal. 2000, 38, 1357–1368. [Google Scholar] [CrossRef]
  6. Ashyralyev, V.; Agirseven, D.; Agarwal, R.P. Stability estimates for delay parabolic differential and difference equations. Appl. Comput. Math. 2020, 19, 175–204. [Google Scholar]
  7. Ashyralyev, A.; Erdogan, A.S.; Tekalan, S.N. An investigation on finite difference method for the first order partial differential equation with the nonlocal boundary condition. Appl. Comput. Math. 2019, 18, 247–260. [Google Scholar]
  8. Harikrishnan, S.; Kanagarajan, K.; Elsayed, E.M. Existence and stability results for differential equations with complex order involving Hilfer fractional derivative. TWMS J. Pure Appl. Math. 2019, 10, 94–101. [Google Scholar]
  9. Qalandarov, A.A.; Khaldjigitov, A.A. Mathematical and numerical modeling of the coupled dynamic thermoelastic problems for isotropic bodies. TWMS J. Pure Appl. Math. 2020, 11, 119–126. [Google Scholar]
  10. Tunc, T.; Sarikaya, M.Z.; Yaldiz, H. Fractional hermite Hadamards type inquality for the co-ordinated convex functions. TWMS J. Pure Appl. Math. 2020, 11, 3–29. [Google Scholar]
  11. Mickens, R.E. Exact solutions to a finite difference model of a nonlinear reaction-advection equation: Implications for numerical analysis. Numer. Methods Part. Differ. Equ. 1989, 5, 313–325. [Google Scholar] [CrossRef]
  12. Mickens, R.E. Nonstandard Finite Difference Models of Differential Equations; World Scientific: Singapore, 1994. [Google Scholar]
  13. Anguelov, R.; Lubuma, J.M. Contributions to the mathematics of the nonstandard finite difference method and applications. Numer. Methods Part. Differ. Equ. 2001, 17, 518–543. [Google Scholar] [CrossRef]
  14. Baleanu, D.; Mehdizadeh Khalsaraei, M.; Shokri, A.; Kaveh, K. On the boundedness stepsizes-coefficients of A-BDF methods. AIMS Math. 2022, 7, 1562–1579. [Google Scholar] [CrossRef]
  15. Mehdizadeh Khalsaraei, M.; Shokri, A.; Ramos, H.; Heydari, S. A positive and elementary stable nonstandard explicit scheme for a mathematical model of the influenza disease. Math. Comput. Simul. 2020, 182, 397–410. [Google Scholar] [CrossRef]
  16. Mehdizadeh Khalsaraei, M.; Shokri, A.; Mohammadnia, Z.; Sedighi, H.M. Qualitatively Stable Nonstandard Finite Difference Scheme for Numerical Solution of the Nonlinear Black–Scholes Equation. J. Math. 2021, 2012, 6679484. [Google Scholar] [CrossRef]
  17. Mickens, R.E. A nonstandard finite difference scheme for a Fisher PDE having nonlinear diffusion. Comput. Math. Appl. 2003, 45, 429–436. [Google Scholar] [CrossRef] [Green Version]
  18. Mickens, R.E. Nonstandard finite difference schemes for reaction diffusion equations having linear advection. Numer. Methods Part. Differ. Equ. 2000, 4, 361–364. [Google Scholar] [CrossRef]
  19. Mickens, R.E.; Jordan, P.M. A new positivity-preserving nonstandard finite difference scheme for the DWE. Numer. Methods Part. Differ. Equ. 2005, 21, 976–985. [Google Scholar] [CrossRef]
  20. Mickens, R.E.; Jordan, P.M. A positivity-preserving nonstandard finite difference scheme for the Damped Wave Equation. Numer. Methods Part. Differ. Equ. 2004, 20, 639–649. [Google Scholar] [CrossRef]
  21. Shokri, A.; Mehdizadeh Khalsaraei, M.; Molayi, M. Nonstandard Dynamically Consistent Numerical Methods for MSEIR Model. J. Appl. Comput. Mech. 2022, 8, 196–205. [Google Scholar]
  22. Cen, Z.; Le, A.; Xu, A. Exponential time integration and second-order difference scheme for a generalized Black-Scholes equation. J. Appl. Math. Art. 2012, 2012, 796814. [Google Scholar] [CrossRef]
  23. Kadalbajoo, M.K.; Tripathi, L.P.; Kumar, A. A cubic B-spline collocation method for a numerical solution of the generalized Black-Scholes equation. Math. Comput. Model. 2012, 55, 1483–1505. [Google Scholar] [CrossRef]
  24. Valkov, R. Fitted finite volume method for a generalized Black-Scholes equation transformed on finite interval. Numer. Algorithms 2014, 65, 195–220. [Google Scholar] [CrossRef] [Green Version]
  25. Huang, J.; Cen, Z. Cubic spline method for a generalized Black-Scholes equation. Math. Probl. Eng. Art. 2014, 2014, 484362. [Google Scholar] [CrossRef]
  26. Mohammadi, R. Quintic B-spline collocation approach for solving generalized Black-Scholes equation governing option pricing. Comput. Math. Appl. 2015, 69, 777–794. [Google Scholar] [CrossRef]
  27. Chandra Sekhara Rao, S.; Manisha, S. High-order numerical method for generalized Black-Scholes model. Procedia Comput. Sci. 2016, 80, 1765–1776. [Google Scholar]
  28. Clevenhaus, A.; Ehrhardt, M.; Günther, M.; Ševcovic, D. Pricing American Options with a Non-Constant Penalty Parameter. J. Risk Financ. Manag. 2020, 13, 124. [Google Scholar] [CrossRef]
  29. Boire, F.M.; Mark Reesor, R.; Stentoft, L. American Option Pricing with Importance Sampling and Shifted Regressions. J. Risk Financ. Manag. 2021, 14, 340. [Google Scholar] [CrossRef]
  30. Sim, J. The economic and environmental values of the R and D investment in a renewable energy sector in South Korea. J. Clean. Prod. 2018, 189, 297–306. [Google Scholar] [CrossRef]
  31. Larissa, B.; Maran, R.M.; Ioan, B.; Anca, N.; Mircea-Iosif, R.; Horia, T.; Gheorghe, F.; Speranta, M.E.; Dan, M.I. Adjusted Net Savings of CEE and Baltic Nations in the Context of Sustainable Economic Growth: A Panel Data Analysis. J. Risk Financ. Manag. 2020, 13, 234. [Google Scholar] [CrossRef]
  32. Agaton, C.B. Application of real options in carbon capture and storage literature: Valuation techniques and research hotspots. Sci. Total Environ. 2021, 795, 148683. [Google Scholar] [CrossRef] [PubMed]
  33. Sisodia, G.S.; Alshamsi, R.; Sergi, B.S. Business valuation strategy for new hydroponic farm development—A proposal towards sustainable agriculture development in United Arab Emirates. Br. Food J. 2021, 123, 1560–1577. [Google Scholar] [CrossRef]
  34. Windish, G. M-Matrices in Numerical Analysis. In Teubner-Texte Zur Mathematik; Vieweg+Teubner Verlag: Wiesbaden, Germany, 1989; Volume 115. [Google Scholar]
  35. Smith, G.D. Numerical Solution of Partial Differential Equation: Finite Difference Method; Clarendon Press: Oxford, UK, 1985. [Google Scholar]
Figure 1. Computed values of U ( x , t ) for Example 6 with N = M = 64 and ε = 10 4 .
Figure 1. Computed values of U ( x , t ) for Example 6 with N = M = 64 and ε = 10 4 .
Symmetry 14 00141 g001
Figure 2. Computed values of U ( x , t ) for Example 7 with N = M = 64 and ε = 10 4 .
Figure 2. Computed values of U ( x , t ) for Example 7 with N = M = 64 and ε = 10 4 .
Symmetry 14 00141 g002
Table 1. Values of the volatility and interest rate parameters.
Table 1. Values of the volatility and interest rate parameters.
σ r
Example 1 0.2 + 0.2 ( 1 t ) ( x 25 1.2 ) 2 / ( ( x 25 ) 2 + 1.44 )  [22] 0.04  [27]
Example 2 0.2 1 + 0.1 ( 1 t ) ( x 1 + x )  [22] 0.04  [27]
Example 3 0.4 ( 2 + ( T t ) sin ( x ) )  [27] 0.06  [2]
Example 4 0.4 ( 2 + t sin ( x ) )  [23] 0.06  [2]
Example 5 0.4 ( 1 + t e x )  [23] 0.1 + 0.02 sin ( 10 T t )  [24]
Example 6 0.2 ( 1 + t e x )  [2] 0.06  [2]
Example 7 0.4 ( 2 + sin x )  [2] 0.06  [2]
Table 2. Errors and rates of convergence obtained by NSFD method for Example 1.
Table 2. Errors and rates of convergence obtained by NSFD method for Example 1.
M = 1024N = 32N = 64N = 128N = 256N = 512N = 1024
E N , M 1.54 × 10 1 7.43 × 10 2 3.54 × 10 2 1.53 × 10 2 5.76 × 10 3 1.36 × 10 3
R N , M 1.0561.0681.2071.4152.083-
Table 3. Errors and rates of convergence obtained by NSFD method for Example 2.
Table 3. Errors and rates of convergence obtained by NSFD method for Example 2.
M = 1024N = 32N = 64N = 128N = 256N = 512N = 1024
E N , M 1.54 × 10 1 7.43 × 10 2 3.54 × 10 2 1.53 × 10 2 5.76 × 10 3 1.36 × 10 3
R N , M 1.0561.0681.2071.4152.083-
Table 4. Errors and rates of convergence obtained by NSFD method for Example 3.
Table 4. Errors and rates of convergence obtained by NSFD method for Example 3.
M = 1024N = 32N = 64N = 128N = 256N = 512N = 1024
E N , M 4.39 × 10 1 5.75 × 10 2 2.27 × 10 2 7.61 × 10 3 1.76 × 10 3 3.51 × 10 4
R N , M 2.9351.3381.5792.1072.332-
Table 5. Errors and rates of convergence obtained by NSFD method for Example 4.
Table 5. Errors and rates of convergence obtained by NSFD method for Example 4.
M = 1024N = 32N = 64N = 128N = 256N = 512N = 1024
E N , M 4.75 × 10 1 8.48 × 10 2 3.29 × 10 2 9.65 × 10 3 2.67 × 10 3 5.38 × 10 4
R N , M 2.4851.3661.7701.8502.314-
Table 6. Errors and rates of convergence obtained by NSFD method for Example 5.
Table 6. Errors and rates of convergence obtained by NSFD method for Example 5.
M = 1024N = 32N = 64N = 128N = 256N = 512N = 1024
E N , M 1.90 × 10 1 9.26 × 10 2 4.23 × 10 2 1.73 × 10 2 5.15 × 10 3 1.10 × 10 3
R N , M 1.0421.1311.2891.7492.227-
Table 7. Comparison of errors and rates of convergence obtained by NSFD method and the method by Cen and Le [2] for Example 6.
Table 7. Comparison of errors and rates of convergence obtained by NSFD method and the method by Cen and Le [2] for Example 6.
M = 1024 N = 128N = 256N = 512N = 1024
Results for NSFD method E N , M 3.33 × 10 2 1.44 × 10 2 5.30 × 10 3 1.25 × 10 3
R N , M 1.2041.4492.082-
Results for Cen and Le [2] E N , M 6.80 × 10 2 3.12 × 10 2 1.25 × 10 2 3.41 × 10 3
R N , M 1.1221.3121.880-
Table 8. Comparison of errors and rates of convergence obtained by NSFD method and the method by Cen and Le [2] for Example 7.
Table 8. Comparison of errors and rates of convergence obtained by NSFD method and the method by Cen and Le [2] for Example 7.
M = 1024 N = 128N = 256N = 512N = 1024
Results NSFD method E N , M 3.53 × 10 2 1.03 × 10 2 2.71 × 10 3 5.38 × 10 4
R N , M 1.7731.9332.331-
Results for Cen and Le [2] E N , M 5.71 × 10 2 1.94 × 10 2 5.15 × 10 3 1.05 × 10 3
R N , M 1.5511.9192.284-
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mehdizadeh Khalsaraei, M.; Rashidi, M.M.; Shokri, A.; Ramos, H.; Khakzad, P. A Nonstandard Finite Difference Method for a Generalized Black–Scholes Equation. Symmetry 2022, 14, 141. https://doi.org/10.3390/sym14010141

AMA Style

Mehdizadeh Khalsaraei M, Rashidi MM, Shokri A, Ramos H, Khakzad P. A Nonstandard Finite Difference Method for a Generalized Black–Scholes Equation. Symmetry. 2022; 14(1):141. https://doi.org/10.3390/sym14010141

Chicago/Turabian Style

Mehdizadeh Khalsaraei, Mohammad, Mohammad Mehdi Rashidi, Ali Shokri, Higinio Ramos, and Pari Khakzad. 2022. "A Nonstandard Finite Difference Method for a Generalized Black–Scholes Equation" Symmetry 14, no. 1: 141. https://doi.org/10.3390/sym14010141

APA Style

Mehdizadeh Khalsaraei, M., Rashidi, M. M., Shokri, A., Ramos, H., & Khakzad, P. (2022). A Nonstandard Finite Difference Method for a Generalized Black–Scholes Equation. Symmetry, 14(1), 141. https://doi.org/10.3390/sym14010141

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop