Exact Solution of Fractional Black-Scholes European Option Pricing Equations

We introduce two algorithms in order to find the exact solution of the nonlinear Time-fractional Partial differential equation, in this research work. Those algorithms are proposed in the following structure: The Modified Homotopy Perturbation Method (MHPM), The Homotopy Perturbation and Sumudu Transform Method. The results achieved using the both methods are the same. However, we calculate the approached theoretical solution of the Black-Scholes model in the form of a convergent power series with a regularly calculated element. Finally, we propose a descriptive example to demonstrate the efficiency and the simplicity of the methods.


Introduction
In recent years, fractional calculus has been increasingly used for numerous applications in many scientific and technical fields such as medical sciences, biological research, as well as various chemical, biochemical and physical fields.
Fractional calculus can be, for instance, employed to solve a lot of problems within the biomedical research field.Such an important application is studying membrane biophysics and polymer viscoelasticity [1].
Numerous methods and approaches have been presented in recent years.Some of these methods are analytical such as Fourier transform method [2], the fractional Green function method [3], the popular Laplace method [4] [5], the iteration method [6], the Mellin transform method and the method of the orthogonal polynomial [4].
Among these numerical methods, the VIM and the ADM are the most popular ones that are used to solve differential and integral equations of integer and fractional order.The HPM is a universal approach which can be used to solve both fractional and ordinary differential equations FODs as well as fractional partial differential equations FPDEs.The HPM method was originally proposed by He [8] [17].The HPM is a coupling of homotopy in the topology and the perturbation method.The method is used to solve various types of equations such us the Hellmholtz equation, the fifth order KdV [6], the Kleein-Gorden Equation [4], the Fokker-Plank equation [9] [18], the nonlinear Kolmogorov-Petrovskii-Piskunov Equation [3] as well as other types of equations as proposed and used in [10] [19].
In this paper, we present two approaches to derive the exact solution of various types of FPDEs.Those methods are: The Modified Homotopy Perturbation Method (MHPM) and the Homotopy Perturbation and Sumudu Transform Method (HPSTM).The MHPM is a fast approach that's based on designing and utilizing a proper initial approximation which satisfies the initial condition of the HPM.However, the HPSTM is a combination of the homotopy perturbation method and Sumudu transform.The paper is structured in six sections.In Section 2, we begin with an introduction to some necessary definitions of fractional calculus theory.In Section 3, we describe the basic idea of the HPM.In Section 4, we describe the MHPM.In Section 5, we describe the HPSTM.In Section 6, we present two examples to show the efficiency of using the MHPM and HPSTM to solve the fractional Black-Scholes equation.Finally, relevant conclusions are drawn in Section 7.

Basic Definitions of Fractional Calculus
In this section, we present the basic definitions and properties of the fractional calculus theory, which are used further in this paper.
The left sided Riemann-Louiville fractional integral of order 0 α ≥ , of a func- tion , 1 where 0 α > , 0 t > , and ( ) is the Gamma function.Properties of the op- , can be founded in [20], and are defined as follow 1) ( ) ( ) The left-Sided Caputo fractional of f in the Ca- puto sense is defined by [4] as follows According to the [20], Equations ( 1) and (2) becomes, for Definition 2.4 The single parameter and the two parameters variants of the Mitting-Leffler function are denoted by ( ) ) , respectively, which are relevant for their connection with fractional calculus, and are defined as: , , 0, the k-th derivatives are some special cases of the Mitting-Leffler function are as follows 1) ( ) Other properties of the Mitting-Leffler functions can be found in [20].These functions, are generalizations of the exponential function, because, most linear differential equations of fractional order have solutions that are expressed in terms of these functions.
Definition 2.5 Sumudu transform over the following set of functions is defined by where , u τ τ ∈ . Some special properties of the Sumudu transform are as follows 1) [ ] Other properties of the Sumudu transform can be found in [21].
Definition 2.6 ( ) G u is the Sumudu transform of ( ) f t .And according to [22] we have: 2) There exists a circular region Γ with radius R and positive constants, M and k, with ( ) The Sumudu transform  , of the Caputo functional integral is defined as [22] ( ) ( ) ( ) ( ) then it can be easily understood that Theorem 2.1 Consider the following n-term liner fractional differential equation [7] theorem ( ) ( ) with the constant initial condition , Then we see that the analytical general solution of Equation ( 14) is where n G is the Green function and it's defined by ( ) ; , , , !
In a special case of the latter theorem, the following relaxation-oscillation [7] is solved: where i b are real constants, and 1 n n α − < ≤ .By utilizing theorem 2.1, we obtain the solution of Equation ( 14) as follows . It is easy to see that if 0 1 α < ≤ , then the solu- tion of the Equation ( 14) is given as follows which will be used in the coming examples, discussed in this work.

The Basic Idea of the HPM
In this section, we will briefly present the main idea of the HPM.At first, we will consider the following nonlinear differential equation , , A B f x and Gamma are a general differential function operator, a boundary operator, a known analytical function and the boundary of the domain Ω , respectively.The operator A can be decomposed into a linear operator denoted by L, and nonlinear operator denoted by N. Hence, the equation can be written as follows As a result, we construct a homotopy ( ) [ ] when the value of p is varied from zero to one, we can see easily that If the parameter p is assumed as small, then the solution of the Equation ( 23) can be expressed as a power series in p as follows The best approximation for the solution of Equation ( 23) is

The Modified Homotopy Perturbation Method (MHPM)
Meanwhile, we are able to apply the HPM to solve the class of time-fractional partial differential equations defined by ( subject to the initial condition where 0 1 α < ≤ .The constructed homotopy satisfies By simplifying Equation (30) we get where the embedding parameter p is considered to be small, and applied to the classical perturbation technique.The next step is to use the homotopy parameter p to expand the solution into the following form For the sake of clarifying the solution procedure of this method, we consider a general nonlinear time-fractional partial differential equation: By substituting Equation (32) into Equation (31) and equating the terms with identical power of p, we can obtain a series of equations as the following Applying the Riemann-Louiville fractional integral of order 0, t J α α ≥ in the both sides Substitute the results of the Equation (35) into the Equation (32), and applying the Perturbation technique [15] [22] we obtain an authentic n-th approximation of the exact solution, given by ( ) ( ) 0 , , If there exist some terms 0, 0 vn n = > , then the exact solution can be written in the following form We Assume that the initial approximation of Equation ( 28) is given by where ( ) ,0 x ω is the initial condition of the Equation (28), and ( ) ( ) The aim of this algorithm is to find the terms 1 c and 2 c .We assume that 0 n ω ≡ , 0 n ∀ ≥ , which means that the exact solution is given by As ( ) , x t ω satisfies also the initial condition, we get and ( ) ( ) By substitution the Equation (38) into the Equation (28), we obtain Therein, the FPDEs are changed to a fractional order differential equation, which simplifies the problem.When α is an integer, the equation is trans- formed to an ordinary differential equation.

The Perturbation Sumudu Transform Method (HPSTM)
We illustrate the basic idea of HPSTM by considering the general time-fractional nonlinear non-homogeneous partial differential equation with the initial condition of the general form defined by: where 0 α > and where ( ) Here by applying Sumudu Transform on the both sides of the Equation (43), we get which, upon using a property of the Sumudu transform, yields where f is a function of x.
Taking the inverse of Sumudu on the both sides, we obtain The nonlinear term can be decomposed as follows where n H are the He's polynomials given by, ( ) ( ) ( ) By substituting Equation (48) and Equation (49) into Equation (47), we obtain the following equation In the meantime, we have completed the coupling of the Sumudu transform and the HPM using He's polynomials.By equating the terms of the same power of p in the Equation (51), we obtain the given approximation for Finally, the approximate and analytical solution of the Equation ( 43) is given by truncating the following series ( ) ( ) ( )

Applications
In order to apply the MHPM and the HPSTM in Sections 4 and 5, respectively, consider the following fractional Black-Scholes option pricing equations as follows: , 0 max e 1, 0

Conclusion
It is widely known that various scientific models in the life sciences end up with a PDEs or FPDEs.Solving these equations is very important to understand many phenomena in life sciences.In this work, two methods are applied: MHPM and HPSTM, in order to find the solution of the Fractional Black-Scholes model.The both methods are based on the HPM.However, they are effectively employed for getting the solution.To conclude, MPHM and HPSTM are exceptionally effective and productive methods to discover the approached solution, also the numerical solution for the Time-Fractional Partial Differential Equations.

1 α 6 . 1 . 2 .
(60)  and Equation (61) by applying the theorem 2.→ in Equation (63) or solving Equation (60) and Equation (61) with 1 α = , we obtain the exact solution Use the HPSTM Applying the Sumudu transform to the both sides of the Equation (54), we get ) This is the exact solution of the given option pricing equation Equation (72).