On the solution of two-dimensional fractional Black–Scholes equation for European put option

The purpose of this paper was to investigate the dynamics of the option pricing in the market through the two-dimensional time fractional-order Black–Scholes equation for a European put option. The Liouville–Caputo derivative was used to improve the ordinary Black–Scholes equation. The analytic solution is a powerful tool for describing the behavior of the option price in the European style market. In this study, analytic solution is carried out by the Laplace homotopy perturbation method. Moreover, the obtained solution showed that the Laplace homotopy perturbation method was an efficient method for finding an analytic solution of two-dimensional fractional-order differential equation.

Actually, one cannot find exact solutions to most fractional-order differential equations, so the approximate solutions are investigated to solve linear and nonlinear fractionalorder differential equations. There are many researches dealing with fractional order differential equations in many different fields [20][21][22][23].
The main feature of fractional order differential equations is memory. The variables in financial problems have long memories. Therefore, the fractional order differential equations are fitted to describe the financial problems. Many researches have been done to investigate the fractional-order differential equations in financial problems [24][25][26].
The Black-Scholes equation was proposed by Fisher Black and Myron Scholes in 1973 [27] to express the behavior of the option price in the European style market. Several papers have investigated how the Black-Scholes equation describes the behavior of the market [28][29][30][31][32]. Moreover, the Black-Scholes equation has been extended to express the behavior in another market, such as the American-style market and Asian-style market [33][34][35][36].
In general, the two-dimensional Black-Scholes equations for a European put option can be written as follows: with the terminal condition where E = max{E 1 , E 2 }, and the boundary conditions: where P denotes the value of a put option of the underlying stock prices {x 1 , x 2 } at time τ , q i denotes the dividend yield on the ith underlying stock, ρ ij denotes the correlation between the ith and jth underlying stock prices, T denotes the expiration date, r denotes the risk-free interest rate, σ i denotes the volatility of the ith underlying stock, E i denotes the strike price of the ith underlying stock, andβ i denotes a coefficient so that prices of all the risky assets are at the same level. There are many researchers who investigated analytical and approximate solutions of the Black-Scholes equation. Various effective methods have been used to solve the Black-Scholes equation, for example, the finite difference method [37], finite element method [15,38], homotopy perturbation method [39,40], the Mellin transform method [41], Adomian decomposition method [42], the variational iteration method [43,44], radial basis function partition of unity method (RBF-PUM) [45,46], and adaptive moving mesh method [47].
In this paper, we apply the Laplace homotopy perturbation method (LHPM) to get the analytical solution of the fractional-order Black-Scholes equation. This method is a combination of the Laplace transform and homotopy perturbation method [48]. It provides an analytical solution in the form of a convergent series. From the mathematical point of view, the analytical solution is a useful tool to describe the behavior of the solution, particularly the financial behavior of the solution to fractional-order differential equations.
The paper is structured as follows: the basic knowledge about fractional calculus is introduced in Sect. 2, and the two-dimensional time-fractional Black-Scholes equation is ( 2 )

Definition 3
The Liouville-Caputo-type fractional derivative of order 0 < α ≤ 1 for a function f : (0, ∞) → R is defined as [19]: In this paper, we use Liouville-Caputo derivative as the time-fractional derivative because the initial condition for the fractional order derivation is similar to the traditional derivative [49].

Generalized time fractional-order Black-Scholes model
The standard two-dimensional Black-Scholes equation for a European put option was presented in Sect. 1. Throughout this paper, we assume that σ 1 , σ 2 , ρ, and r are constants. By transforming coordinates and time to a forward time coordinate t = Tτ and applying the Liouville-Caputo fractional derivative to the two-dimensional Black-Scholes equation, we obtain the following two-dimension time-fractional order Black-Scholes equation with α ∈ (0, 1] and the initial and boundary conditions as follows: subject to the initial condition: and the boundary conditions: ω(x, y, t) = 0 as (x, y) → -∞, and where β 1 =β 1 e (r-1 2 σ 2 1 )T and β 2 =β 2 e (r-1 2 σ 2 2 )T .

Solving two-dimensional time-fractional Black-Scholes equation by LHPM
In this section, we apply LHPM techniques for finding the solution of two-dimensional time-fractional Black-Scholes model (5) subject to the initial condition (6) and the boundary conditions (7).

Conclusions
In this paper, we have found the analytical solution of the two-dimensional time-fractional Black-Scholes equation for a European put option using the Liouville-Caputo derivative. The Laplace homotopy perturbation method has been used to obtain the solution which can be written in the form of the generalized Mittag-Leffer function. The main feature of the analytical solution is that it is convenient for finding the option price.