Implicit cubic B-spline scheme for the fractional Black-Scholes model with Caputo-Hadamard derivative

Document Type : Research Paper

Author

Department of Mathematics, Payame Noor University (PNU), P.O.Box 19395-4697, Tehran, Iran

Abstract

In this study, we introduce a novel numerical scheme for solving the Black–Scholes equation endowed with a Caputo-Hadamard fractional time derivative. The temporal derivative is discretized via a finite-difference approach, ensuring both stability and high accuracy in capturing the model’s dynamics. For spatial discretization, we employ a cubic B-spline quasi-interpolation technique, which delivers a robust and efficient framework. A rigorous error analysis is carried out, establishing theoretical bounds that confirm the method’s reliability. Finally, we validate the proposed scheme through a series of numerical experiments, demonstrating its effectiveness in practical financial modeling and its ability to address the challenges posed by fractional differential equations.

Keywords

Main Subjects


[1] Abdi, N., Aminikhah, H., & Sheikhani, A. R. (2022). High-order compact  nite di erence schemes for the time-fractional Black-Scholes model governing European options. Chaos, Solitons & Fractals, 162, 112423. https://doi.org/10.1016/j.chaos.2022.112423
[2] Akram, T., Abbas, M., Abualnaja, K. M., Iqbal, A., & Majeed, A. (2022). An ecient numerical technique based on the extended cubic B-spline functions for solving time fractional Black-Scholes model. Engineering with Computers, 38(Suppl 2), 1705-1716. https://doi.org/10.1007/s00366-021-01436-1
[3] Alavi, J., & Aminikhah, H. (2021). Numerical study of the boundary layer problem over a at plate by orthogonal cubic spline basis functions. The European Physical Journal Plus, 136(7), 780. https://doi.org/10.1140/epjp/s13360-021-01788-z
[4] Alavi, J., & Aminikhah, H. (2021). Orthogonal cubic spline basis and its applications to a partial integro-di erential equation with a weakly singular kernel. Computational and Applied Mathematics, 40(2), 55. https://doi.org/10.1007/s40314-021-01442-5
[5] Alazemi, F., Alsena , A., & Naja , A. (2025). A spectral approach using fractional Jaiswal functions to solve the mixed time-fractional Black-Scholes European option pricing model with error analysis. Numerical Algorithms, 98(1), 347-371.
https://doi.org/10.1007/s11075-024-01797-w
[6] Aljethi, R. A., & Klcman, A. (2023). Analysis of fractional di erential equation and its application to realistic data. Chaos, Solitons & Fractals, 171, 113446. https://doi.org/10.1016/j.chaos.2023.113446
[7] Aminikhah, H., & Alavi, J. (2016). Numerical solution of convection-di usion equation using cubic B-spline quasi-interpolation. Thai Journal of Mathematics, 14(3), 599-613. http://thaijmath.in.cmu.ac.th/index.php/thaijmath/article/view/255
[8] Aminikhah, H., & Alavi, J. (2017). B-spline collocation and quasi-interpolation methods for boundary layer 
ow and convection heat transfer over a at plate. Calcolo, 54, 299-317. https://doi.org/10.1007/s10092-016-0188-x
[9] Aminikhah, H., & Alavi, S. J. (2019). A new approach to using the cubic B-spline functions to solve the Black-Scholes equation. https://sanad.iau.ir/fa/Article/797725 [In Persian]
[10] Bian, L., & Li, Z. (2021). Fuzzy simulation of European option pricing using sub-fractional Brownian motion. Chaos, Solitons & Fractals, 153, 111442. https://doi.org/10.1016/j.chaos.2021.111442
[11] Black, F., & Scholes, M. (1973). The pricing of options and corporate liabilities. Journal of Political Economy, 81(3), 637-654. https://www.journals.uchicago.edu/doi/abs/10.1086/260062
[12] Deng, B., Zhang, J., & Zhang, Z. (2019). Superconvergence points of integer and fractional derivatives of special Hermite interpolations and its applications in solving FDEs. ESAIM: Mathematical Modelling and Numerical Analysis, 53(3), 1061-1082. https://doi.org/10.1051/m2an/2019012
[13] Diethelm, K., & Ford, N. J. (2002). Analysis of fractional di erential equations. Journal of Mathematical Analysis and Applications, 265(2), 229-248. https://doi.org/10.1006/jmaa.2000.7194
[14] Guo, B., Pu, X., & Huang, F. (2015). Fractional Partial Di erential Equations And Their Numerical Solutions. World Scienti c.
[15] Karami, S., Heydari, M. H., & Baleanu, D. (2024). Fractional logarithmic Chebyshev cardinal functions: Application for Caputo-Hadamard fractional optimal control problems. Journal of Vibration and Control, 10775463241272962.
https://doi.org/10.1177/10775463241272962
[16] Lee, Y., Kim, Y., & Lee, J. (2020). Pricing various types of power options under stochastic volatility. Symmetry, 12(11), 1911. https://doi.org/10.3390/sym12111911
[17] Liu, Y., Zhu, Y., & Lu, Z. (2021). On Caputo-Hadamard uncertain fractional di erential equations. Chaos, Solitons & Fractals, 146, 110894. https://doi.org/10.1016/j.chaos.2021.110894
[18] Nagy, A. M., Makhlouf, A. B., Alsena , A., & Alazemi, F. (2021). Combination synchronization of fractional systems involving the Caputo-Hadamard derivative. Mathematics, 9(21), 2781. https://doi.org/10.3390/math9212781
[19] Nikan, O., Rashidinia, J., & Jafari, H. (2025). Numerically pricing American and European options using a time fractional Black-Scholes model in  nancial decision-making. Alexandria Engineering Journal, 112, 235-245. https://doi.org/10.1016/j.aej.2024.10.083
[20] Pan, G., & Zhang, S. (2023). A meshless multiquadric quasi-interpolation method for time fractional Black-Scholes model. International Journal of Financial Engineering, 10(02), 2350008. https://doi.org/10.1142/S2424786323500081
[21] Podlubny, I. (1998). Fractional Di erential Equations: An Introduction To Fractional Derivatives, Fractional Di erential Equations, To Methods Of Their Solution And Some Of Their Applications (Vol. 198). Elsevier.
[22] Samko, S. G. (1993). Fractional Integrals And Derivatives: Theory And Applications. Gordan and Breach Science.
[23] Zhang, H., Zhang, M., Liu, F., & Shen, M. (2024). Review of the fractional Black-Scholes equations and their solution techniques. Fractal and Fractional, 8(2), 101. https://doi.org/10.3390/fractalfract8020101