An Operational Spectral Framework for Estimating Solutions of Systems of Fractional Integro-Differential Equations with Variable-order Derivatives
Received: 09-11-2025
Revised: 08-02-2026
Accepted: 12-02-2026
Published in Issue 26-05-2026
Published Online: 23-02-2026
Copyright (c) 2025 Shima Mizani, Hamideh Ebrahimi, Mahmoud Shafiee (Author)

This work is licensed under a Creative Commons Attribution 4.0 International License.
Abstract
This paper proposes a novel and efficient spectral numerical framework for approximating solutions of Volterra–Fredholm fractional integro-differential equations with variable-order derivatives. The method is based on the use of normalized and shifted Vieta–Lucas orthogonal polynomials as basis functions to represent the unknown solution. New operational matrices are derived for both the Caputo fractional derivative and the Riemann–Liouville fractional integral operators, enabling the transformation of the original fractional system into a system of linear or nonlinear algebraic equations. The proposed scheme exhibits favorable properties such as simple implementation, strong numerical stability, and reliable convergence behavior. Rigorous theoretical results concerning the existence, uniqueness, and convergence of the approximate solution are established. Several numerical examples are presented to demonstrate the accuracy and computational efficiency of the proposed approach.
Keywords
- Volterra-Fredholm integro-differential equations,
- Vieta-Lucas polynomials,
- Variable-order fractional integral and derivatives operators,
- Operational matrix
References
- Podlubny, I. (1999). Fractional Differential Equations, Academic Press. SanDiago.
- EI-Shahed, M., Nieto, J., & Ahmed, A. (2017). Fractional-Order Model for Biocontrol of the Lesser Date Moth in Palm Trees and Its Discretization. Advances in Difference Equations, 295, 193-209. https://doi.org/10.1186/s13662-017-1349-1
- Li, M. (2018). Three Classes of Fractional Oscillators. Symmetry-Basel, 10(2), 40-56. https://doi.org/10.3390/sym10020040
- Li, M., (2020). Multi-Fractional Generalized Cauchy Process and Its Application to Teletraffic. Physica A: Statistical Mechanics and its Applications, 550, 105-127. https://doi.org/10.1016/j.physa.2019.123982
- Li, M., & Zhao, Y. L. (2018). A Fast Energy Conserving Finite Element Method for the Nonlinear Fractional Schrodinger Equation with Wave Operator. Applied Mathematics and Computation, 338, 758-773. https://doi.org/10.1016/j.amc.2018.06.010
- Li, M., Gu, X. M., Huang, C., Fei, M., & Zhang, G. (2018). A fast linearized Conservative Finite Element Method for the Strongly Coupled Nonlinear Fractional Schrodinger Equations, Computational Physics, 358, 256-282. https://doi.org/10.1016/j.jcp.2017.12.044.
- Li, M., Huang, C., & Zhao, Y. L. (2020). Fast conservative numerical Algorithm for the Coupled Fractional Klein-Gordon-Schrodinger Equation. Numerical Algorithms, 84, 1081-1119. https://doi.org/10.1007/s11075-019-00793-9.
- Mirzaee, F., Samadyar, N. (2019). Numerical solution of Time Fractional stochastic Korteweg-de Vries Equation via Implicit Meshless Approach. Iranian Journal of Science and Technology, 43: 2905-2912. https://doi.org/10.1007/s40995-019-00763-9
- Mirzaee, F., & Samadyar, N. (2019). Numerical Solution based on Two-dimensional Orthonormal Bernstein Polynomials for Solving Some Classes of Two-Dimensional Nonlinear Integral Equations of Fractional Order. Applied Mathematics and Computation, 344-345, 191-203. https://doi.org/10.1016/j.amc.2018.10.020
- Mirzaee, F., & Samadyar, N. (2021). Implicit Meshless Method to Solve 2D Fractional Stochastic Tricomi-Type Equation Defined on Irregular Domain Occurring in Fractal Transonic Flow. Numerical Methods for Partial Differential Equations, 37 1781-1799. https://doi.org/10.1002/num.22608
- Mirzaee, F., & Samadyar, N. (2021). Combination of Finite Difference Method and Meshless Method based on Radial Basis Functions to Solve Fractional Stochastic Advection-Diffusion Equations. Engineering with Computers,36 1673-1686. https://doi.org/10.1007/s00366-019-00789-y.
- Mirzaee, F., Sayevand, K., Rezaeei, S., & Samadyar, N. (2020). Finite Difference and Spline Approximation for Solving Fractional Stochastic Advection-Diffusion Equation. Iranian Journal of Science and Technology, 45, 607-617. https://doi.org/10.1007/s40995-020-01036-6
- Saemi, F., Ebrahimi, H., & Shafiee, M. (2020). An Effective Scheme for Solving System of Fractional Volterra-Fredholm Integro-Differential Equations based on the Müntz-Legendre Wavelets. Computational and Applied Mathematics, 374, 273-289. https://doi.org/10.1016/j.cam.2020.112773
- Wang, L., & Chen, Y. M. (2020). Shifted-Chebyshev-polynomial-based numerical algorithm for fractional order polymer visco-elastic rotating beam. Chaos Solutions and Fractals, 132, 109585-109601. https://doi.org/10.1016/j.chaos.2019.109585
- Xu, T., Chen, W., Lu, S., & Chen, H. (2016). An Implicit Midpoint Difference Scheme for the Fractional Ginzburg-Landau Equation. Computational Physics, 312, 31-49. https://doi.org/10.1016/j.jcp.2016.02.018
- Babaei, A., Jafari, H., & Banihashemi, S. (2020). Numerical Solution of Variable Order Fractional Nonlinear Quadratic integro-Differential Equations based on the Sixth-kind Chebyshev Collocation Method. Computational and Applied Mathematics, 377, 112908-112931. https://doi.org/10.1016/j.cam.2020.112908
- Babaei, A., Parsa Moghadam, B., Banihashemi, S. , & Tenreiro, J. (2020). Numerical Solution of Variable-Order Fractional Integro-Partial Differential Equations via Sinc Collocation Method based on Dingle and double Exponential Transformations. Communications in Nonlinear Science and Numerical Simulation, 82, 104985-105010.
- https://doi.org/10.1016/j.cnsns.2019.104985
- Hosseininia, M., Heydari, M. H., Maalek Ghaini, F. M., & Avazzadeh, Z. (2019). A Wavelet Method to Solve Nonlinear Variable-Order Time Fractional 2D Klein–Gordon Equation. Computers and Mathematics with Applications, 78, 3713-3730. https://doi.org/10.1016/j.camwa.2019.06.008
- Hosseininia, M., Heydari, M. H., & Avazzadeh, Z. (2021). Numerical Study of the Variable-Order Fractional Version of the Nonlinear Fourth-Order 2D Diffusion-Wave Equation via 2D Chebyshev Wavelats. Engineering with Computers, 37, 3319-3328. https://doi.org/10.1007/s00366-020-00995-z
- Xie, J., Huang, Q., & Zhao, F. (2017). Numerical Solution of Nonlinear Volterra-Fredholm-Hammerstein Integral Equations in Two-Dimensional Space Based on Block Pulse Functions. Computational and Applied Mathematics, 317, 565-572. https://doi.org/10.1016/j.cam.2016.12.028.
- Xie, J., & Yi, M. (2019). Numerical Research of Nonlinear System of Fractional Volterra-Fredholm Integral Equations via Block-Pulse Functions and Error Analysis. Computational and Applied Mathematics, 345, 159-167. https://doi.org/10.1016/j.cam.2018.06.008. Xie, J., Ren, Z., & Li, Y. (2019). Numerical Scheme for Solving System of Fractional Partial Differential Equations with Volterra-Type Integral Term Through Two-Dimensional Block-Pulse Functions. Numerical Methods for Partial Differential Equations, 35, 1890-1903. https://doi.org/10.1002/num.22383.
- Nagy, A. M., & EI-Sayed, A. A. (2019). An Accurate Numerical Technique for Solving Two-Dimensional Time Fractional Order Diffusion Equation. International Journal Modelling and Simulation, 39, 214-221. https://doi.org/10.1080/02286203.2019.1592285.
- Nagy, A. M., Sweilam, N. H., & EI-Sayed, A. A. (2018). New Operational Matrix for Solving Multi-Term Variable Order Fractional Differential Equations. Journal of The Computational and Nonlinear Dynamics, 13, 11001-11021. https://doi.org/10.1115/1.4037922.
- Sweilam, N. H., Khader, M., & Al-Bar, R. F. (2008). Homotopy Perturbation Method for linear and Nonlinear System of Fractional Integro-Differential Equations. International Journal of Modelling and Simulation, 1, 73-87. https://doi.org/10.1016/j.enganabound.2015.02.002
- Ganji, D., & Rajabi, A. (2006). Assessment of Homotopy-Perturbation and Perturbation Methods in Heat Radiation Equation. International Communications in Heat and Mass Transfer, 33, 391-400. https://doi.org/10.1016/j.icheatmasstransfer.2005.11.001
- Ganji, D., & Sadighi, A. (2006). Applications of He's Homotopy-Perturbation Method to Nonlinear Coupled Systems of Reaction-Diffusion Equations. International Journal of Nonlinear Sciences and Numerical Simulation, 7, 411-418. https://doi.org/10.1515/IJNSNS.2006.7.4.411
- Riahi Beni, N., (2024). Vieta-Lucas operational matrix Technique for Fractional variable-order Integro-Differential Equations. Journal of Mahani Mathematical Research, 14(2), 147-173. https://doi.org/10.22103/JMMR.2025.24074.1699
- Noor, Z. A., Talib, I., Abdljewad, T., & Alqudah, M. A. (2022). Numerical study of Caputo Fractional-order Differential Equations by Developing New Operational Matrices of Vieta-Lucas Polynomials. Fractal and Fractional, 6(2), 79. https://doi.org/10.3390/fractalfract6020079
- Agarwal, P., & El-Sayed, A. A. (2023). Operational Matrix Approach for Solving Variable-Order Fractional Integro-Differential Equations. Fractional Order Systems and Applications in Engineering, 301-317. https://doi.org/10.1016/b978-0-32-390953-2.00025-6
- Gokcen Kocer, E. (2016). Bivariate Vieta-Fibonacci and Bivariate Vieta-Lucas Polynomials. Mathematics, 12(04), 44-50. https://doi.org/10.9790/5728-1204024450
- Heydari, M. H., Avazzadeh, Z., & Razzaghi, M. (2021). Vieta-Lucas Polynomials for the Coupled Nonlinear Variable-Order Fractional Ginzburg-Landau Equation. Applied Numerical Mathematics, 165, 442-458. https://doi.org/10.1016/j.apnum.2021.03.007
- Agarwal, P., & EI-Sayed, A. A. (2020). Vieta-Lucas Polynomials for Solving a Fractional-Order Mathematical Physics Model. Advances in Difference Equations, 626, 754-772. https://doi.org/10.1186/s13662-020-03085-y
10.57647/ijm2c.2026.1604.24
