Abstract

In this paper, radial basis functions (RBFs) method was used to solve a fractional Black-Scholes-Schrodinger equation in an option pricing of financial problems. The RBFs method is applied in discretizing a spatial derivative process. The approximation of time fractional derivative is interpreted in the Caputo’s sense by a simple quadrature formula. This RBFs approach was theoretically proved with different problems of two numerical examples: time step arbitrage bubble case and time linear arbitrage bubble case. Then, the numerical results were compared with the semiclassical solution in case of fractional order close to 1. As a result, both numerical examples showed that the option prices from RBFs method satisfy the semiclassical solution.

1. Introduction

An option is one of the most important and popular financial derivatives in financial market. There are various types of mathematical model for option pricing. The Black-Scholes equation, introduced by Black and Scholes [1], provided an approximate description of underlying asset price behavior. This equation becomes popular in various kinds of studies such as economics, physics, and financial mathematics since it can be simply solved with a short time in conversion into the solutions. The Black-Scholes equation is a well-known financial model in option pricing which is constructed under strict assumptions. In fact, arbitrage exists in real financial markets; however, one of the key assumptions in this equation has no arbitrage. Thus, Classical Black-Scholes equation was extended for arbitrage possibilities by Contreras et al. [2]. The Black-Scholes equation with arbitrage can be interpreted using quantum mechanic’s view point in a sense of an imaginary time from Schrodinger equation of a free particle. Therefore, the Black-Scholes equation including arbitrage possibilities was proposed by Contreras et al. [3] which the equation was solved by the semiclassical method. Although the Black-Scholes-Schrodinger equation can be used to describe the analysis of option pricing in financial markets, this equation cannot be completely described in the physical meaning of the actual financial market.

For a decade, fractional differential equations (FDEs) have been further used as a tool to describe the phenomena in applied sciences and engineering. Problems related to biology, chemistry, physics, mechanics, and engineering (e.g., surface and subsurface hydrology [46], finance [7, 8], epidemiology [9], and ecology [10, 11]) were efficiently explained by fractional differential and integral equations. As a result, FDEs are more suitably compared to the integer-order models [12]. Moreover, the fractional derivative can be used to describe some occurrences that integer cannot. There are many different types of fractional derivative. A most popular of fractional derivatives was proposed by Caputo [13] that it was deeply explained referred to Oldham and Spanier [14], Miller and Ross [15], Podlubny [16], and Kilbas et al. [17]. Sometimes, FDEs are not suitable for some analytical methods; the numerical methods that contribute to an abundance of approaches are used to overcome this disadvantage. Finite difference and meshless method are commonly applied in numerical methods. Cen and Le [18] presented a numerical method based on central difference spatial discretization for a generalized Black-Scholes equation for option pricing. Song and Wang [19] also solved the put option pricing problem based on the fractional Black-Scholes equation by finite difference method. The results proved that applications of this technique are efficiency and less requirement for computational work to solve the fractional Black-Scholes equation. Kumar et al. [20] provide the numerical algorithm called homotopy perturbation and homotopy analysis method for time fractional Black-Scholes equation to solve a European option problem. Phaochoo et al. [21, 22] proposed a numerical method based on the meshless local Petrov-Galerkin (MLPG) to solve a Black-Scholes equation and fractional Black-Scholes equation via moving kriging interpolation for financial problems. Likewise, Phramrung et al. [23] applied the numerical method of meshless local Petrov-Galerkin (MLPG) to approximate the problem of the fractional HIV model. In addition, Cen et al. [24] applied the central difference spatial discretization for time fractional Black-Scholes equation. Chen et al. [25] introduced a new operator splitting method for numerical approach in an American option under fractional Black-Scholes model. Numerical scheme was also operated by Uddin and Taufiq [26] using radial kernels and Laplace transform to approximate the time fractional Black-Scholes model governing European options.

Radial basis functions (RBFs), firstly introduced in 1971, is a new technique for the numerical solution of partial differential equations (PDEs) in Hardy research [27]. The RBFs are high-dimensional and highly accurate meshless computational algorithm with a number of distinct advantages. It is widely applied in field of applied science and engineering such as diffusion equations [28], reconstruction of corrupted images [29, 30], and surface reconstruction [31]. The generation of a grid impacted is not required in the RBFs method because it is particularly efficient in solving such kind of free boundary and convection dominated problems [32]. Furthermore, the RBFs approximation technique is based on collocation in a set of scattered nodes. This method is independent with respect to the dimension of the space. However, there are few studies that applied the RBFs method in option pricing. Hon and Mao [33] proposed a radial basis functions (RBFs) method for solving options pricing model with RBFs interpolation by converting Black-Scholes equation into a system of ordinary differential equations (ODEs). Besides, Hon [34] combined the quasi-interpolation and RBFs method, called a quasi-radial basis function method, to solve the option pricing model. The result showed a high accuracy in the computations for European and American options. Zhang [35] applied a radial basis functions method for valuing options with multinomial tree approach. The study claimed that the RBFs method is highly efficient for both European option and American option.

In this study, the fractional Black-Scholes-Schrodinger equation is solved by using RBFs method for an option pricing. Spatial derivative was discretized through this method, and the approximation of time fractional derivative is interpreted in the Caputo’s senses by a simple quadrature formula. Then, numerical solution is compared with the semiclassical solution. The procedure in this study was demonstrated in Section 2, where the basic concept of the Black-Scholes-Schrodinger equation was briefly described as a problem formulation. Spatial discretization of the fractional Black-Scholes-Schrodinger equation was analyzed by using the RBFs method as it was shown in Section 3, and temporal discretization of fractional Black-Scholes-Schrodinger equation was investigated with a simple quadrature formula in Section 4. Then, the RBFs method was confirmed by stability analysis to ensure the suitability of this method in Section 5. After that, in Section 6, numerical solutions were examined and discussed to validate the proposed method to summarize the conclusion of this work in Section 7.

2. Problem Formulation

The Black-Scholes-Schrodinger equation is a quantum financial model which is used for analyzing fair prices of options in real financial market. This equation interprets the Black-Scholes equation with arbitrage possibilities in quantum mechanic’s view point in the senses of the Schrodinger equation. The Black-Scholes equation with arbitrage possibilities is transformed to the Black-Scholes-Schrodinger equation. The Black-Scholes equation with arbitrage possibilities in the domain is presented in the form of Equation (1). where is the option price at underlying asset price with time , represents an underlying asset price, represents time variable, represents the volatility of underlying asset price, is the risk free interest rate, is the expiration date, and the is called the arbitrage bubble function. In case of , Equation (1) is reduced to the original Black-Scholes equation with arbitrage possibilities. Consequently, the degeneration will occur in approximation when converges to 0. In order to solve this problem, the changing variable technique is applied by , to obtain

In 2010, a new variable, and , is first introduced by Contreras et al. [3]. Therefore, Equation (2) transforms to the Black-Scholes-Schrodinger equation in the domain as follows: where represents a wave function at time , is a potential function, , and is called a space variable.

However, the Black-Scholes-Schrodinger equation is not completely consistent with the actual financial market. Therefore, the fractional calculus is used to apply in the Black-Scholes-Schrodinger equation to describe occurrences in financial market especially in field of log-price probability and to specify the variability in prices. The fractional Black-Scholes-Schrodinger equation in the domain can be expressed as where is a fractional order, 0 It is also considered as a model parameter, in which each model will provide its solution. Equation (4) is reduced to the original Black-Scholes-Schrodinger equation when Furthermore, this study investigates the solution in each different value of , where it can be any number in (0, 1] and examine how it affects the solution of the model.

3. Spatial Discretization

In this section, the radial basis functions (RBFs) method is applied in the process of discretizing a space variable, because the RBFs interpolation is stable and accurate [36]. The RBFs interpolation formulation can be written as where is the RBFs, is an unknown coefficient at time, , and is the number of support nodes in the interpolation domain of point .

Firstly, Equation (5) is substituted into the fractional Black-Scholes-Schrodinger equation (Equation (4)). where and for each node

Secondly, Equation (7) can be written in the matrix form as where

From a spatial discretization of fractional Black-Scholes-Schrodinger equation, the ordinary differential equation system of each point in the space is obtained. Therefore, these systems are in form of time-dependent equation. In the next section, a simple quadrature formula is applied for discretization of time variable.

4. Temporal Discretization

In this section, the systems of ODEs from previous section will be discretized on time variable by a simple quadrature formula. First of all, the time fractional derivative, in Equation (4) is defined by Caputo’s viewpoint of order as where denote the gamma function.

The variable is defined as , where is a number of time step and is the step size of a time variable. A simple quadrature formula from [37] is applied as where and .

Hence, and the first-order approximation method for the computation of Caputo’s fractional derivative is given by

Applying Equation (12) in Equation (8),

For , Equation (13) can be rewritten as and for ,

The formula in Equation (15) is applied to approximate at the time level and also used in Equation (16) for . Then, Equation (16) is substituted into Equation (5) to obtain the solution at each time level . Equation (16) can be expressed as where .

5. Stability Analysis

In this section, we analyze the stability of the RBFs method. Let be a small perturbation at the time level , where is the exact solution and is an approximate solution. Therefore, the equation of error can be written as

Equation (17) would be stable if boundary of in is increased indefinitely by exists as a positive number, that , then . Hence, , . Consider Equation (18) in the case of and , which if for . In case of , this equation can be done by mathematical induction technique. Equation (18) is first rearranged as follows:

Take the norm on both sides of Equation (19) and then apply the triangle inequality to get

For , the first step, base case, of mathematical induction intends to give which is obviously fulfilled. The second step, the inductive case, is proved that for any positive integer. Supposed that and a term of is obtained. From an inequality (Equation (20)), it can be written as which completes the proof. The proof from Equation (22) shows that the errors made at each time level of calculation will be no more than the errors made in the initial step as long as . Since contains a part of matrix and , their parameter was chosen to satisfy the condition of .

6. Numerical Experiments and Results

This section consists of a time step arbitrage bubble and time linear arbitrage bubble cases obtained from [3]. Both cases are represented for confirming the accuracy of the proposed numerical method, since the fractional Black-Scholes-Schrodinger equation has no analytical solution. Therefore, the solution of fractional Black-Scholes-Schrodinger equation is verified by comparing with the semiclassical solution given by [3]. The semiclassical solution in the presence of a time-dependent arbitrage bubble can be computed as where is the arbitrage-free Black-Scholes solution for the specific option with contract, , and is the factor. The arbitrage bubble function, , takes part of the function. In this way, the function renormalizes the bare arbitrage-free Black-Scholes solution. The pure Black-Scholes solution is given by where is the normal distribution function and with a strike price, . The contract function, , is given by

The fractional Black-Scholes-Schrodinger equation in the domain is considered as

The initial and boundary conditions can be obtained by changing variables of the contract function (Equation (25)) and the analytical solution (Equation (23)). Therefore, the initial and boundary conditions are as follows: where and are the end points boundary of the spatial domain.

Example 1. Consider the following time step arbitrage bubble case of with factor condition is determined by

In this paper, a binary put option is analyzed. Since the semiclassical solution is an analytical solution, the result can be computed from any point on a mesh. However, the aim of this paper is to compare the RBFs solution with semiclassical solution. Accordingly, semiclassical solution that is on the same node with RBFs solution was chosen. A set of node is defined as and where is the number of support nodes in the interpolation domain of point , and is a number of time step. In Example 1, the parameter values are, , , , , , , and respectively. Figure 1 shows the wave function solved by RBFs method, where dots on direction represent computational node, direction represents a time level, and vertical direction represent a wave function. Figure 2 shows the option price solved by RBFs method and semiclassical method, where direction represents an underlying asset price, direction represents a time level, and vertical direction represents the option price. The option price from RBFs method is shown as dots, while the semiclassical solution is shown as a mesh. As it is shown in Figure 2, the dots almost overlapped with the mesh; as a result, the option price of fractional Black-Scholes-Schrodinger equation (in case of close to 1) from RBFs method satisfied the semiclassical solution. On the other hand, Figure 3(a) only shows the option price solved by RBFs method with the initial and boundary conditions; a graph of option price remains steady at the initial time. The option price decrease continuously over time in each time step and remains steady in the final of the graph. In Figure 3(b), the option price plotted at time 1.0, 0.9, 0.8, 0.7, and 0.6 shows no difference compared with the graph at time 0.5, 0.4, and 0.3 which show a huge different. Figure 4 shows the option price from RBFs method when the values are varied. The result shows that the option price decrease more rapidly with the lower values; for example, the option price at decrease with higher slope than the option price at . Moreover, the relative error was analyzed in this paper. The values of the relative error between the RBFs solution of fractional Black-Scholes-Schrodinger equation with the semiclassical solution is . It was found that the numerical result agrees as the semiclassical solution.

In the following example, the option price of the fractional model by RBFs method for time linear arbitrage bubble case is examined.

Example 2. Consider the following time linear arbitrage bubble case with factor condition is determined by

According to Example 1, a set of node is defined as and where is the number of support nodes in the interpolation domain of point , and is a number of time step. In Example 2, the parameter values are, , , , , , , and , respectively. Figure 5 shows the wave function solved by RBFs method, where dot on direction represents computational nodes, direction represents a time level, and vertical direction represents a wave function. Figure 6 shows the option price graph solved by RBFs method and semiclassical method, where direction represents an underlying asset price, direction represents a time level, and vertical direction represents the option price which is shown as dots and mesh similar to the graph in previous example. According to Example 1, the dots and mesh almost overlap. Therefore, the option prices of fractional Black-Scholes-Schrodinger equation (in case of close to 1) by using RBFs method satisfied the semiclassical solution. In contrast, Figure 7(a) only shows the option price solved by RBFs method, a graph of option price remains steady at the initial time. The option price decrease continuously over time in each time step and remains steady in the final of the graph. Figure 7(b) shows that the option price at times 1.0, 0.9, 0.8, 0.7, and 0.6 is almost no difference in each time, while at time 0.5and 0.4, the option price drops considerably. At time 0.5 and 0.4, the option price fall slightly at time 0.3, 0.2, 0.1, and 0. Figure 8 shows the option price in different values of . When, value decreases, the option price changes dramatically. In other hand, the option price changes slightly when the values of increases. According to Example 1, the values of the relative error between the RBFs solution of fractional Black-Scholes-Schrodinger equation with the semiclassical solution is . It was found that the numerical result agrees as the semiclassical solution.

From previous examples, it concluded that the results from RBFs method (in case of close to 1) satisfied the semiclassical solution. The results of the fractional model for various values show that when value is decreased, the option price is changing rapidly while when value increases close to 1, the solutions are changing slowly.

7. Conclusion

The numerical method in solving the fractional Black-Scholes-Schrodinger equation based on the radial basis functions (RBFs) method is successfully developed. The spatial derivative is discretized by the RBFs method. A simple quadrature formula is applied in terms of time fractional derivative discretizing. The Gaussian function is selected as a radial basis function in RBFs method. The numerical solutions of fractional Black-Scholes-Schrodinger equation based on RBFs method in case of fractional order close to 1 is compared with the semiclassical solution. The results presented as time step arbitrage bubble case and time linear arbitrage bubble case, which show that the option price from RBFs method satisfied the semiclassical solution. Therefore, RBFs method can be used to solve the fractional Black-Scholes-Schrodinger equation. Because the effect in changing the value of alpha on the fractional model is considered, the option price changes rapidly when the value of alpha is decreased. On the other hands, when alpha was increasing close to 1, the option price changed slowly. Therefore, it can be concluded that the varying fractional order affected to the model solutions. The value of the relative error between the RBFs solution of fractional Black-Scholes-Schrodinger equation with the semiclassical solution was analyzed in both example cases. It shows that the numerical result agrees as the semiclassical solution. Furthermore, the RBFs method is verified by the stability analysis. The results of solving the fractional Black-Scholes-Schrodinger equation by applied RBFs method can be ensured by stability analysis. Thus, this study provided an alternative useful approach in solving the fractional Black-Scholes-Schrodinger equation.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no competing interests.

Authors’ Contributions

The idea of this research was introduced by AL and NN. All authors contributed to the main results and numerical simulations. AL contributed to revising the manuscript. All authors read and approved the final manuscript.

Acknowledgments

The authors wish to express their gratitude to the editors and the reviewers for the helpful comments. This research is supported by the Department of Mathematics, Faculty of Science, King Mongkut’s University of Technology Thonburi (KMUTT) and Uttaradit Rajabhat University.