Homotopy Perturbation Method for Fractional Black-Scholes European Option Pricing Equations Using Sumudu Transform

The homotopy perturbation method, Sumudu transform, and He’s polynomials are combined to obtain the solution of fractional Black-Scholes equation. The fractional derivative is considered in Caputo sense. Further, the same equation is solved by homotopy Laplace transform perturbation method. The results obtained by the two methods are in agreement. The approximate analytical solution of Black-Scholes is calculated in the form of a convergence power series with easily computable components. Some illustrative examples are presented to explain the efficiency and simplicity of the proposed method.


Introduction
Fractional differential equations have attracted much attention, recently, see, for instance [1][2][3][4][5].This is mostly due to the fact that fractional calculus provides an efficient and excellent instrument for the description of many practical dynamical phenomena arising in engineering and scientific disciplines such as, physics, chemistry, biology, economy, viscoelasticity, electrochemistry, electromagnetic, control, porous media, and many more, see, for example, [6][7][8][9].
Many partial differential equations of fractional order have been studied and solved.For example many researchers studied the existence of solutions of the Black-Scholes model using many methods, see [10][11][12][13][14].
The homotopy perturbation method was first introduced and applied by He [15][16][17].This method has been applied by many authors in many fields, for example, it is applied to nonlinear oscillator [18], nonlinear wave equation [19], nonlinear partial differential equations [20], integro-differential equation of fractional order [21], fuzzy differential equation [22], and other fields [23,24].Further homotopy perturbation methods are combined with Laplace transform to solve many problems such as one dimensional nonhomogeneous partial differential equations with a variable coefficient [25], Black-Scholes of fractional order [26], and parabolic partial differential equations [27].The homotopy perturbation method coupled with Sumudu transform basically illustrates how Sumudu transform can be used to approximate the solutions of the linear and nonlinear differential equations by manipulating the homotopy perturbation method.In [28] Singh et al. studied the solution of linear and nonlinear partial differential equations by using the homotopy perturbation method coupled with Sumudu transform.Further, in [29] the authors proposed the homotopy perturbation method coupled with Sumudu transform to solve nonlinear fractional gas dynamics equation.
The Black-Scholes equation is one of the most significant mathematical models for a financial market.It is a secondorder parabolic partial differential equation that governs the value of financial derivatives.This Black-Scholes model for the value of an option is described by the following equation: where V(, ) is the European call option price at asset price  and at time ,  is the maturity, () is the risk free interest rate, and (, ) represents the volatility function of underlying asset.The payoff functions are where V  (, ) and V  (, ) are the value of the European call and put options, respectively,  denotes the expiration price for the option, and the function max(, 0) gives the large value between  and 0.
In this paper, we consider the following fractional Black-Scholes of the form In [29] Singh et al. used homotopy perturbation method coupled with Sumudu transform to solve fractional gas dynamics equation.The aim of this paper is to applied the homotopy perturbation method for fractional Black-Scholes equation by using He's polynomials and Sumudu transform.

Sumudu Transform
The Sumudu transform was first introduced and applied by Watugala [30] in (1998).For further details and properties of Sumudu transform see [31][32][33][34].The Sumudu transform is defined over the set of functions: by the following formula The existence and uniqueness was discussed in [35].For further properties of Sumudu transform and its derivatives, see [36].Some fundamental further established properties of Sumudu transform can be found in [31].
Similarly, this new transform was applied to one-dimensional neutron transport equation [37].In [34] Kılıc ¸man et al.
show that there is a strong relationship between Sumudu and other integral transforms.Further in [33] the Sumudu transform was extended to the distributions, and some of their properties were also studied in [38].Recently Kılıc ¸man et al. applied this transform to solve system of differentials equations, for more details see [34,35,[37][38][39].

Basic Definitions of Fractional Calculus
In this section, we give some basic definitions and properties of fractional calculus theory which will be used in this paper.
For Riemann-Liouville fractional integral, one has Definition 2. The Caputo fractional derivative of  ∈   −1 ,  ∈  is defined as follows: Lemma 3. If  − 1 <  ≤ ,  ∈ ,  ∈    ,  > −1 then the following two properties hold: Definition 4. The Mittag-Leffler function   () with  > 0 is defined by the following series representation, valid in the whole complex plane: Definition 5.The Sumudu transform of the Caputo fractional derivative is defined as follows [40]:

Homotopy Perturbation Method
To illustrate the basic idea of this method, we consider the following nonlinear differential equation: with boundary conditions where  is a general differential operator,  is a boundary operator, () is a known analytic function, and Γ is the boundary of the domain Ω.
In general, the operator  can be divided into two parts  and , where  is linear, while  is nonlinear.Equation (11) therefor can be rewritten as follows: By the homotopy technique [41,42] we construct a homotopy V(, ) : Ω × [0, 1] →  which satisfies or where  ∈ [0, 1] is an embedding parameter, and  0 is an initial approximation of ( 11) which satisfies the boundary conditions.From ( 14) and ( 15) we have The changing in the process of  from zero to unity is just that of V(, ) from  0 () to ().In topology this is called deformation, and (V) − ( 0 ) and (V) − () are called homotopic.Now, assume that the solution of ( 14), ( 15) can be expressed as Setting  = 1 results in the approximate solution of (11):

Homotopy Perturbation Method Coupled with Sumudu Transform
To illustrate the basic idea of this method, we consider the following nonlinear fractional differential equation: where    =   /  is the fractional Caputo derivative of the function (, ),  is the linear differential operator,  is the nonlinear differential operator, and (, ) is the source term.Now, applying the Sumudu transform on both sides of (19), we have Operating with Sumudu inverse on both sides of ( 21) where (, ) represents the term arising from the source term and the prescribed initial conditions.Now, applying the classical homotopy perturbation technique, the solution can be expressed as a power series in  as given below: where the homotopy parameter  is considered as a small parameter ( ∈ [0, 1]).
We can decompose the nonlinear term as where   are He's polynomials of  0 ,  1 ,  2 , . . .,   [43][44][45], and it can be calculated by the following formula: By substituting ( 23) and ( 24) and using HPM [15] we get This is coupling of Sumudu transform and homotopy perturbation method using He's polynomials.By equating the coefficient of corresponding power of  on both sides, the following approximations are obtained as  0 :  0 (, ) =  (, ) ,

Mathematical Problems in Engineering
Proceeding in the same manner, the rest of the components   (, ) can be completely obtained, and the series solution is thus entirely determined.Finlay we approximate the solution (, ) by truncated series This series solutions generally converge very rapidly.

Examples
In this section, we discuss the implementation of the proposed method.
Example 6.We consider the following fractional Black-Scholes option pricing equation as follows: subject to initial condition Applying Sumudu transform on both sides of (29) subject to initial condition (30), we get Operating the inverse Sumudu transform on both sides in (31), we have Now, applying homotopy perturbation method where Equating the corresponding power of  on both sides in (38), we have So that the solution V(, ) of the problem is given by where   () is Mittag-Leffler function in one parameter.For special case  = 1, we get which is an exact solution of the given Black-Scholes equation (29) for  = 1.
The behaviour of V (, ) with respect to  and  when  = 1 is given in Figure 1.
Example 7. We consider the following fractional Black-Scholes option pricing equation as follows: subject to initial condition Firstly, applying Sumudu transform on both sides of (38) subject to initial condition (39), we get Operating the inverse Sumudu transform on both sides in (40), we have Now, applying the homotopy perturbation method we have where Equating the corresponding power of  on both sides in (42), we have So that the solution V(, ) of the problem is given by V (, ) = lim This is the exact solution of the given option pricing equation (38).The solution of (38) at the special case  = 1 is V (, ) =  (1 −  0.06 − 1, 0) + max ( − 25 −0.06 , 0)  −0.06 . ( The behaviour of V(, ) with respect to  and  when  = 1 is given in Figure 2.

Conclusion
In this paper, the homotopy perturbation Sumudu transform method (HPSTM) is successfully applied for getting the analytical solution of the fractional Black-Scholes option pricing equation.Two examples from the literature [26] are presented.The results of the illustrated examples are in agreement with the results of the method presented in [26].In conclusion, HPSTM is a very powerful and efficient method to find approximate solutions as well as numerical solutions.

Figure 1 :
Figure 1: The surface shows the V(, ) for (29) with respect to  and  for  = 1.

Figure 2 :
Figure 2: The surface shows the V(, ) for (38) with respect to  and  for  = 1.