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.

An option is the right to buy or sell a good at a predetermined price in the future. For customers or financial companies, knowing an option’s pricing is crucial. It is well recognized that the Black–Scholes model is an effective tool for estimating the cost of an option. The Black–Scholes equation has an explicit analytical solution known as the Black–Scholes formula. In some cases, such as the fractional-order Black–Scholes equation, there is no closed form expression for the modified Black–Scholes equation. This article shows how to find the approximate analytic solutions for the two-dimensional fractional-order Black–Scholes equation based on the generalized Riemann–Liouville fractional derivative. The generalized Laplace variational iteration method, which incorporates the generalized Laplace transform with the variational iteration method, is the methodology used to discover the approximate analytic solutions to such an equation. The expression of the two-parameter Mittag–Leffler function represents the problem’s approximate analytical solution. Numerical investigations demonstrate that the proposed scheme is accurate and extremely effective for the two-dimensional fractional-order Black–Scholes Equation in the perspective of the generalized Riemann–Liouville fractional derivative. This guarantees that the generalized Laplace variational iteration method is one of the effective approaches for discovering approximate analytic solutions to fractional-order differential equations.