Abstract

The residual power series method is effective for obtaining solutions to fractional-order differential equations. However, the procedure needs the derivative of the residual function. We are all aware of the difficulty of computing the fractional derivative of a function. In this article, we considered the simple and efficient method known as the Laplace residual power series method (LRPSM) to find the analytical approximate and exact solutions of the time-fractional Black–Scholes option pricing equations (BSOPE) in the sense of the Caputo derivative. This approach combines the Laplace transform and the residual power series method. The suggested method just needs the idea of an infinite limit, so the computations required to determine the coefficients are minimal. The obtained results are compared in the sense of absolute errors against those of other approaches, such as the homotopy perturbation method, the Aboodh transform decomposition method, and the projected differential transform method. The results obtained using the provided method show strong agreement with different series solution methods, demonstrating that the suggested method is a suitable alternative tool to the methods based on He’s or Adomian polynomials.

1. Introduction

Fractional calculus (FC) deals with fractional derivatives and integrations. The pioneers of FC were two mathematicians, Leibniz and L’Hospital, and the date September 30, 1695, is regarded as its exact birthday. Many scientists and researchers have been drawn to FC in recent decades because it is commonly used in scientific contexts such as engineering, image processing, physics, biochemistry, biology, fluid mechanics, and entropy theory [15]. Although fractional derivatives can be defined in a variety of ways, not all of them are generally used. The Atangana–Baleanu, Riemann–Liouville (R-L), Caputo–Fabrizio, Caputo, and conformable operators are the most frequently used [612]. In some cases, fractional derivatives are preferable to integer-order derivatives for modeling because they can simulate and analyze complicated systems having complicated nonlinear processes and higher-order behaviors. There are two main causes of this. First, we can select any order for the derivative operator, rather than being restricted to an integer order. Second, depending on both the past and current circumstances, noninteger order derivatives are advantageous when the system has long-term memory. The primary components of FC, which is the generalized version of classical calculus and has piqued the interest of numerous academics and scientists due to its wide range of applications, are fractional order differential equations (FODEs). The FODEs are frequently used for their logical support in the mathematical framework of physical problems, including technology, healthcare, monetary markets, and decision theory. As a result, the solutions provided by the FODEs are significant and useful. Applications regularly face FODEs that are too complex for close-form solutions. Under the specified initial and boundary conditions, numerical methods present a potent alternative tool for solving FODEs. Several numerical techniques, including the Shehu decomposition approach [13], the differential transform method [14], the variational iteration method [15], the operational matrix approach [16], the homotopy analysis technique [17], the Aboodh transform decomposition method [18], the finite difference method [19], the fractional power series method [20], the Chebyshev polynomials method [21], the residual power series method [22], and the natural transform homotopy perturbation method [23], have been developed in recent years for solving FODEs.

It is efficient to obtain approximate analytical solutions to FODEs using the residual power series method (RPSM). However, the algorithm requires the derivative of the residual function. We are all aware of how difficult it is to calculate a function’s fractional derivative. As a result, the use of conventional RPSM is somewhat constrained. To overcome the limitations of the RPSM, Eriqat et al. introduced a new technique called LRPSM [24]. Several FODEs have been solved through the recommended method [2429]. The provided equation is transformed into the Laplace transform (LT) space in accordance with the set of rules for this novel approach; the solution to the new form of the problem is established; the solution to the original equation is then achieved by applying the inverse LT.

The pricing of financial derivatives is a topic that has generated a great deal of interest and literature. A financial derivative is an asset whose price is based on the value of another asset. Frequently, a stock or bond serves as the underlying asset. Financial derivatives are not a completely novel concept. It is generally agreed that Charles Castelly’s work, which was published in 1877, was the first attempt at contemporary derivative pricing, despite some historical controversy regarding the exact year of the birth of financial derivatives [30]. Although it lacked mathematical rigor, Castell’s book served as a general introduction to ideas such as hedging and speculative trading.

The first widely used mathematical method to calculate the theoretical value of an option contract using prevailing equity markets’ predicted dividends, the option’s strike price, projected interest rates, time till the cessation, and expected unpredictability was the Black–Scholes (BS) model, developed in 1997 by Fischer Black, Robert Merton, and Myron Scholes. The Pricing of Options and Corporate Liabilities by Black and Scholes, published in the Journal of Political Economy in 1973, provided the initial formulation of the equation. Robert C. Merton contributed to the editing of that document. Later that year, he wrote his own work, “Theory of Rational Option Pricing,” expanding the model’s mathematical capabilities and applications while also coining the term “BS theory of options pricing” [31]. One of the most significant mathematical representations of a financial market is the BS equation. The value of financial derivatives is controlled by a second-order parabolic partial differential equation. Many different commodities and payout structures have been used in the BS model for pricing stock options. The following equation describes the BS model for the value of an option [32]:where is the order of the Caputo derivative (CD), is the European offer’s value at the fundamental market cap of and time , the fluctuation component of the company’s shares, often referred to as , quantifies the variance of the stock’s return, is the expiration date, and is the risk-free interest rate. and , correspondingly, stand for the call and put values on European options. Then, the payoff functions are given by and , where indicates the termination price for the option and the function gives the larger value between and 0. The important financial view of the BS equation is that it minimizes risk by allowing for the prudent selection of purchasing and selling the stock under scrutiny. It is a sign that, according to the BS financial model, there is only one correct price for the option. The important financial view of the BS equation is that it minimizes risk by allowing for the prudent selection of purchasing and selling the stock under scrutiny. It is a sign that, according to the BS financial model, there is only one correct price for the option. In this article, a fractional model that may be used to model the pricing of various financial derivatives is presented.

Numerous techniques have been used to examine the time-fractional BSOPE [3340]. Each of these techniques has specific restrictions and flaws. These techniques involve a lot of computing work and a long running time. In this article, we considered LRPSM, which is a simple and efficient technique to solve BSOPE. The advantage of the recommended method over the homotopy perturbation method (HPM) and the Adomian decomposition method (ADM) is its strength in handling problems without the use of He’s and Adomian polynomials. The advantage of this approach is that the problem does not involve any physical parameter assumptions, no matter how big or small. For a series such as the RPSM, the coefficients must be determined each time using the fractional derivative. Since LRPSM just needs the idea of an infinite limit, the computations required to determine the coefficients are minimal. The LRPSM results are also compared to those of other approaches, including the projected differential transform method (PDTM) [34], the Aboodh transform decomposition method (ATMD) [35], and the homotopy perturbation method (HPM) [36]. The results obtained using the suggested method show strong agreement with numerous methodologies, proving that the LRPSM is a useful substitute for approaches using He’s or Adomian polynomials. Additionally, error functions are used to compare the exact and approximative solutions graphically and numerically. The higher degrees of accuracy and convergence rates were confirmed by the error analysis, demonstrating the suggested method’s efficacy and reliability. The process is rapid, exact, and simple to use, and it produces excellent results.

The article’s structure is as follows: In Section 2, we use a number of important definitions and conclusions from the theory of FC. The main concept of LRPSM is examined in Section 3 in order to establish and determine the solution of the time-fractional BSOPE. Section 4 investigates the potential, capability, and simplicity of the suggested approach using three numerical models. In Section 5, graphics and tables are used to investigate the numerical outcomes and discussion. Section 6, towards the end, has the conclusion.

2. Preliminaries

In this section, we examine several common definitions, properties, and some useful consequences that we used in this article.

Definition 1 (see [24]). We assume that the function is of exponential order and piecewise continuous. Then, the LT of for is formulated asand the inverse LT is defined bywhere lies in the right half plane of the absolute convergence of the Laplace integral.

Lemma 1 (see [2429]). Let and be piecewise continuous on and be of exponential order. We assume that , , and , are constants. Then, the properties mentioned as follows are valid:(i)(ii)(iii)(iv)(v)

Definition 2 (see [41]). The fractional derivative of I order in the CD sense is defined as follows:where is the R-L integral of .

Theorem 1 (see [24]). We assume that the multiple fractional power series (MFPS) representation for the function is given bythen we havewhere .

The conditions for the convergence of the MFPS are determined in the following theorem.

Theorem 2 (see [24]). Let can be denoted as the new form of MFPS explained in Theorem 1. If , on with , then the remainder of the new form of MFPS satisfies the following inequality:

3. Methodology for the LRPSM for the Time-Fractional BSOPE

This section examines the steps for using the suggested method to find solutions to BSOPEs. Running the LT on the BSOPE and then considering MFPS as the BSOPE’s new space solution constitute the main idea of LRPSM. The way in which the coefficients of this series utilize the limit idea is the main difference between the LRPSM and the RPSM. The generated consequents are then transformed into real space using the inverse LT. The guidelines for using the LRPSM to find solutions are as follows:Step 1: Equation (1) should be changed as follows:Step 2: by considering LT on both sides of equation (8), we obtain the following:where andStep 3: we assume that the solution of equation (9) is the following series:Step 4: we obtained the following as a result of using Lemma 1(iii):Step 5: we define the -truncated expansion of asStep 6: we introduce Laplace residual function (LRF) of equation (9) and the -LRF, respectively, as follows:Step 7: we use the expansion form of into .Step 8: we multiply both sides of with .Step 9: by utilizing the fact in equation (16), we solve the following sequence of algebraic equations for , where , step by step:Step 10: we use the attained values of into the -truncated expansion of for each to attain the -approximate solution of the algebraic equation in equation (9).Step 11: we apply the inverse LT on the final form of to attain -approximate solution, of the suggested problem.

4. Some Illustrated Problems

In this section, three time-fractional BSOPE in the CD sense are solved in order to assess the effectiveness and suitability of the suggested approach.

Problem 1. consider the time-fractional BSOPE that follows [34]:under the following initial conditions:We will examine the case when . Applying LT on both sides of equation (18), we haveMaking use of the process outlined in Section 3, we get the findings from equation (20) as follows:We suppose that the expansion of is as follows:The th-truncated expansion is given as follows:We obtained the following as a result of using Lemma 1(iii):As a result, the -truncate expansion of equation (21) is as follows:The LRF of equation (21) is as follows:The -LRF of equation (21) is as follows:We expand the characteristics of the RPSM to highlight the following details [42, 43]:(i) and (ii)(iii), where . To determine the first unknown co-efficient in (25), we have to use the truncated series into the -LRF, , to obtain is used on both sides of equation (28):We use the fact thatAs a result, we obtained as follows:Similarly, to find out values of the second undefined coefficient , we have to use the -truncates series into the -LRF to obtainUsing on both sides of equation (32), we get the following equation:Again, we use the fact thatAs a result, we obtained the coefficient in the following form:Therefore, the -approximate LRPS solution of equation (21) isTypically, to find the coefficients , first we use the -truncated series in equation (25), then we utilize it in the -LRF, equation (27), we multiply by , and then we solve the algebraic equation as follows:We get the following results by utilizing the previous procedure:The approximate solution of equation (21) is obtained by five iterations as follows:By applying the inverse LT to equation (39), we are able to approximate the fifth step solution in the original feature space:When we use in equation (40), we get the following form:which are the first six terms of the expansion and, thus, is the exact solution of equations (18) and (19) at .

Problem 2. consider the following time-fractional BSOPE [35]:subject to the following initial conditions:First, we perform LT on both sides of equation (42), we use the initial condition from equation (43), and then we format the resulting equation as follows:We describe the expansion solution of the algebraic equation (44). So, we suppose that the series of is as follows:The -truncated series of the expansion of is as follows:We obtained the following as a result of using Lemma 1(iii):So, the -truncated expansion becomes as follows:The LRF of (44) is as follows:The -LRF of equation (44) is designed as follows:To determine the first unknown coefficient in equation (46), we have to use -truncated expansion into the -LRF , then we multiply by on both sides, and then we use the following fact to obtainSimilarly, to establish the value of the second undefined coefficient , we have to utilize the truncated expansion into the -LRF and use the following fact , we haveTherefore, the approximate solution derived from the iteration of equation (44) is as follows:To determine the , , and unknown coefficients, we repeat the same process. We get as follows:Therefore, the approximate solution derived from the iteration of equation (44) is as follows:By utilizing the inverse LT on both sides of equation (55), the approximate solution derived from the iteration by LRPSM of equations (42) and (43) is as follows:When is used in equation (56), we get as follows:As a result, the following is the exact solution of equations (42) and (43) for :

Problem 3. consider the following time-fractional BSOPE [36]:subject to the initial condition:First, we perform LT on both sides of equation (59), using the initial condition from equation (60), and then format the resulting equation as follows:We describe the expansion solution of the algebraic equation (61). Therefore, we assume that the expansion of is as follows:The -truncated expansion of equation (61) is as follows:The LRF of equation (61) takes the following form:Accordingly, the -LRF takes the following form:We substitute the -truncated series equations (63) into (64), multiply the resulting equation by , and then solve the equation , where for givesTherefore, the approximate solution derived from the iteration of equation (61) is as follows:By utilizing the inverse LT on both sides of equation (67), the approximate solution derived from the iteration by LRPSM of equations (59) and (60) is as follows:When is used in equation (68), we get as follows:As a result, the following is the exact solution of equations (59) and (60) for :

5. Numerical Simulation and Discussion

The findings of the results of the models presented in Problems 13 are evaluated graphically and numerically in this section. Error functions can be used to evaluate the numerical method’s correctness and competency. It is important to provide the errors of the approximation analytical solution that the LRPSM offers in terms of an infinite fractional power series. We used the absolute and recurrence error functions to demonstrate the accuracy and strength of LRPSM.

The 2D graphs of the comparative analysis of the exact and approximative solutions derived by the suggested method in Problems 1–3 are shown in Figures 1(a)1(c). These figures display the 2D plots of the precise solution and the approximate solution attained from the fifth iteration attained by LRPSM for Examples 1–3 when and 1.0 in the range of . These graphs indicate that when is applied, the approximative solution converges to the precise solution. The precise result and the approximation overlap at , demonstrating the efficacy and accuracy of the recommended method.

Figures 2(a)2(c) display the 2D curve of calculating the similarity with the help of absolute error of the approximation formed in the fifth step and the precise solution found by the suggested technique for Examples 1–3, respectively, for in the range of . The article has shown that the fifth-step approximation solution of the recommended approach is very close to the exact solution. By showing the absolute error of the precise and approximate outputs on a graph, LRPSM’s precision is proven.

The comparison article using the 3D curve is shown in Figures 3(a)3(c) in the sense of the absolute error of the approximate finding from the fifth iteration and the exact finding found using the suggested method to Examples 1–3, respectively, at in the ranges of and . The article has revealed that the fifth-step approximation of the recommended approach is very similar to the precise result. The absolute error of the precise and approximation findings on 3D graphs serves as a demonstration of the precision of LRPSM.

In Tables 16, the numerical convergence of the approximation to the precise solution has been demonstrated by and in the range . Tables 1–6 shows that and obtained by the suggested method quickly approaches the when . We can see from Tables 1, 3, and 5 that all of the test problems for the fourth stage have very low recurrence errors. The recurrence error will further decrease if we take into account the fifth-step approximation shown in Tables 2, 4, and 6. The approximation is rapidly approaching the exact solution as a result of the accuracy of our suggested strategy being demonstrated by the recurring error process. We arrived at the conclusion that the suggested approach is a feasible and efficient technique for solving particular classes of FODEs with fewer calculations and iteration steps.

For appropriately selected points, and in the range obtained by LRPSM at in Examples 1–3 are displayed in Tables 712 for comparison article in the sense of the absolute error of the approximate and the exact finding. We can see from Tables 7, 9, and 11 that approximate solutions derived from the 6th iteration for all of the test problems have very low errors. The absolute error will further decrease if we take into account the 7th-step approximation shown in Tables 8, 10, and 12. By quantitatively comparing the absolute inaccuracy of the precise and approximative findings, LRPSM’s precision is shown.

Tables 1315 also compares the absolute error of the approximations from the fifth iteration obtained by the LRPSM of Examples 1–3 at plausible short-listed grid points in the range with the absolute error of the fifth-step approximations obtained by various methods, including the PDTM [34], ATDM [35], and HPM [36]. Strong agreement between the results produced using the suggested method and various series solution techniques shows that the LRPSM is a useful substitute for approaches based on He’s or Adomian polynomials.

6. Conclusion

In this article, we used the LRPSM to solve time-fractional BSOPE in the sense of CD. The effectiveness of the LRPSM has been demonstrated by results in both graphs and numerically. The approximate solutions achieved using LRPSM are in perfect agreement with the corresponding exact solutions, as can be seen from the graphs and tables. The numerical evidence for the convergence of the approximative solution to the exact solution is presented in Tables 16. The comparison article is established in Tables 712 in terms of the absolute error of the approximate and exact solutions. The results of these numerical approaches are also compared with the PDTM, ATDM, and HPM in terms of absolute errors in Tables 1315.

The results obtained using the suggested approach, which demonstrates excellent agreement with PDTM, ATDM, and HPM, show that LRPSM is a suitable replacement tool for He’s or Adomian polynomial-based methods used to solve FODEs. When determining the coefficients for a series such as the RPSM, the fractional derivative needs to be determined each time. However, LRPSM just needs the idea of an infinite limit. As a result, the computations required to determine the coefficients are minimal. The results led us to the conclusion that our technique is easy to use, accurate, flexible, and effective.

Data Availability

No data were generated or analyzed during the current study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

The authors declare that the study was realized in collaboration with equal responsibility. All authors read and approved the final manuscript.