A High-Order Numerical Method Based on a Spatial Compact Exponential Scheme for Solving the Time-Fractional Black–Scholes Model

This paper investigates a high-order numerical method based on a spatial compact exponential scheme for solving the time-fractional Black–Scholes model. Firstly, the original time-fractional Black–Scholes model is converted into an equivalent time-fractional advection–diffusion reaction model by means of a variable transformation technique. Secondly, a novel high-order numerical method is constructed with (2−α) accuracy in time and fourth-order accuracy in space based on a spatial compact exponential scheme, where α is the fractional derivative. The uniqueness of solvability of the derived numerical method is rigorously discussed. Thirdly, the unconditional stability and convergence of the derived numerical method are rigorously analyzed using the Fourier analysis technique. Finally, numerical examples are presented to test the effectiveness of the derived numerical method. The proposed numerical method is also applied to solve the time-fractional Black–Scholes model, whose exact analytical solution is unknown; numerical results are illustrated graphically.