It is well known that the Black-Scholes model is used to establish the behavior of the option pricing in the financial market. In this paper, we propose the modified version of Black-Scholes model with two assets based on the Liouville-Caputo fractional derivative. The analytical solution of the proposed model is investigated by the Laplace transform homotopy perturbation method.

In the finance market, the Black–Scholes equation is used to model the price change of the underlying fractal transmission system. Moreover, the fractional differential equations recently are accepted by researchers that fractional differential equations are a powerful tool in studying fractal geometry and fractal dynamics. Fractional differential equations are used in modeling the various important situations or phenomena in the real world such as fluid flow, acoustics, electromagnetic, electrochemistry and material science. There is an important question in finance: “Can the fractional differential equation be applied in the financial market?”. The answer is “Yes”. Due to the self-similar property of the fractional derivative, it can reply to the long-range dependence better than the integer-order derivative. Thus, these advantages are beneficial to manage the fractal structure in the financial market. In this article, the classical Black–Scholes equation with two assets for the European call option is modified by replacing the order of ordinary derivative with the fractional derivative order in the Caputo type Katugampola fractional derivative sense. The analytic solution of time-fractional Black–Scholes European call option pricing equation with two assets is derived by using the generalized Laplace homotopy perturbation method. The used method is the combination of the homotopy perturbation method and generalized Laplace transform. The analytic solution of the time-fractional Black–Scholes equation is carried out in the form of a Mittag–Leffler function. Finally, the effects of the fractional-order in the Caputo type Katugampola fractional derivative to change of a European call option price are shown.

In the finance market, it is well known that the price change of the underlying fractal transmission system can be modeled with the Black-Scholes equation. This article deals with finding the approximate analytic solutions for the time-fractional Black-Scholes equation with the fractional integral boundary condition for a European option pricing problem in the Katugampola fractional derivative sense. It is well known that the Katugampola fractional derivative generalizes both the Riemann–Liouville fractional derivative and the Hadamard fractional derivative. The technique used to find the approximate analytic solutions of the time-fractional Black-Scholes equation is the generalized Laplace homotopy perturbation method, the combination of the generalized Laplace transform and homotopy perturbation method. The approximate analytic solution for the problem is in the form of the generalized Mittag-Leffler function. This shows that the generalized Laplace homotopy perturbation method is one of the most effective methods to construct approximate analytic solutions of the fractional differential equations. Finally, the approximate analytic solutions of the Riemann–Liouville and Hadamard fractional Black-Scholes equation with the European option are also shown.

This article deals with the fractional multi-dimensional Burgers equation in the sense of the Caputo fractional derivative. An approximate analytical solution of the problem is established by the homotopy perturbation method (HPM). Furthermore, the convergence analysis and the error estimation derived by the HPM are shown.