DDNSDiscrete Dynamics in Nature and Society1607-887X1026-0226Hindawi10.1155/2020/80473478047347Research ArticleNovel Approaches for Getting the Solution of the Fractional Black–Scholes Equation Described by Mittag-Leffler Fractional Derivativehttps://orcid.org/0000-0002-8664-6464SeneNdolane1SèneBabacar2NdiayeSeydou Nourou3TraoréAwa4Al-MdallalQasem M.1Laboratoire LmdanDépartement de Mathématiques de la DécisionUniversité Cheikh Anta Diop de DakarFaculté des Sciences Economiques et GestionBP 5683 Dakar FannSenegalucad.sn2Centre de Recherches Economiques AppliquéesDirecteur du Laboratoire de Finances pour le Developpement (LAFIDEV)Faculté des Sciences Economiques et GestionUniversité Cheikh Anta Diop de DakarBP 5683 Dakar FannSenegalucad.sn3Centre de Recherches Economiques AppliquéesLaboratoire de Finances pour le Development (LAFIDEV)Faculté des Sciences Economiques et GestionUniversité Cheikh Anta Diop de DakarBP 5683 Dakar FannSenegalucad.sn4Centre de Recherches Economiques AppliqúeesFacultédes Sciences Economiques et GestionUniversité Cheikh Anta Diop de DakarBP 5683 Dakar FannSenegalucad.sn202025720202020060520202705202025720202020Copyright © 2020 Ndolane Sene 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.

The value of an option plays an important role in finance. In this paper, we use the Black–Scholes equation, which is described by the nonsingular fractional-order derivative, to determine the value of an option. We propose both a numerical scheme and an analytical solution. Recent studies in fractional calculus have included new fractional derivatives with exponential kernels and Mittag-Leffler kernels. These derivatives have been found to be applicable in many real-world problems. As fractional derivatives without nonsingular kernels, we use a Caputo–Fabrizio fractional derivative and a Mittag-Leffler fractional derivative. Furthermore, we use the Adams–Bashforth numerical scheme and fractional integration to obtain the numerical scheme and the analytical solution, and we provide graphical representations to illustrate these methods. The graphical representations prove that the Adams–Bashforth approach is helpful in getting the approximate solution for the fractional Black–Scholes equation. Finally, we investigate the volatility of the proposed model and discuss the use of the model in finance. We mainly notice in our results that the fractional-order derivative plays a regulator role in the diffusion process of the Black–Scholes equation.

1. Introduction

There are many mathematical models  used in finance to predict the values of cost, revenue, and options. In this paper, we address the application of fractional calculus in economics and finance. Fractional derivatives occupy an important place in fractional calculus, so this paper investigates the use of fractional derivatives for modeling financial and economic models. The Black–Scholes model is an important tool used in finance to predict the value of an option . There are many styles of options: European options [3, 4], American options [5, 6], and Asian options. The pricing of options is a subject that has been very intensely debated in economics. So far, the economic and financial literature continues to take an interest in this subject. Today, the questions arise concerning alternative methods for pricing European options, which are derivatives that can only be exercised at maturity. The literature is also interested in the valuation of American options (there is not yet a consensus model). American options are derivatives that can be exercised at any time. The methodology proposed in this paper, therefore, makes it possible to price a European option. It could be extended to American and Asian options and knock-in and knock-out barrier options. Another possible application is the determination of the default risk of bonds listed on the financial market. The model could possibly be applied to determine the unknowns of the Merton–Black–Scholes-based model. Note that predicting the exact value of the American option is an ongoing problem that is generally not common knowledge. For European options, the model proposed by Black and Scholes  has an exact analytical solution and a numerical scheme and that analytical solution uses the normal distribution. The use of the normal distribution in the formula for an option is not always suitable; thus, our work proposes an analytical solution to this problem using recursive approximation that avoids the need for a normal distribution in the formula for the value of a European option.

Some studies in the context of fractional calculus have investigated the use of the Black and Scholes equation for European options. Much of this research has involved the use of the Caputo–Liouville and the Riemann–Liouville derivatives. In , Fall et al. offered a new work on fractional Black–Scholes equations described by the generalized fractional derivative. In , the authors present a new procedure to obtain the analytical solutions of the fractional Black–Scholes model, which is called the homotopy perturbation method. Sawangston et al. proposed an analytical solution for the fractional Black–Scholes equation with the Caputo–Liouville derivative in . The research in  offered an analytical solution for the fractional Black–Scholes described by the conformable derivative. For more numerical schemes proposed for the Black–Scholes equation, see [1, 6, 913]. For other numerical procedures in other differential equations, the readers can refer to .

Modeling the physical phenomena using the fractional-order derivative has many advantages. First of all, it permits us to study the differential models with arbitrary noninteger-order derivatives. Second, it permits us to take into account the memory effect; that is, the next behavior of the dynamic is explained by the past behavior of the dynamic. As we will observe in this paper, the fractional-order derivative can play a regulator role in differential dynamics. In this paper, we investigate and introduce the fractional Black–Scholes equation described by a new fractional-order derivative: the Mittag-Leffler fractional derivative . We prove this can be used in finance to evaluate the value of an option. Note that the new model uses the Atangana–Baleanu derivative and takes into account the integer-order time derivative. We propose a numerical scheme of the introduced model using the Adams–Bashforth method and an analytical solution using the recursive procedure proposed by Liao in . Finally, we use the analytical solution to analyse the volatility of the option to a change in the price of the underlying security.

This manuscript is structured as follows: In Section 2, we define the fractional derivatives with Mittag-Leffler kernels. In Section 3, we introduce and discuss the fractional Black–Scholes equation described by the Mittag-Leffler fractional derivative. In Section 4, we prove the existence and uniqueness of our introduced model. In Section 4, we describe the Adams–Bashforth method. In Section 5, we propose the numerical discretization of the fractional Black–Scholes equation. In Section 6, we propose an analytical solution for the Black–Scholes equation described by the Atangana–Baleanu fractional derivative. In Section 7, we analyse the volatility of the fractional Black–Scholes equation, and in Section 8, we present graphical representations and a discussion. Section 9, includes our final remarks and conclusions.

2. Fractional Operators without Singular Kernels

In this section, we recall the definitions of fractional derivatives with nonsingular kernels. Fractional calculus began with the Riemann–Liouville fractional derivative and the Caputo–Liouville fractional derivative; for definitions, refer to . These two derivatives were proposed in response to Leibniz’s question in 1695. Many other fractional derivatives were introduced in the literature after this date and have been applied in physics, mechanics, and mathematical modeling . Some properties were lost in modeling real-world phenomena, specifically the physical aspects. Recently, Caputo and Fabrizio proposed a new fractional derivative with an exponential kernel  and gave it an associated fractional integral. The main advantage of this fractional derivative is its lack of singularity, and it can be used in modeling many physical phenomena [26, 33]. The definitions are given in the following.

Definition 1 (see [<xref ref-type="bibr" rid="B32">32</xref>]).

The Caputo–Fabrizio fractional derivative of the function u:×0,+ of order α is defined in the form(1)DαCFuy,t=Mα1α0tuys,sexpα1αtsds,where t>0 and the order α0,1, with the normalization term satisfying M0=M1=1.

Definition 2 (see [<xref ref-type="bibr" rid="B32">32</xref>]).

The Caputo–Fabrizio integral for a given function u:×0,+, of order α0,1, is defined in the form(2)IαCFuy,t=21α2αMαuy,t+2α2αMα0tuys,sds,for all t>0 and the order α0,1, with the normalization term satisfying M0=M1=1.

The exponential form is a particular case of the Mittag-Leffler function. Motivated by the fact that the Cauchy problem with the Caputo–Fabrizio derivative generates a solution with an exponential function, Atangana and Baleanu proposed another fractional derivative with a Mittag-Leffler kernel in 2016 . In other words, we can translate the Cauchy equation using Caputo–Fabrizio as the first-order equation with an integer-order derivative.

Definition 3 (see [<xref ref-type="bibr" rid="B21">21</xref>]).

The Atangana–Baleanu–Caputo derivative for a function u:×0,+, of order α0,1, is defined in the form(3)DαABCuy,t=Bα1α0tuys,sEαα1αtsαds,for all t>0, where the function Γ is the Euler gamma function and Eα denotes the Mittag-Leffler function, with the normalization term satisfying B0=B1=1.

Definition 4 (see [<xref ref-type="bibr" rid="B21">21</xref>]).

The Atangana–Baleanu integral for a given function u:×0,+, of order α0,1, is defined as the form(4)IαABuy,t=1αBαuy,t+αBα1Γα0ttsα1uys,sds,for all t>0, where the function Γ is the Euler gamma function and Eα denotes the Mittag-Leffler function, with the normalization term satisfying B0=B1=1.

We were motivated to consider these two fractional derivatives because of their successful application in modeling real-life phenomena. In this paper, we apply the Atangana–Baleanu fractional derivative in modeling the value of options and investigate the fractional Black–Scholes equation described by the Atangana–Baleanu–Caputo fractional derivative.

3. Black–Scholes in the Context of Mittag-Leffler Fractional Derivative

In this section, we introduce the Black–Scholes equation in the context of the Atangana–Baleanu fractional derivative. We begin by recalling the classical model proposed by Black and Scholes . Let the asset price be S at time, the constant volatility of an underlining asset be represented by the parameter σ, and μ be the expected rate of return. Myron and Fischer stipulated that the stock price follows a Brownian motion denoted by the parameter w; thus, we have the following:(5)dS=μSdt+σSdw.

Note that equation (5) describes the asset price S as Brownian motion and represents a particular case of Ito’s lemma. There are two types of derivations in the literature for pricing options: in this paper, we use Fisher and Myron’s derivation. Black and Scholes expressed the value of the portfolio denoted by P in the following form:(6)dP=dVVSdS,where V represents the value of an option. We have decomposed the portfolio into the value of the option and the asset price. V/SdS denotes the asset price obtained per year. Using Ito’s lemma, equation (6) can be represented in the following form:(7)dP=Vτ+μSVS+σ22S2V2S2VSSμdt+σSVSσVSSdw.

Taking into account the interest rate r into the value of the portfolio, we express equation (6) in the following form:(8)dP=rVdtrVSSdt.

Combining equations (7) and (8), the differential equation which calculates the value of a European option is given by the following equation:(9)Vt+σ22S22VS2+rSVSrV=0.

We made the following assumptions related to equation (9): it considers a European option, the risk-less interest rate r is constant, and there are no transaction costs, and we authorize the possibility to buy and to sell any number of stocks with no restriction to short selling at the last moment. The boundary conditions for equation (9) are defined as V0,t=0, VS,TS as S, and the terminal condition is given by(10)VS,T=maxSE,0,where the parameter E denotes the strike price of the underlying stock and T represents the expiration time. For the European option, we have the possibility to buy and to exercise the option with no obligation at time T. That means we can sell the risky asset to a seller at a strike price E. Equation (10) can be explained as follows: we exercise the option at time T when the condition E<S is held. That is, VS,T=maxSE,0=SE; in other words, the buyer receives the playoff SE. The benefit is in selling the asset to the seller of the contract rather than on the financial market. When the condition S<E is held, the contract is not good for the buyer, and the buyer can sell the risky asset for a larger price on the financial market.

The use of equation (9) for the analytical solution or the numerical scheme is not trivial, but equation (9) is a diffusion equation and can be rewritten more simply. We use the following changes to the variables described by the relationships:(11)S=Eex,t=T2τσ2,V=Eux,t.

From which it follows the classical Euler equation given by(12)uτ=2ux2+k1uxku,with the initial boundary condition defined by(13)ux,0=maxex1,0,where k=2r/σ2 denotes the balance between the free interest rate and the volatility of the stocks. The difference between equation (9) and equation (12) shows the transition between the finance model and the physical model. Note that it was provided in fractional calculus; many models described by the classical derivative as in equation (9) cannot describe the real behavior of the modeled phenomena. It has been proved that many real-world problems follow fractional phenomena. Equation (12) is a diffusion equation, and there are many diffusion processes in physics, such as subdiffusion, superdiffusion, ballistic diffusion, and hyperdiffusion, which equation (12) does not take into account. This problem requires the introduction of a fractional derivative that takes into account all types of diffusion processes. In this paper, we replace the ordinary time derivative with the fractional-order time derivative, and we describe the fractional differential equation which we consider using the following fractional differential equation:(14)DαABCu=2ux2+k1uxku,with the initial boundary condition defined by(15)ux,0=maxex1,0.

In the next section, we try to prove our new model is well defined, admit a unique solution, and use numerical and analytical methods to approach it. For the readers and more understanding of the paper, we summarize the description of the parameters used in this paper in Table 1.

Parameters of the Black–Scholes equation.

ParametersDescription of the parameters.
V and TThe value of an option and the expiration time, respectively.
σThe volatility of the underlining stock.
kThe balance between the free interest rate and the volatility of the stock.
rThe risk-less interest rate.
E and SThe strike price of the underlying stock and the asset price, respectively.
4. Adams–Bashforth Numerical Approach

In this section, we describe the procedure of discretization used in this paper. The method is called the Adams–Bashforth numerical scheme and was introduced in fractional calculus by Atangana in . Adams–Bashforth is a useful method that involves the following fundamental theorem .

Theorem 1 (see [<xref ref-type="bibr" rid="B34">34</xref>]).

The solution of the fractional differential equation described by DαABCv=fv,t with initial boundary condition v0 satisfies the following relationship:(16)vtv0=1αBαfv,t+αBα0ttsα1fv,sds.

The approximation of the function f uses a Lagrange polynomial, which is the main novelty of Atangana’s proposed numerical approximation. The following relationship describes the Lagrange polynomial:(17)pt=ttn1tntn1fvn,tn+ttn1tn1tnfvn1,tn1.

Using equation (17), the discretized approximation of equation (16) at time tn+1 and tn, considering h=tntn1 as described in , is obtained by the following equation:(18)vtn+1vtn=θα,1+θα,2,where the function θα,1 is given by the expression described in the following equation:(19)θα,1=fvn,tn1αBα+αBαΓαh2htn+1α+1αtn+1α+1α+1αBαΓαhhtnααtα+1α+1,and where the function θα,2 is given by the expression described in the following equation:(20)θα,2=fvn1,tn1α1BααBαΓαhhtn+1ααtn+1α+11+α+tα+1BαΓαh.

The above discretization proposed by Atangana in  is very useful in numerical solutions for fractional differential equations described by certain fractional derivatives, such as the Caputo fractional derivative and the fractional derivative with the Mittag-Leffler kernel. The application of the Adams–Bashforth scheme uses the discretization of the terms fvn,tn and fvn1,tn1. We can do the standard discretization procedures to discretize them. Here, we use the central difference schemes for the second-order space derivative and numerical approximation of the first-order space derivative. Before moving on, we look at the stability of the method used and see that stability is obtained when the function f is Lipschitzian. Using equations (19) and (20), the following relationship is obtained:(21)vtn+1vtn=1αBαfvn,tnfvn1,tn1+Ωα,1Ωα,2,where the function Ωα,1 is given by the relationship(22)Ωα,1=αBα0tn+1tn+1sα1fv,sds,and the function Ωα,2 is given by the relationship(23)Ωα,2=αBα0tntnsα1fv,sds.

We will find a threshold for the function Ωα,1Ωα,2 by applying the Euclidean norm:(24)Ωα,1Ωα,2αBα0tn+1tn+1sα1fv,sdsαBα0tntnsα1fv,sdsαBα0tn+1tn+1sα1fv,sds+αBα0tntnsα1fv,sdsαBα0tn+1tn+1sα1fv,sds+αBα0tntnsα1fv,sdsαMBα0tn+1tn+1sα1ds+αMBα0tntnsα1dsαMBαΓα+1tn+1α+tnα=αMhαBαΓα+1n+1α+nα.

Using the assumption tn=nh, we get the following relationship:(25)Ωα,1Ωα,2αMhαBαΓα+1n+1α+nα.

Applying the norm to both sides of equation (21), we obtain the following relationships:(26)vtn+1vtn1αBαfvn,tnfvn1,tn1+Ωα,1Ωα,2,1αBαfvn,tnfvn1,tn1+αMhαBαΓα+1n+1α+nα.

From equation (26), it can be seen that when the function f is locally Lipschitz and h converges to zero, we obtain the following relationship:(27)vtn+1vtn0.

We can conclude the Adams–Bashforth numerical scheme is unconditionally stable. In the next section, we apply the Adams–Bashforth numerical scheme to the numerical approximation of the fractional Black–Scholes equation described by the fractional derivative with Mittag-Leffler.

5. Numerical Approach for Fractional Black–Scholes Equation

In this section, we describe the Adams–Bashforth numerical scheme for the fractional Black–Scholes equation represented by the Atangana–Baleanu fractional derivative. Let us begin the numerical approximation of the fractional Black–Scholes equation. Let tn=nh. The Adams–Bashforth numerical scheme for the fractional Black–Scholes equation described by the Atangana–Baleanu fractional derivative takes the following form:(28)utn+1utn=θα,1θα,2,where the function θα,1 is given by the expression described in the following equation:(29)θα,1=fun,tn1αBααBαΓαhα2n+1ααn+1α+1α+1αBαΓαhαnααnα+1α+1,and where the function θα,2 is given by the expression described in the following equation:(30)θα,2=fun1,tn11αBααBαΓαhαn+1ααn+1α+1α+1+nα+1BαΓαh.

The next step consists of finding the discretization of the functions fvn,tn and fvn1,tn1. In the Black–Scholes equation, the function f is given by(31)fu,t=2ux2+k1uxku.

Using the central difference approximation for the second-order derivative with respect to the space coordinate and the numerical approximation for the space derivative, we obtain the following discretization at the points un and tn for the function f:(32)fun,tn=uj+1n2ujn+uj1nΔx2+k1uj+1nuj1n2Δxkujn+OΔx.

Using the central difference approximation again for the second-order derivative with respect to the space coordinate and the numerical approximation for the space derivative, we obtain the following discretization at the points un1 and tn1 for the function f:(33)fun1,tn1=uj+1n12ujn1+uj1n1Δx2+k1uj+1n1uj1n12Δxkujn1+OΔx.

Numerical discretization using the Adams–Bashforth method for the Black–Scholes equation is obtained by combining equations (28)–(33). That is,(34)ujn+1=ujn+θα,1θα,2.

For the computation of our numerical schemes, we make some changes to the variables, such that the terms depending on n, h, and α are constants. Thus, we rewrite θα,1 and θα,2, respectively, as in the following relationships:(35)Hn,h,α,1=1αBααBαΓαhα2n+1ααn+1α+1α+1αBαΓαhαnααnα+1α+1,(36)Hn,h,α,2=1αBααBαΓαhαn+1ααn+1α+1α+1+nα+1BαΓαh.

Finally, the numerical scheme using the Adams–Bashforth method for the fractional Black–Scholes equation is given by(37)ujn+1=ujn+Hn,h,α,1uj+1n2ujn+uj1nΔx2+k1uj+1nuj1n2ΔxkujnHn,h,α,2uj+1n12ujn1+uj1n1Δx2+k1uj+1n1uj1n12Δxkujn1.

To complete the numerical discretization given in equation (37), we recall the discretized form of the initial boundary condition described by the following expression:(38)u0nx=maxex1,0.

6. Analytical Solution for Fractional Black–Scholes Equation

In this section, we use the fractional integrator to propose the analytical solution of the fractional Black–Scholes equation described by the Atangana–Baleanu fractional derivative. The method is described in the following theorem.

Theorem 2.

The solution of the fractional differential equation described by DαABCv=fv,τ with initial boundary condition v0 satisfies the following relationship:(39)vn+1τvn+10=1αBαfvn,τ+αBα0τfs,vnds,where n=0,1,2, Furthermore, the solution is given by(40)vx,τ=v00+v1τ+v2τ+.

Consider the fractional Black–Scholes equation (14). Using Theorem 2, with assumption u1x,0=0, we have the following solution, which is the first step:(41)u1x,τ=u1x,0+1αBαkmaxex,0maxex1,0+kmaxex,0maxex1,0ταBαΓ1+α=1αBαkmaxex,0maxex1,0+αkmaxex,0maxex1,0ταBαΓ1+α=kmaxex,0maxex1,01αBα+αταBαΓ1+α.

Under assumption u2x,0=0, again using equation (39), the second step gives the following solution:(42)u2x,τ=k2maxex,0maxex1,01α1αBα+αταBαΓ1+αk2maxex,0maxex1,0α1αταBαΓ1+α+ατ2αBαΓ1+2α.

We adopt the procedure described in Theorem 2 for the rest of the steps. The following expression gives the approximate solution of the fractional Black–Scholes equation in the context of the Mittag-Leffler fractional derivative:(43)ux,τ=ux,0+u1x,0+u2x,0+.

We recover the approximate solution of the classical Black–Scholes equation when α=1. In the context of the Mittag-Leffler fractional derivative, when we suppose α=1, we get the following form:(44)ux,τ=ux,0+u1x,0+u2x,0+=maxex1,0+kmaxex,0maxex1,0τk2maxex,0maxex1,0τ2+=maxex1,01ekτ+maxex,0ekτ.

Due to space limitation, all term are not written. Finally, we describe the analytical solution for the classical Black–Scholes equation (12) in the following form:(45)ux,τ=maxex,01ekτ+maxex1,0ekτ.

The method adopted here gives a solution which is in good agreement with the classical solution of the Black–Scholes equation (12). Let us prove our solution (45) is another representation for the traditional solution of the Black–Scholes equation using the following:(46)VS,E=SNd1EekτNd2,where N designs the normal distribution function, d1=logS/E+r+σ2/2Tt/σTt, and d2=d1σTt. Note in the case “in the money,” that is, E<S, we have high volatility; in other words, we have the following equation:(47)Nd1=Nd2=1.

From which, we rewrite equation (46) as the following form:(48)VS,E=SEekτ.

Equation (48) represents the value of the option “in the money.” We recover this value with the solution proposed in equation (45) using the following reasoning. Note that in the case “in the money,” there exists the following relationship:(49)maxex,0=ex,maxex1,0=ex1.

Now, we replace ex=S/E into equation (45), and using the assumption posed in equation (11), we obtain the following relationship:(50)ux,τ=SE1ekτSE1ekτ=SEekτ.

Multiplying the function u by E as in equation (11), we obtain(51)VS,E=SEekτ.

We can see that the solution represented in equation (48) and the solution in equation (51) are the same. Thus, our analytical solution can be used in finance to determine the value of an option, satisfying the Black–Scholes equation.

Note that, in the money “call” is when the price of the underlying asset is higher than the strike price. It is in the interest of the holder of the option to exercise it. He/she has made good anticipations (S>E); in this case, the value of the option is given by the formula V=maxSE,0=SE. Out of the money “call” is when the price of the underlying asset is lower than the strike price. It is not in the interest of the option holder to exercise the option. He/she has made wrong expectations (S<E); in this case, the value of the option is given by V=maxSE,0=0. In the money “put” is when the price of the underlying asset is lower than the strike price. It is in the interest of the option holder to exercise the option. He/she has made good anticipations (S<E); in this case, the value of the option is given by V=maxES,0=ES. Out of the money “put” is when the price of the underlying asset is higher than the strike price. It is not in the interest of the option holder to exercise the option. He/she has made wrong anticipations (S>E); in this case, the value of the option is given by V=maxES,0=0.

7. The Volatility of the Fractional Black–Scholes Equation

In this section, we analyse the volatility of the fractional Black–Scholes equation using an analytical solution. In the market, we can buy the call and put with different strike prices and maturity. The volatility analyses the liquidity of the cost of the call and the put in the market. Many types of volatility can be generated by the Black–Scholes equation. In general, volatility measures risk in the financial markets. Volatility is used to control both upward and downward movements. It is calculated from log returns. The delta is part of the Greek letters of options. It is the derivative of the price of the call or put option in relation to the price of the underlying asset. It is used for trading, arbitrage, or hedging operations on options. It is an important indicator of market risk management. The Basel Committee on International Banking Regulation recommends that banks use delta for exposures to options. In reality delta is a sensitivity factor.

Here, we recall the formula of the volatility delta of the fractional Black–Scholes equation. The delta measures the sensitivity of the option price to a change in the price of the underlying security. Under variable changes, this is expressed as follows:(52)δ=ux.

Using the approximate solution (52), it is clear that the sensitivity of the option price to a change in the price of the underlying security does not depend on the fractional-order and is given by(53)δ=maxex,0.

Given the conditions in equation (11), the volatility of the option price to a change in the price of the underlying security is given by(54)δ=maxSE,0.

Earlier, we stated that the volatility studied in this paper does not depend on the order of the fractional derivative or the balance between the free interest rate and the volatility of the stocks k. This remark can be explained simply from the fact that, except for the first term ux,0, we have the following relationship:(55)uix=0,for all i=1,2,3,.. Figure 1 shows the volatility surface of the option price to a change in the price of the underlying security in Figure 1.

Volatility surface of the Black–Scholes equation.

The volatility surface gives the sensitivity of the option price when the asset price S and the strike price of the underlying stock E both vary in time. Let the asset price S be fixed, and let the strike price of the underlying stock E vary in time. The behavior of the volatility can be seen in Figure 2. In other words, we can see that when S<E, the volatility decreases rapidly and converges to zero. Let the underlying stock E and the asset price S vary in time. The behavior of the volatility can be seen in Figure 2. Thus, we can see in Figure 2 that when E<S, the volatility quickly increases linearly and converges to infinity. We note high volatility.

Volatility when S<E and E<S.

8. Graphical Representations and Discussion

In this section, we illustrate our results graphically. Specifically, we depict the behavior of the solutions of the fractional Black–Scholes equation obtained with numerical schemes and the recursive method. We begin with the approximate solution generated by the recursive method previously described in Section 7. In this context, we consider the solution with three iterations given by equation (43). Figure 3 shows the behavior of the approximate solution for the fractional Black–Scholes equation with the Atangana–Baleanu fractional derivative with order α=0.5 and the balance between the free interest rate and the volatility of the stocks k=2. The constant k is obtained with the risk-free interest rate to expiration r=0.04 and the volatility of the stocks σ=0.2 .

Behavior of the Black–Scholes equation with α=0.5.

Figure 4 shows the behavior of the approximate solution for the fractional Black–Scholes equation with the Atangana–Baleanu fractional derivative with order α=1 and the balance between the free interest rate and the volatility of the stocks k=2.

Behavior of the Black–Scholes equation with α=1.

Figure 5 depicts the analytical solution with the recursive method and the approximate solution obtained with the Adams–Bashforth numerical scheme. We suppose t=1 and α=0.5 and observed the numerical solution, and the approximate analytical solution is in good agreement.

Analytical solution vs numerical approximation with α=0.5.

The main question now is how to find the optimal order of α. In our previous example, we depicted the figures by choosing the order α=0.5 or α=1. A complicated method is required to find the optimal order α in fractional calculus. First, we must collect the data in the considered market. Second, we create a figure generated by the obtained data, and finally, we proceed by interpolation to find an order α.

To support the numerical discretization, we describe in Table 2 the different values of the options generated by our numerical scheme, and we compare them with values obtained with the homotopy perturbation method. In other words, the robustness of the used numerical scheme is authenticated by comparing the numerical and analytical results in Table 2. In Table 2, we consider α=0.5, the risk-free interest rate to expiration α=0.5, the volatility of the stocks σ=0.2, the volatility of the stocks k=2, and the strike price of the underlying stock E=10.

Option values when E<S.

xτuNumericaluHomotopyuclassic, α=1
0.15410.010.32400.32421.1467
0.46770.030.84550.84561.5381
0.24010.050.57450.57471.1762
0.16970.020.49240.49251.1457

We mainly observe that the results in Table 2 are in good agreement with the results on the Black–Scholes equation studied in terms of the Caputo derivative in . In general, we also notice that the order of the fractional derivative has a significant impact on the value of the options. In this investigation, the order α has a regulator impact, and we notice in Table 2 that the values obtained using the fractional order are more beneficial rather than the values obtained with the classical derivative in the financial market.

9. Conclusion

In this paper, we have discussed the numerical scheme and the analytical solution for the fractional Black–Scholes equation described by the Atangana–Baleanu derivative. As observed, the analytical solution of the fractional Black–Scholes equation with the Atangana–Baleanu fractional derivative is not trivial. We have used the Adams–Bashforth numerical scheme to approach the solution, as it is useful and straightforward for proposing the approximate solutions of the fractional Black–Scholes equation. We have also considered the liquidity of the cost of the call and the put in the market, namely, the volatility. The graphical representations have proved the good agreements between the analytical solution and the numerical solutions for the fractional Black–Scholes equation.

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

AkramiM. H.ErjaeeG. H.Numerical solutions for fractional black-scholes option pricing equationGlobal Analysis and Discrete Mathematics20151914BlackF.ScholesM.The pricing of options and corporate liabilitiesJournal of Political Economy197381363765410.1086/260062KumarS.KumarD.SinghJ.Numerical computation of fractional Black-Scholes equation arising in financial marketEgyptian Journal of Basic and Applied Sciences201413-417718310.1016/j.ejbas.2014.10.003FallA. N.NdiayeS. N.SeneN.Black-Scholes option pricing equations described by the Caputo generalized fractional derivativeChaos, Solitons & Fractals201912510811810.1016/j.chaos.2019.05.0242-s2.0-85066141016ChenW.DuK.QiuX.Analytic properties of american option prices under a modified Black-Scholes equation with spatial fractional derivativePhysica A: Statistical Mechanics and its Applications2017491374410.1016/j.physa.2017.08.0682-s2.0-85032023053SongL.WangW.Solution of the fractional black-scholes option pricing model by finite difference methodAbstract and Applied Analysis201320131019428610.1155/2013/1942862-s2.0-84880168354SawangtongP.TrachooK.SawangtongW.WiwattanapatapheeB.The analytical solution for the black-scholes equation with two assets in the liouville-caputo fractional derivative senseMathematics20186812910.3390/math60801292-s2.0-85052818915YavuzM.ÖzdemirN.A different approach to the European option pricing model with new fractional operatorMathematical Modelling of Natural Phenomena20181311210.1051/mmnp/20180092-s2.0-85045213243CenZ.LeA.A robust and accurate finite difference method for a generalized Black-Scholes equationJournal of Computational and Applied Mathematics2011235133728373310.1016/j.cam.2011.01.0182-s2.0-79955570032PhaochooP.LuadsongA.AschariyaphothaN.The meshless local Petrov-Galerkin based on moving kriging interpolation for solving fractional Black-Scholes modelJournal of King Saud University—Science201628111111710.1016/j.jksus.2015.08.0042-s2.0-84955185013YavuzM.ÖzdemirN.European vanilla option pricing model of fractional order without singular kernelFractal and Fractional201821310.3390/fractalfract2010003YavuzM.ÖzdemirN.A quantitative approach to fractional option pricing problems with decomposition seriesKonuralp Journal of Mathematics201861102109YavuzM.ÖzdemirN.OkurY. Y.Generalized differential transform method for fractional partial differential equation from financeProceedings of the International Conference on Fractional Differentiation and its Applications2016Novi Sad, Serbia778785WakifA.A novel numerical procedure for simulating steady MHD convective flows of radiative casson fluids over a horizontal stretching sheet with irregular geometry under the combined influence of temperature-dependent viscosity and thermal conductivityMathematical Problems in Engineering2020202020167535010.1155/2020/1675350WakifA.ChamkhaA.ThummaT.AnimasaunI. L.SehaquiR.Thermal radiation and surface roughness effects on the thermo-magneto-hydrodynamic stability of alumina-copper oxide hybrid nanofluids utilizing the generalized Buongiorno’s nanofluid modelJournal of Thermal Analysis and Calorimetry2020In press10.1007/s10973-020-09488-zWakifA.BoulahiaZ.MishraS. R.RashidiM. M.SehaquiR.Influence of a uniform transverse magnetic field on the thermo-hydrodynamic stability in water-based nanofluids with metallic nanoparticles using the generalized Buongiorno’s mathematical modelThe European Physical Journal Plus201813318110.1140/epjp/i2018-12037-72-s2.0-85046625848NayakM. K.WakifA.AnimasaunI. L.AlaouiM. S. H.Numerical differential quadrature examination of steady mixed convection nanofluid flows over an isothermal thin needle conveying metallic and metallic oxide nanomaterials: a comparative investigationArabian Journal for Science and Engineering2020In press10.1007/s13369-020-04420-xZaibA.KhanU.WakifA.ZaydanM.Numerical entropic analysis of mixed MHD convective flows from a non-isothermal vertical flat plate for radiative tangent hyperbolic blood biofuids conveying magnetite ferroparticles: dual similarity solutionsArabian Journal for Science and Engineering2020In press10.1007/s13369-020-04393-xWakifA.QasimM.AfridiM. I.SaleemS.Al-QarniM. M.Numerical examination of the entropic energy harvesting in a magnetohydrodynamic dissipative flow of Stokes’ second problem: utilization of the gear-generalized diferential quadrature methodJournal of Non-Equilibrium Thermodynamics201944438540310.1515/jnet-2018-00992-s2.0-85069756358ThummaT.Wakif A.AnimasaunI. L.Generalized differential quadrature analysis of unsteady three‐dimensional MHD radiating dissipative Casson fluid conveying tiny particlesHeat Transfer20204952595262610.1002/htj.21736AtanganaA.BaleanuD.New fractional derivatives with nonlocal and non-singular kernel: theory and application to heat transfer modelThermal Science201620276376910.2298/tsci160111018a2-s2.0-84971631213LiaoS.On the homotopy analysis method for nonlinear problemsApplied Mathematics and Computation2004147249951310.1016/s0096-3003(02)00790-72-s2.0-0141961626OwolabiK. M.AtanganaA.Robustness of fractional difference schemes via the Caputo subdiffusion-reaction equations, chaosSolitons & Fractals201811111912710.1016/j.chaos.2018.04.0192-s2.0-85045613657KilbasA. A.SrivastavaH. M.TrujilloJ. J.Theory and Applications of Fractional Differential Equations2006Amsterdam, NetherlandsNorth-Holland Mathematics Studies, Elsevier204SeneN.Second-grade fluid model with Caputo-Liouville generalized fractional derivativeChaos, Solitons & Fractals202013310963110.1016/j.chaos.2020.109631AbroK. A.Gomez-AguilarJ. F.A comparison of heat and mass transfer on a Walter’s-B fluid via Caputo-Fabrizio versus Atangana-Baleanu fractional derivatives using the Fox-H functionThe European Physical Journal Plus2019134311010.1140/epjp/i2019-12507-42-s2.0-85063084677AtanganaA.MekkaouiT.Trinition the complex number with two imaginary parts: fractal, chaos and fractional calculusChaos, Solitons & Fractals201912836638110.1016/j.chaos.2019.08.0182-s2.0-85070856107SeneN.Stokes' first problem for heated flat plate with Atangana-Baleanu fractional derivativeChaos, Solitons & Fractals2018117687510.1016/j.chaos.2018.10.0142-s2.0-85054825744SeneN.SIR epidemic model with Mittag-Leffler fractional derivativeChaos, Solitons & Fractals202013710983310.1016/j.chaos.2020.109833Al-MdallalQ.AbroK. A.KhanI.Analytical solutions of fractional walter’s B fluid with applicationsComplexity2018201810813132910.1155/2018/81313292-s2.0-85044167226AbdeljawadT.HajjiM. A.Al-MdallalQ.JaradF.Analysis of some generalized ABC-Fractional logistic modelsAlexandria Engineering Journal2020In press10.1016/j.aej.2020.01.030CaputoM.FabrizioM.A new definition of fractional derivative without singular kernelProgress in Fractional Differentiation and Applications201512115AbdeljawadT.AminR.ShahK.Al-MdallalQ.JaradF.Efficient sustainable algorithm for numerical solutions of systems of fractional order differential equations by Haar wavelet collocation methodAlexandria Engineering Journal2020In press10.1016/j.aej.2020.02.035AtanganaA.OwolabiK. M.New numerical approach for fractional differential equationsMathematical Modelling of Natural Phenomena2018131310.1051/mmnp/20180102-s2.0-85045204146ÖzdemirN.YavuzM.Numerical solution of fractional black-scholes equation by using the multivariate padé approximationActa Physica Polonica A20171321050105310.12693/aphyspola.132.10502-s2.0-85033391178