^{1}

^{*}

^{1}

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.

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 [

Numerous methods and approaches have been presented in recent years. Some of these methods are analytical such as Fourier transform method [

Numerical methods and approaches are also popular and used to obtain approximate solutions of FPDEs. Examples of such a numerical method for solving FPDEs are the Homotopy Perturbation Method (HPM) [

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 [

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.

In this section, we present the basic definitions and properties of the fractional calculus theory, which are used further in this paper.

Definition 2.1

A real function f ( t ) , t > 0 is said to be in the space C σ , σ ∈ ℝ , if there exists a real number p > σ such that f ( t ) = t p f 1 ( t ) where f 1 ( t ) ∈ C [ 0 , ∞ ] and it is said to be in the space C σ m if f m ∈ C σ , m ∈ ℕ .

Definition 2.2

The left sided Riemann-Louiville fractional integral of order α ≥ 0 , of a function f ∈ C σ , σ ≥ − 1 is defined as:

J t α f ( t ) = 1 Γ ( α ) ∫ 0 t ( t − α ) α − 1 f ( ζ ) d ζ (1)

where α > 0 , t > 0 , and Γ ( m ) is the Gamma function. Properties of the operator J t α for f ∈ C μ , μ ≥ − 1 , α , β ≥ 0 , γ ≥ − 1 , can be founded in [

1) J t 0 f ( t ) = f ( t )

2) J t α f ( t ) J t β f ( t ) = J t α + β f ( t )

3) J t α f ( t ) J t β f ( t ) = J t β f ( t ) J t α f ( t )

4) J t β t γ = Γ ( γ + 1 ) Γ ( α + γ + 1 ) t α + γ + 1

Definition 2.3

Let f ∈ C μ m , n ∈ ℕ ∪ { 0 } . The left-Sided Caputo fractional of f in the Caputo sense is defined by [

D t α f ( t ) = f ( x ) = { 1 Γ ( n − α ) ∫ 0 t ( t − α ) n − 1 − α f n ( ζ ) d ζ , n − 1 < α ≤ n D t n f ( t ) , α = n (2)

According to the [

J t α f ( t ) = 1 Γ ( α ) ∫ 0 t ( t − ζ ) α − 1 f ( ζ ) d ζ (3)

D t α f ( t ) = 1 Γ ( n − α ) ∫ 0 t ( t − α ) n − 1 − α f n ( ζ ) d ζ (4)

Definition 2.4

The single parameter and the two parameters variants of the Mitting-Leffler function are denoted by E α ( t ) and E α , β ( t ) , respectively, which are relevant for their connection with fractional calculus, and are defined as:

E α ( t ) = ∑ j = 0 ∞ t j Γ ( α j + 1 ) , α > 0 , t ∈ ℂ (5)

E α , β ( t ) = ∑ j = 0 ∞ t j Γ ( α j + β ) , α , β > 0 , t ∈ ℂ (6)

the k-th derivatives are

E α ( t ) = d k d t k E α , β ( t ) = ∑ j = 0 ∞ ( k + j ) ! t j j ! Γ ( α j + α k + 1 ) , k = 0 , 1 , 2 , ⋯ , (7)

E α , β ( t ) = d k d t k E α , β ( t ) = ∑ j = 0 ∞ ( k + j ) ! t j Γ ( α j + α k + β ) , k = 0 , 1 , 2 , ⋯ , (8)

some special cases of the Mitting-Leffler function are as follows

1) E 1 ( t ) = e t

2) E α , 1 ( t ) = E α ( t )

3) d k d t k [ t β − 1 E α , β ( a t α ) ] = t β − α − 1 E α , β − k ( a t α )

Other properties of the Mitting-Leffler functions can be found in [

Definition 2.5

Sumudu transform over the following set of functions

A = { f ( t ) | ∃ M > 0 , τ 1 , τ 2 > 0 , | f ( t ) | > M e | t | τ j if t ∈ ( − 1 ) j × [ 0 , ∞ ] } (9)

is defined by

S [ f ( t ) ] = G ( u ) ∫ 0 ∞ f ( u t ) e − t d t (10)

where u ∈ ( τ 1 , τ 2 ) . Some special properties of the Sumudu transform are as follows

1) S [ 1 ] = 1

2) S [ t n Γ ( n + 1 ) ] = u n , n > 0

3) S [ e a t ] = 1 1 − a u

4) S [ α f ( x ) + β g ( x ) ] = α S [ f ( x ) ] + β S [ g ( x ) ]

Other properties of the Sumudu transform can be found in [

Definition 2.6

G ( u ) is the Sumudu transform of f ( t ) . And according to [

1) G ( 1 / s ) s is a meromorphic function, with singularities having R e ( s ) < γ .

2) There exists a circular region Γ with radius R and positive constants, M and k, with

| G ( 1 / s ) s | < M R − k

then the function f ( t ) is given by

S − 1 [ G ( s ) ] = 1 2 π i ∫ γ − i ∞ γ + i ∞ e s t G ( 1 s ) d s s = ∑ residuse [ e s t G ( 1 s ) s ] (11)

Definition 2.6

The Sumudu transform S , of the Caputo functional integral is defined as [

S [ D t α f ( t ) ] = G ( u ) u α − ∑ k = 0 n − 1 f ( k ) ( 0 ) u α − k (12)

then it can be easily understood that

S [ D t α f ( t ) ] = S [ f ( x , t ) ] u α − ∑ k = 0 n − 1 f ( k ) ( 0 ) u α − k , n − 1 ≤ α < n (13)

Theorem 2.1

Consider the following n-term liner fractional differential equation [

( a n D t β n + a n − 1 D t β n − 1 + ⋯ + a 0 D t β 0 ) ω t = f ( t ) (14)

with the constant initial condition

ω j i = C i j i , i = 0 , 1 , ⋯ , n , j = 1 , 2 , ⋯ , l i

where

a i , C i j i ∈ ℝ , n i − 1 < β i ≤ n i , n i ∈ ℕ and β 0 < β 1 < β 2 < ⋯ < β n − 1 < n ≤ β n < n + 1 .

Then we see that the analytical general solution of Equation (14) is

ω ( t ) = ∫ 0 t G n ( t − ζ ) f ( ζ ) d ζ + ∑ i = 0 ∞ ∑ j i = 0 l i − 1 a i C i j i G n B i − j i − 1 ( t ) (15)

where G n is the Green function and it’s defined by

G n ( t ) = 1 a n ∑ m = 0 ∞ ( − 1 ) m m ! × ∑ k 0 , k 1 , ⋯ , k n − 2 ≤ 0 , k 0 + k 1 + ⋯ + k n − 2 = m ( m ; k 0 , k 1 , ⋯ , k n − 2 ) × ∏ p = 0 n − 2 ( a p a n ) k p t ( β n − β n − 1 ) m + β n + ∑ j = 0 n − 2 ( β n − 1 − β j ) k j − 1 × E β n − β n − 1 , β n + ∑ j = 0 n − 2 ( β n − 1 − β j ) k j − 1 m ( − a n − 1 a n D β n − β n − 1 ) (16)

where ( m ; k 0 , k 1 , ⋯ , k n − 2 ) = m ! k 0 ! k 01 ! ⋯ k n − 2 ! .

And E ( . ) , ( . ) m is the m-th derivative of the Mitting-Leffler function.

In a special case of the latter theorem, the following relaxation-oscillation [

D t α ω ( t ) + A ω ( t ) = f ( t ) , t > 0 , ω i ( 0 ) = b i , i = 1 , 2 , ⋯ , n − 1 (17)

where b i are real constants, and n − 1 < α ≤ n . By utilizing theorem 2.1, we obtain the solution of Equation (14) as follows

ω ( t ) = ∫ 0 t G 2 ( t − ζ ) f ( ζ ) + ∑ j = 0 n − 1 b j D t α − j − 1 G 2 ( t ) (18)

where G 2 ( t ) = t α E α , α ( − A t α ) . It is easy to see that if 0 < α ≤ 1 , then the solution of the Equation (14) is given as follows

ω ( t ) = ∫ 0 t G 2 ( t − ζ ) f ( ζ ) d ζ + b 0 D t α − 1 G ( t ) (19)

which will be used in the coming examples, discussed in this work.

In this section, we will briefly present the main idea of the HPM. At first, we will consider the following nonlinear differential equation

A ( u ) − f ( x ) = 0 , x ∈ Ω (20)

where 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

L ( u ) + N ( u ) − f ( x ) = 0 , x ∈ Ω (21)

As a result, we construct a homotopy v ( x , p ) : Ω × [ 0 , 1 ] → ℝ that satisfies:

H ( v , p ) = ( 1 − p ) [ L ( u ) + N ( u 0 ) ] + p [ A ( u ) − f ( x ) ] = 0 , 0 ≤ p < 1 (22)

which equivalent to

H ( v , p ) = L ( v ) − N ( u 0 ) + p L ( u 0 ) + p [ N ( v ) + f ( x ) ] = 0 , 0 ≤ p < 1 (23)

when the value of p is varied from zero to one, we can see easily that

H ( v , p ) = L ( v ) − N ( u 0 ) = 0 (24)

H ( v , p ) = L ( v ) − N ( u ) − f ( x ) = A ( u ) − f ( x ) = 0 (25)

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

v = v 0 + p v 1 + p 2 v 2 + p 3 v 3 + ⋯ (26)

The best approximation for the solution of Equation (23) is

ω = lim p → 1 v = ∑ i = 0 ∞ v = v 0 + v 1 + v 2 + v 3 + ⋯ (27)

Meanwhile, we are able to apply the HPM to solve the class of time-fractional partial differential equations defined by

D t α ( ω ( x , t ) ) = L ( ω ( x , t ) ) + N ( ω ( x , t ) ) + f ( ω ( x , t ) ) (28)

subject to the initial condition

ω ( x , 0 ) = h ( x ) (29)

where 0 < α ≤ 1 . The constructed homotopy satisfies

H ( ω ( x , t ) ) = ( 1 − p ) [ D t α ( ω ( x , t ) ) − D t α ( ω 0 ( x , t ) ) ] + p [ D t α ( ω ( x , t ) ) − L ( ω ( x , t ) ) − N ( ω ( x , t ) ) − f ( ω ( x , t ) ) ] (30)

By simplifying Equation (30) we get

D t α ( ω ( x , t ) ) = D t α ( ω 0 ( x , t ) ) + p [ D t α ( ω ( x , t ) ) − L ( ω ( x , t ) ) − N ( ω ( x , t ) ) − f ( ω ( x , t ) ) ] (31)

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

ω ( x , t ) = ∑ k = 0 ∞ p k ω k ( x , t ) (32)

For the sake of clarifying the solution procedure of this method, we consider a general nonlinear time-fractional partial differential equation:

D t α ω ( x , t ) = L ω ( x , t ) + N ω ( x , t ) + f ( x , t ) (33)

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

p 0 : D t α ( ω 0 ( x , t ) ) = D t α ( h ( x ) ) p 1 : D t α ( ω 1 ( x , t ) ) = D t α ( ω 0 ( x , t ) ) − L ω 0 ( x , t ) − N ω 0 ( x , t ) − f ( x , t ) p 2 : D t α ( ω 2 ( x , t ) ) = D t α ( ω 1 ( x , t ) ) − L ω 1 ( x , t ) − N ω 1 ( x , t ) − f ( x , t ) ⋮ (34)

Applying the Riemann-Louiville fractional integral of order α ≥ 0 , J t α in the both sides

ω 0 ( x , t ) = h ( x ) ω 1 ( x , t ) = J t α [ D t α ( ω 0 ( x , t ) ) − L ω 0 ( x , t ) ( x , t ) − N ω 0 ( x , t ) − f ( x , t ) ] ω 2 ( x , t ) = J t α [ D t α ( ω 1 ( x , t ) ) − L ω 1 ( x , t ) ( x , t ) − N ω 1 ( x , t ) − f ( x , t ) ] ⋮ (35)

Substitute the results of the Equation (35) into the Equation (32), and applying the Perturbation technique [

ω ( x , t ) = ∑ i = 0 n ω i ( x , t ) (36)

If there exist some terms v n = 0 , n > 0 , then the exact solution can be written in the following form

ω ( x , t ) = ∑ i = 0 n − 1 ω i ( x , t ) (37)

We Assume that the initial approximation of Equation (28) is given by

ω 0 ( x , t ) = ω ( x , 0 ) c 1 ( t ) + ω ( x ) c 2 ( t ) (38)

where ω ( x , 0 ) is the initial condition of the Equation (28), and ω ( x ) = ∂ ω ∂ t ( x , 0 ) .

The aim of this algorithm is to find the terms c 1 and c 2 . We assume that ω n ≡ 0 , ∀ n ≥ 0 , which means that the exact solution is given by

ω ( x , t ) = ω 0 ( x , t ) (39)

As ω ( x , t ) satisfies also the initial condition, we get

ω ( x , 0 ) = ω 0 ( x , 0 ) = ω ( x , 0 ) c 1 ( 0 ) + ω ( x ) c 2 ( 0 ) = h ( x ) (40)

and c 1 ( 0 ) = 1 , c 2 ( 0 ) = 0 , on the other part we have,

D t α ( ω 1 ( x , t ) ) = D t α ( ω 0 ( x , t ) ) − L ω 0 ( x , t ) − N ω 0 ( x , t ) − f ( x , t ) ≡ 0 (41)

By substitution the Equation (38) into the Equation (28), we obtain

ω ( x , 0 ) D t α ( c 1 ( t ) ) + ω ( x ) D t α ( c 2 ( t ) ) = L ( ω ( x , 0 ) c 1 ( t ) + ω ( x ) c 2 ( t ) ) + N ( ω ( x , 0 ) c 1 ( t ) + ω ( x ) c 2 ( t ) ) + f ( x , t ) (42)

Therein, the FPDEs are changed to a fractional order differential equation, which simplifies the problem. When α is an integer, the equation is transformed to an ordinary differential equation.

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:

D t α ω ( x , t ) = L ω ( x , t ) + N ω ( x , t ) + g ( x , t ) (43)

where α > 0 and n − 1 ≤ α < n .

With the initial condition

D m ( ω ) ( 0 ) = { f m ( x ) , ( m = 0 , 1 , ⋯ , n − 1 ) 0 , m = n ≡ [ α ] (44)

where ω : = ω ( x , t ) .

Here by applying Sumudu Transform on the both sides of the Equation (43), we get

S [ D t α ω ( x , t ) ] = S [ L ω ( x , t ) ] + S [ N ω ( x , t ) ] + S [ g ( x , t ) ] (45)

which, upon using a property of the Sumudu transform, yields

S [ ω ] = f ( x ) + u α S [ L ω ] + u α S [ N ω ] + u α S [ g ( x , t ) ] (46)

where f is a function of x.

Taking the inverse of Sumudu on the both sides, we obtain

ω ( x , t ) = S − 1 [ ∑ k = 0 n − 1 u k ω ( k ) f ( x , 0 ) ] + p S − 1 [ u α S [ L ω ( x , t ) + N ω ( x , t ) + g ( x , t ) ] ] (47)

Now, we apply the HPM to the Equation (47) to get the solution

ω ( x , t ) = ∑ n = 0 ∞ p n ω n ( x , t ) (48)

where p ∈ [ 0 , 1 ] is the embedding (or homotopy) parameter. ω 0 is an initial approximation which satisfies the boundary conditions, and ω n , n ∈ ℕ are the nth order approximations which are functions yet to be determined. It is noted that setting p = 1 in the Equation (48) gives an approximate solution of the given nonlinear differential equation.

The nonlinear term can be decomposed as follows

N ( ω ( x , t ) ) = ∑ n = 0 ∞ p n H n ( ω 0 , ω 1 , ⋯ , ω n ) (49)

where H n are the He’s polynomials given by,

H n ( ω 0 , ω 1 , ⋯ , ω n ) = 1 n ! ∂ n ∂ p n N ( ∑ k = 0 n p k ω k ( x , t ) ) (50)

By substituting Equation (48) and Equation (49) into Equation (47), we obtain the following equation

∑ n = 0 ∞ p n ω n ( x , t ) = S − 1 [ ∑ k = 0 n − 1 u k f k ( x ) ] + p S − 1 [ u α S [ L ( ∑ n = 0 ∞ p n ω n ( x , t ) ) + ∑ n = 0 ∞ p n H n + g ( x , t ) ] ] (51)

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 ω n , n ∈ ℕ .

p 0 : ω 0 ( x , t ) = S − 1 [ ∑ k = 0 n − 1 u k f k ( x ) ] p 1 : ω 1 ( x , t ) = S − 1 [ u α S [ L ( ω 0 ) + H 0 ( ω 0 ) ] ] p 2 : ω 2 ( x , t ) = S − 1 [ u α S [ L ( ω 1 ) + H 1 ( ω 0 , ω 1 ) ] ] ⋮ p n : ω n ( x , t ) = S − 1 [ u α S [ L ( ω n − 1 ) + H 1 ( ω 0 , ω 1 , ω 2 , ⋯ , ω n − 1 ) ] ] (52)

Finally, the approximate and analytical solution of the Equation (43) is given by truncating the following series

ω ( x , t ) = lim p → 1 ∑ n = 0 ∞ p n ω n ( x , t ) = ∑ n = 0 ∞ ω n ( x , t ) (53)

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:

∂ α ψ ∂ x α = ∂ 2 ψ ∂ x 2 + ( k − 1 ) ∂ ψ ∂ x − k ψ (51)

ψ ( x , 0 ) = max ( e x − 1 , 0 ) (55)

By applying the MHPM, the initial approximation will be given by

ψ 0 ( x , t ) = ψ ( x , 0 ) c 1 ( t ) + ψ ( x ) c 2 ( t ) (56)

where

∂ ψ ∂ x ( x , 0 ) = max ( e x , 0 ) (57)

therefore

ψ 0 ( x , t ) = max ( e x − 1 , 0 ) c 1 ( t ) + max ( e x , 0 ) c 2 ( t ) (58)

Hence,

D t α ψ 1 = D t α ( max ( e x − 1 , 0 ) c 1 ( t ) + max ( e x , 0 ) c 2 ( t ) ) − ∂ 2 ∂ x 2 ( max ( e x − 1 , 0 ) c 1 ( t ) + max ( e x , 0 ) c 2 ( t ) ) + ( k − 1 ) ∂ ∂ x ( max ( e x − 1 , 0 ) c 1 ( t ) + max ( e x , 0 ) c 2 ( t ) ) − k ( max ( e x − 1 , 0 ) c 1 ( t ) + max ( e x , 0 ) c 2 ( t ) ) ≡ 0. (59)

We obtain the fractional system

{ D t α c 1 ( t ) + k c 1 ( t ) = 0 c 1 ( 0 ) = 0 (60)

{ D t α c 1 ( t ) + D t α c 2 ( t ) = 0 c 2 ( 0 ) = 0 (61)

Solving the Equation (60) and Equation (61) by applying the theorem 2.1, we obtain

{ c 1 ( t ) = E α ( − k t α ) c 2 ( t ) = 1 − E α ( − k t α ) (62)

and the exact solution is

ψ ( x , t ) = max ( e x − 1 , 0 ) E α ( − k t α ) + max ( e x , 0 ) ( 1 − E α ( − k t α ) ) (63)

If we put α → 1 in Equation (63) or solving Equation (60) and Equation (61) with α = 1 , we obtain the exact solution

ψ ( x , t ) = max ( e x − 1 , 0 ) e ( − k t α ) + max ( e x , 0 ) ( 1 − e ( − k t α ) ) (64)

Applying the Sumudu transform to the both sides of the Equation (54), we get

S [ ψ ( x , t ) ] = max ( e x − 1 , 0 ) + u α S [ ∂ 2 ψ ∂ x 2 + ( k − 1 ) ∂ ψ ∂ x − k ψ ] (65)

Operating the inverse Sumudu transform on both sides in the Equation (65), we have

ψ ( x , t ) = max ( e x − 1 , 0 ) + p S − 1 [ u α S [ ∂ 2 ψ ∂ x 2 + ( k − 1 ) ∂ ψ ∂ x − k ψ ] ] (66)

Now, applying the HPM

∑ n = 0 ∞ p n ψ n ( x , t ) = max ( e x − 1 , 0 ) + p ( S − 1 [ u α S [ ∑ n = 0 ∞ p n Ψ n ( x , t ) ] ] ) (67)

where

Ψ n ( x , t ) = ∂ 2 ψ n ∂ x 2 + ( k − 1 ) ∂ ψ n ∂ x − k ψ n (68)

Equating the corresponding power where, of p on both sides in Equation (67), we obtain

p 0 : ψ 0 ( x , t ) = max ( e x − 1 , 0 )

p 1 : ψ 1 ( x , t ) = S − 1 [ u α S [ Ψ 0 ( x , t ) ] ] = S − 1 [ u α S [ k max ( e x , 0 ) − k max ( e x − 1 , 0 ) ] ] = k t α Γ ( α + 1 ) max ( e x , 0 ) − k t α Γ ( α + 1 ) max ( e x − 1 , 0 )

p 2 : ψ 2 ( x , t ) = S − 1 [ u α S [ Ψ 1 ( x , t ) ] ] = ( k t ) 2 α Γ ( 2 α + 1 ) max ( e x , 0 ) − k ( − k t ) 2 α Γ ( 2 α + 1 ) max ( e x − 1 , 0 )

⋮

p n : ψ n ( x , t ) = S − 1 [ u α S [ Ψ n − 1 ( x , t ) ] ] = ( k t ) n α Γ ( n α + 1 ) max ( e x , 0 ) − k ( − k t ) n α Γ ( n α + 1 ) max ( e x − 1 , 0 ) (69)

So that the solution ψ ( x , t ) of the problem is given by

ψ ( x , t ) = lim p → 1 ∑ n = 0 ∞ p n ψ n ( x , t ) = max ( e x − 1 , 0 ) E α ( − k t α ) + max ( e x , 0 ) ( 1 − E α ( − k t α ) ) (70)

where E α ( x ) is the Mittag-Leffler function in one parameter. For special case, α = 1 , we get

ψ ( x , t ) = max ( e x − 1 , 0 ) e ( − k t α ) + max ( e x , 0 ) ( 1 − e ( − k t α ) ) (71)

∂ α ψ ∂ x α + 0.08 ( 2 + sin x ) 2 x ∂ 2 ψ ∂ x 2 + 0.06 x ∂ ψ ∂ x − 0.06 ψ = 0 , 0 < α ≤ 1 (72)

Subject to the initial condition,

ψ ( x , 0 ) = max ( x − 25 e x , 0 ) (73)

By applying the MHPM, the initial approximation will be given by

ψ 0 ( x , t ) = ψ ( x , 0 ) c 1 ( t ) + ψ ( x ) c 2 ( t ) (74)

where,

∂ ψ ∂ x ( x , 0 ) = 1 (75)

Therefore,

ψ 0 ( x , t ) = max ( x − 25 e x , 0 ) c 1 ( t ) + c 2 ( t ) (76)

Thus,

max ( x − 25 e x , 0 ) + D t α c 1 ( t ) + D t α c 2 ( t ) + 0.06 c 1 ( t ) − 0.06 max ( x − 25 e x , 0 ) c 1 ( t ) − 0.06 c 2 ( t ) = 0 (77)

We obtain the final system

{ D t α c 1 ( t ) − 0.06 c 1 ( t ) = 0 c 1 ( 0 ) = 0 (78)

{ D t α c 2 ( t ) + 0.06 x c 1 ( t ) − 0.06 c 2 ( t ) = 0 c 2 ( 0 ) = 0 (79)

Solving the Equation (49) and Equation (50), we get

{ c 1 ( t ) = E α ( 0.06 t α ) c 2 ( t ) = x ( 1 − E α ( 0.06 t α ) ) (80)

Hence,

ψ ( x , t ) = max ( x − 25 e − 0.06 , 0 ) E α ( 0.06 t α ) + x ( 1 − E α ( 0.06 t α ) ) (81)

This is the exact solution of the given option pricing Equation (72). The solution at the special case α = 1 , is given as follows

ψ ( x , t ) = max ( x − 25 e − 0.06 , 0 ) e 0.06 t + x ( 1 − e 0.06 t ) (82)

Applying the Sumudu transform on both sides of the Equation subject to initial condition given in Equation , we get

S [ ψ ( x , t ) ] = max ( x − 25 e − 0.06 , 0 ) − u α S [ 0.08 ( 2 + sin x ) 2 x ∂ 2 ψ ∂ x 2 + 0.06 x ∂ ψ ∂ x − 0.06 ψ ] (83)

Operating the inverse Sumudu transform on both sides in the Equation (83), we have

ψ ( x , t ) = max ( x − 25 e − 0.06 , 0 ) − p S − 1 { u α S [ 0.08 ( 2 + sin x ) 2 x ∂ 2 ψ ∂ x 2 + 0.06 x ∂ ψ ∂ x − 0.06 ψ ] } (84)

Now, applying the homotopy perturbation method we have

∑ n = 0 ∞ p n ψ n ( x , t ) = max ( x − 25 e − 0.06 , 0 ) − p S − 1 { u α S [ ∑ n = 0 ∞ p n Ψ n ( x , t ) ] } (85)

where,

Ψ n ( x , t ) = 0.08 ( 2 + sin x ) 2 x ∂ 2 ψ n ∂ x 2 + 0.06 x ∂ ψ n ∂ x − 0.06 ψ n (86)

Equating the corresponding power of p on both sides in Equation (85), we obtain

p 0 : ψ 0 ( x , t ) = max ( x − 25 e − 0.06 , 0 )

p 1 : ψ 1 ( x , t ) = S − 1 [ u α S [ Ψ 0 ( x , t ) ] ] = S − 1 [ u α S [ 0.06 x − 0.06 max ( x − 25 e − 0.06 , 0 ) ] ] = x − 0.06 t α Γ ( α + 1 ) + 0.06 t α Γ ( α + 1 ) max ( x − 25 e − 0.06 , 0 )

p 2 : ψ 2 ( x , t ) = S − 1 [ u α S [ Ψ 1 ( x , t ) ] ] = x ( − 0.06 t ) 2 α Γ ( 2 α + 1 ) + ( 0.06 t ) 2 α Γ ( 2 α + 1 ) ( 0.06 t ) 2 α Γ ( 2 α + 1 ) max ( x − 25 e − 0.06 , 0 )

⋮

p n : ψ n ( x , t ) = S − 1 [ u α S [ Ψ n − 1 ( x , t ) ] ] = x ( − 0.06 t ) n α Γ ( n α + 1 ) + ( 0.06 t ) n α Γ ( n α + 1 ) ( 0.06 t ) n α Γ ( n α + 1 ) max ( x − 25 e − 0.06 , 0 ) (87)

So that the solution ψ ( x , t ) of the problem is given by

ψ ( x , t ) = lim p → 1 ∑ n = 0 ∞ p n ψ n ( x , t ) = max ( x − 25 e − 0.06 , 0 ) E α ( 0.06 t α ) + x ( 1 − E α ( 0.06 t α ) ) (88)

This is the exact solution of the given option pricing equation Equation (72). The exact solution at the special case α = 1 is

ψ ( x , t ) = max ( x − 25 e − 0.06 , 0 ) e 0.06 t + x ( 1 − e 0.06 t ) (89)

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.

Ouafoudi, M. and Gao, F. (2018) Exact Solution of Fractional Black-Scholes European Option Pricing Equations. Applied Mathematics, 9, 86-100. https://doi.org/10.4236/am.2018.91006