Fractional Black Scholes Model Research Articles

Series form solutions of time–space fractional Black–Scholes model via extended He-Aboodh algorithm

The objective of the current study is analyze linear and nonlinear time–space fractional Black–Scholes models via modified homotopy perturbation method (m-HPM). In current investigation, memory effects in financial markets are explored through fractional derivative in Caputo sense. The effectiveness of proposed methodology is checked numerically by finding residual errors and presented in tables. These tables also provide a benchmark for the comparison with already existing results in literature. Furthermore, solutions are graphically analyzed via 3D and contour plots across a range of parameters under varying market conditions. Analysis confirms the efficiency of m-HPM for predicting solutions of time–space fractional Black–Scholes models. The current study can contribute in understanding the applications of fractional calculus in finance, and can be a valuable computational tool for pricing financial derivatives in fractional environments.

Stability and Convergence Analysis of a Numerical Method for Solving a $$\zeta$$-Caputo Time Fractional Black–Scholes Model via European Options

Stability and Convergence Analysis of a Numerical Method for Solving a $$\zeta$$-Caputo Time Fractional Black–Scholes Model via European Options

Galerkin-finite difference method for fractional parabolic partial differential equations

The fractional form of the classical diffusion equation embodies the super-diffusive and sub-diffusive characteristics of any flow, depending on the fractional order. This study aims to approximate the solution of parabolic partial differential equations of fractional order in time and space. For this, firstly, we briefly discussed on some existing methods to solve Partial Differential Equations (PDEs) and fractional Differential Equations (DEs), then introduce a combined technique such as the Galerkin weighted residual method for the space fractional term with modified Bernoulli polynomials as basis functions, and the finite difference approximation for the time fractional term, respectively. The mathematical formulation of the proposed method is explained elaborately. Then we describe the order of convergence for the time fractional term only, as the convergence of the Galerkin method is obvious. We impose this technique on the fractional Black-Scholes model subsequently. Finally, we experimented with our proposed technique on some numerical problems. All the results are depicted in both tabular and 3D visualizations as well. We compare our results with the available methods in the literature, and our accuracy is considerable. To summarize: •The paper introduces an approach that integrates the Galerkin weighted residual method with modified Bernoulli polynomials to handle space fractional terms, alongside employing a finite difference approximation for time fractional terms.•The convergence analysis is focused.•The technique is implemented on the fractional Black-Scholes model and other numerical problems, with outcomes depicted through tables and 3D visualizations.

Numerical solutions for the time fractional Black-Scholes model governing European option by using double integral transform decomposition method

Numerical solutions for the time fractional Black-Scholes model governing European option by using double integral transform decomposition method

Analytical solution of time fractional Black-Scholes equation with two assets through new Sumudu Transform iterative method

There is a scopious rise in the study of financial derivatives over the past two or three decades. Mathematical model proposed by Black and Scholes expounds financial derivatives in a more momentous way. The Black-Scholes model on a single asset is a partial differential equation characterizing the behavior of European options. In this article, we introduce the new Sumudu transform iterative method (NSTIM) as a new technique to obtain the analytical solution of time fractional Black-Scholes model involving European options with two assets. The proposed model is the advanced version of the regular Black-Scholes model. Explicit solution of the problem has been obtained with the help of generalized Mittag-Leffer function. The numerical analysis prove that this method is efficacious in solving various problems of financial theory.

PRICING AMERICAN OPTION USING A MODIFIED FRACTIONAL BLACK–SCHOLES MODEL UNDER MULTI-STATE REGIME SWITCHING

The American option pricing problem is examined in this work using a regime switching finite moment log-stable model. The option prices under this model are governed by a coupled system of fractional partial differential equations. Combination of the coupled system and the spatial-fractional derivative makes it extremely difficult to find an analytic solution. We have constructed a numerical algorithm to numerically solve such problems. The developed predictor-corrector type method is highly efficient and reliable in solving coupled system in each regime having different volatility and interest rates. Two-sided Riesz space fractional diffusion term is approximated using fractional finite difference scheme whereas the classical space derivative term is approximated using central difference formula. Splitting technique is utilized to construct a highly efficient scheme which can also be implemented on parallel processors. Stability and error analysis of the scheme is proved analytically and demonstrated through numerical experiments. Effect of the order of the fractional derivative (also called tail index) on the option prices is shown through graphs by performing numerical experiments for different values of the tail index.

Analytical and Numerical Solution for the Time Fractional Black-Scholes Model Under Jump-Diffusion

In this work, we study the numerical solution for time fractional Black-Scholes model under jump-diffusion involving a Caputo differential operator. For simplicity of the analysis, the model problem is converted into a time fractional partial integro-differential equation with a Fredholm integral operator. The L1 discretization is introduced on a graded mesh to approximate the temporal derivative. A second order central difference scheme is used to replace the spatial derivatives and the composite trapezoidal approximation is employed to discretize the integral part. The stability results for the proposed numerical scheme are derived with a sharp error estimation. A rigorous analysis proves that the optimal rate of convergence is obtained for a suitable choice of the grading parameter. Further, we introduce the Adomian decomposition method to find out an analytical approximate solution of the given model and the results are compared with the numerical solutions. The main advantage of the fully discretized numerical method is that it not only resolves the initial singularity occurred due to the presence of the fractional operator, but it also gives a higher rate of convergence compared to the uniform mesh. On the other hand, the Adomian decomposition method gives the analytical solution as well as a numerical approximation of the solution which does not involve any mesh discretization. Furthermore, the method does not require a large amount of computer memory and is free of rounding errors. Some experiments are performed for both methods and it is shown that the results agree well with the theoretical findings. In addition, the proposed schemes are investigated on numerous European option pricing jump-diffusion models such as Merton’s jump-diffusion and Kou’s jump-diffusion for both European call and put options.

A meshless multiquadric quasi-interpolation method for time fractional Black–Scholes model

Based on multiquadric quasi-interpolation, this study presents a meshless numerical method for time fractional Black–Scholes (B–S) model. The method is highly accurate and flexible. The stability and convergence of the proposed scheme are discussed in detail. Numerical results show the accuracy and efficiency of the proposed method.

Numerical Investigation of Fractional Step-Down ELS Option

In this paper, we use the finite difference methods to explore step-down Equity Linked Securities (ELS) options under the fractional Black-Scholes model. We establish Crank-Nicolson scheme under one asset and study the impact of Hurst exponent (H) on return of repayment under fixed stock price. We also explore the impact of stock price on return of repayment under different H. Through numerical experiments, it is found that the return of repayment of options is related to H, and the result of difference scheme will increase with the increase of H. In the case of two assets, we establish implicit scheme, and in the case of three assets, we use operator splitting method (OSM) method to establish semi-implicit scheme. We get the result that the H also influences the return of repayment in two and three assets. We also conduct Greeks analysis. Through Greeks analysis, we find that the long-term correlation of stocks has a huge impact on investment gains or losses. Therefore, we take historical volatility (fractal exponents) into account which can significantly reduce risk and increase revenue for investors.

Error and stability estimates of a time-fractional option pricing model under fully spatial–temporal graded meshes

To price vanilla European and American options via the fractional Black–Scholes model, first a (2−α)-order discretization scheme for the Caputo fractional derivative based upon graded meshes along time is presented. This is fruitful for problems having nonsmooth data at the initial time. Second, a nonuniform discretization of space is also considered and the coefficient weights based upon meshless radial basis function generated finite difference (RBF-FD) methodology are contributed under the inverse quadratic function. The fully space–time nonuniform solution method is then constructed via the presented discretization formulas. Numerical tests are given to show effectiveness and accuracy of our numerical scheme.

Numerical Method for American Option Pricing under the Time‐Fractional Black–Scholes Model

The fractional Black–Scholes model has had limited applications in financial markets. Instead, the time‐fractional Black–Scholes equation has attracted much research interest. However, it is difficult to obtain the analytic expression for American option pricing under the time‐fractional Black–Scholes model. This paper will present an operator‐splitting method to price the American options under the time‐fractional Black–Scholes model. The fractional partial differential complementarity problem (FPDCP) that the American option price satisfied is split into two subproblems: a linear boundary value problem and an algebraic system. A high‐order compact (HOC) scheme and a grid stretching (GS) method are considered for the linear boundary problem. Furthermore, numerical results show that the HOC scheme with a GS method gives an accurate numerical solution for American options under the time‐fractional Black–Scholes model.

Parameter Estimation for a Fractional Black–Scholes Model with Jumps from Discrete Time Observations

We consider a stochastic differential equation (SDE) governed by a fractional Brownian motion (BtH) and a Poisson process (Nt) associated with a stochastic process (At) such that: dXt=μXtdt+σXtdBtH+AtXt−dNt,X0=x0>0. The solution of this SDE is analyzed and properties of its trajectories are presented. Estimators of the model parameters are proposed when the observations are carried out in discrete time. Some convergence properties of these estimators are provided according to conditions concerning the value of the Hurst index and the nonequidistance of the observation dates.

Inverse Multiquadric Function to Price Financial Options under the Fractional Black–Scholes Model

The inverse multiquadric radial basis function (RBF), which is one of the most important functions in the theory of RBFs, is employed on an adaptive mesh of points for pricing a fractional Black–Scholes partial differential equation (PDE) based on the modified RL derivative. To solve this problem, discretization along space is carried out on a non-uniform grid in order to focus on the hot area, at which the initial condition of the pricing model, i.e., the payoff, has discontinuity. The L1 scheme having the convergence order 2−α is used along the time fractional variable. Then, our proposed numerical method is built by matrices of differentiations to be as efficient as possible. Computational pieces of evidence are brought forward to uphold the theoretical discussions and show how the presented method is efficient in contrast to the exiting solvers.

Application of Lie symmetries to solving modified Black-Scholes equation

We perform Lie symmetry analysis on the modified Black-Scholes model described as a partial differential equation (PDE). As a result, a new complete Lie symmetry group and infinitesimal generators of the one-dimensional modified fractional Black-Scholes model are derived. Furthermore, we compute a family of exact invariant solutions that constitute the modified fractional Black-Scholes model using the associated infinitesimal generators and the corresponding similarity reduction equations. Using known solutions, more solutions are generated via group point transformations.

An Efficient Method for Solving Fractional Black-Scholes Model with Index and Exponential Decay Kernels

The Black-Scholes equation (BSe) is fascinating in the business world for predicting the performance of financial investment valuation systems. The Caputo fractional derivative (CFD) and Caputo-Fabrizio fractional derivative operators are used in this research to analyze the BSe. The Adomian decomposition method (ADM) and the new iterative transform (NIM) approach are combined alongside the Yang transform. In addition, the convergence and uniqueness results for the aforementioned framework have been calculated. The existence and uniqueness results have been established and frequently accompanied innovative aspects of the prospective system in fixed point terminologies. To provide additional insight into such concepts, a variety of illustrations and tabulations are used. Additionally, the provided techniques regulate and modify the obtained analytical results in a really productive fashion, allowing us to modify and regulate the converging domains of the series solution in a pragmatic manner.

Shadow Price Approximation for the Fractional Black Scholes Model

In this work, we used Tran Hung Thao’s approximation of fractional Brownian motion to approximate the shadow price of the fractional Black Scholes model. In the case to maximize expectation of the utility function in a portfolio optimization problem under transaction cost, the shadow price is approximated by a Markovian process and semimartingale.

Numerical Solution of Time Fractional Black–Scholes Model Based on Legendre Wavelet Neural Network with Extreme Learning Machine

In this paper, the Legendre wavelet neural network with extreme learning machine is proposed for the numerical solution of the time fractional Black–Scholes model. In this way, the operational matrix of the fractional derivative based on the two-dimensional Legendre wavelet is derived and employed to solve the European options pricing problem. This scheme converts this problem into the calculation of a set of algebraic equations. The Legendre wavelet neural network is constructed; meanwhile, the extreme learning machine algorithm is adopted to speed up the learning rate and avoid the over-fitting problem. In order to evaluate the performance of this scheme, a comparative study with the implicit differential method is constructed to validate its feasibility and effectiveness. Experimental results illustrate that this scheme offers a satisfactory numerical solution compared to the benchmark method.

Spectral Solutions for Fractional Black–Scholes Equations

This paper presents a numerical method to solve accurately the fractional Black–Scholes model of pricing evolution. A fully spectral collocation technique for the two independent variables is derived. The shifted fractional Jacobi–Gauss–Radau and shifted fractional Jacobi–Gauss–Lobatto collocation techniques are utilized. Firstly, the independent variables are interpolated at the shifted fractional Jacobi nodes, and the solution of the model is approximated by means of a sequence of shifted fractional Jacobi orthogonal functions. Then, the residuals at the shifted fractional Jacobi quadrature locations are estimated. As a result, an algebraic system of equations is obtained that can be solved using any appropriate approach. The accuracy of the proposed method is demonstrated using two numerical examples. It is observed that the new technique is more accurate, efficient, and feasible than other approaches reported in the literature. Indeed, the results show the exponential convergence of the method, both for smooth and nonsmooth solutions.

Pricing European Double Barrier Option with Moving Barriers Under a Fractional Black–Scholes Model

Pricing European Double Barrier Option with Moving Barriers Under a Fractional Black–Scholes Model

Numerical Valuation of European and American Options under Fractional Black-Scholes Model

In this paper, we investigate the numerical valuation of European and American options under the time fractional Black-Scholes model. We first apply a coordinate stretching transformation to the asset price so that the spatial region can focus on the vicinity of singularities, which are usually found in the payoff function. The radial basis function finite difference method is used for the spatial discretization, and the improved L1 method is used to deal with the reduced order of convergence for the nonsmooth initial data. We use the operator splitting method for solving the linear complementary problem of American options. The proposed scheme leads to a sparse linear system which is trivial to solve. Moreover, the stability of the proposed numerical scheme is analyzed using Fourier analysis. Numerical experiments demonstrate the accuracy and efficiency of the proposed method.