A Bernoulli Operational Matrix Method for Solving Nonlinear Multi-Term Fractional Variable-Order Delay Differential Equation

Document Type : Research Article

Authors

1 Department of Applied Mathematics, Imam Khomeini International University, Qazvin, 34148-96818, Iran

2 Department of Mathematical Sciences, Montana State University, Bozeman, MT 59717, USA

3 School of Computing, University of North Florida, Jacksonville, FL 32224, USA

Abstract

This paper investigates the generalized non-linear multi-term delay fractional differential equation of variable order. The Bernoulli operational matrix method is utilized to address a category of these equations, transforming the original problem into a system of algebraic equations amenable to numerical solutions. Sufficient and complete numerical tests are presented to showcase the accuracy, generality, and efficiency of the presented technique, as well as the adaptability of this approach. The numerical results of this method are compared with the exact solution. A comparison of this scheme's results with the exact solution illustrates the method's effectiveness and validity.

Keywords

Main Subjects


 

Article PDF

[1] D. Bojovic and J. Boško, Fractional order convergence rate estimates of finite difference method on nonuniform meshes, J. Comput. ´ Meth. Appl. Math., 1(3), 213–221, (2001).
[2] D. Baleanu, R. L. Magin, S. Bhalekar, and V. Daftardar-Gejji, Chaos in the fractional order nonlinear Bloch equation with delay, Commun. Nonlinear. Sci. Numer. Simul., 25(1-3), 41–49, (2015).
[3] K. Diethelm, N. J. Ford, and A. D. Freed, Detailed error analysis for a fractional Adams method, Numer. Algorithms, 36(1), 31–52, (2004).
[4] Y. Kuang, Delay differential equations: with applications in population dynamics, vol. 191, Academic Press, 1993.
[5] M. Saedshoar Heris and M. Javidi, On fractional backward differential formulas for fractional delay differential equations with periodic and anti-periodic conditions, Appl. Numer. Math., 118, 203–220, (2017).
[6] E. Shivanian, A. Jafarabadi, T. G. Chegini, and A. Dinmohammadi, Analysis of a time-dependent source function for the heat equation with nonlocal boundary conditions through a local meshless procedure, Computational and Applied Mathematics, 44(6), 282, (2025).
[7] N. Hosseinzadeh, E. Shivanian, M. Z. Fairooz, and T. G. Chegini, A robust RBF-FD technique combined with polynomial enhancements for valuing European options in jump-diffusion frameworks, International Journal of Dynamics and Control, 13(6), 212, (2025).
[8] C. Lubich, Discretized fractional calculus, SIAM J. Math. Anal., 17(3), 704–719, (1984).
[9] L. Galeone and R. Garrappa, On multistep methods for differential equations of fractional order, Mediterr. J. Math., 3(3-4), 565–580, (2006).
[10] R. Garrappa, Trapezoidal methods for fractional differential equations: theoretical and computational aspects, Math. Comput. Simul., 110, 96–112, (2015).
[11] B. Azarnavid, M. Nabati, and M. Emamjomeh, The iterative reproducing kernel Hilbert space method based on the Fibonacci polynomials for the nonlinear fractional differential equations with fractional integral conditions, Journal of Advanced Mathematical Modeling, 13(3), 344–354, (2023).
[12] V. Daftardar-Gejji, Y. Sukale, and S. Bhalekar, Solving fractional delay differential equations: a new approach, Frac. Calculus Appl. Anal., 18(2), 400–418, (2015).
[13] A. Jhinga and V. Daftardar-Gejji, A new numerical method for solving fractional delay differential equations, J. Comput. Appl. Math., 38, 166, (2019).
[14] M. Ghasemi, M. Fardi, and R. Khoshsiar Ghaziani, Numerical solution of nonlinear delay differential equations of fractional order in reproducing kernel Hilbert space, Appl. Math. Comput., 268, 815–831, (2015).
[15] M. L. Morgado, N. J. Ford, and P. Lima, Analysis and numerical methods for fractional differential equations with delay, J. Comput. Appl. Math., 252, 159–168, (2013).
[16] J. Cermák, J. Horní ˇ cek, and T. Kisela, Stability regions for fractional differential systems with a time delay, Commun. Nonlinear Sci. ˇ Numer. Simul., 31(1), 108–123, (2016).
[17] M. P. Lazarevic and A. M. Spasi ´ c, Finite-time stability analysis of fractional order time-delay systems: Gronwall’s approach, Math. ´ Comput. Model., 49(3), 475–481, (2009).
[18] V. Daftardar-Gejji and H. Jafari, An iterative method for solving non linear functional equations, J. Math. Anal. Appl., 316, 753–763, (2006).
[19] V. Daftardar-Gejji, Y. Sukale, and S. Bhalekar, A new predictor-corrector method for fractional differential equations, Appl. Math. Comput., 244, 158–182, (2014).
[20] S. Bhalekar and V. Daftardar-Gejji, A predictor-corrector scheme for solving non-linear delay differential equations of fractional order, J. Frac. Calculus Appl., 1(5), 1–8, (2011).
[21] Z. Wang, A numerical method for delayed fractional-order differential equations, J. Appl. Math., 2013, 256071, (2013).
[22] K. Diethelm, N. J. Ford, and A. D. Freed, A predictor-corrector approach for the numerical solution of fractional differential equations, Nonlinear Dynamics, 29, 3–22, (2002).
[23] K. Diethelm and N. J. Ford, Analysis of fractional differential equations, J. Math. Anal., 265(2), 229–248, (2002).
[24] H. R. Khodabandehlo, E. Shivanian, and S. Abbasbandy, Numerical solution of nonlinear delay differential equations of fractional variable-order using a novel shifted Jacobi operational matrix, Eng. Comp., 38(Suppl 3), S2593–S2607, (2022).
[25] H. R. Khodabandelo, E. Shivanian, and S. Abbasbandy, A novel shifted Jacobi operational matrix method for nonlinear multi-terms delay differential equations of fractional variable-order with periodic and anti-periodic conditions, Math. Meth. Appl. Sci., 45(16), 10116–10135, (2022).
[26] H. R. Khodabandehlo, E. Shivanian, and S. Abbasbandy, A Novel Shifted Jacobi Operational Matrix for Solution of Nonlinear Fractional Variable-Order Differential Equation with Proportional Delays, Int. J. Indust. Math., 14(4), 415–432, (2022).
[27] H. R. Khodabandehlo, E. Shivanian, and S. Abbasbandy, A novel shifted Jacobi operational matrix method for linear multi-terms delay differential equations of fractional variable-order with periodic and anti-periodic conditions, Kragujevac Journal of Mathematics, 50(1), 39–69, (2026).
[28] S. G. Samko and B. Ross, Integration and differentiation to a variable fractional order, Integr. Transf. Spec. Funct., 1, 277–300, (1993).
[29] T. Odzijewicz, A. B. Malinowska, and D. F. M. Torres, Noethers theorem for fractional variational problems of variable order, Cent. Eur. J. Phys., 11, 691–701, (2013).
[30] S. Chen, F. Liu, and K. Burrage, Numerical simulation of a new two-dimensional variable-order fractional percolation equation in non-homogeneous porous media, Comput. Math. Appl., 67, 1673–1681, (2014).
[31] D. Tavares, R. Almeida, and D. F. M. Torres, Caputo derivatives of fractional variable order: numerical approximations, Commun. Nonlinear Sci. Numer. Simul., 35, 69–87, (2016).
[32] J. Liu, X. Li, and L. Wu, An operational matrix of fractional differentiation of the second kind of Chebyshev polynomial for solving multi-term variable order fractional differential equation, Math. Probl. Eng., 2016, 7126080, (2016).
[33] A. H. Bhrawy, E. Tohidi, and F. Soleymani, A new Bernoulli matrix method for solving high-order linear and nonlinear Fredholm integro-differential equations with piecewise intervals, Appl. Math. Comput., 219, 482–497, (2012).
[34] B. Azarnavid, M. Emamjomeh, M. Nabati, and A. Dinmohammadi, An efficient iterative method for multi-order nonlinear fractional differential equations based on the integrated Bernoulli polynomials, Computational and Applied Mathematics, 43(1), 43–68, (2024).
[35] B. Azarnavid, The Bernoulli polynomials reproducing kernel method for nonlinear Volterra integro-differential equations of fractional order with convergence analysis, Computational and Applied Mathematics, 42(1), (2022).
[36] M. Fardi, M. A. Zaky, and B. Azarnavid, Fractional Convection-Diffusion Systems in Complex 2D and 3D Geometries: A Bernoulli Polynomial-Based Kernel Method, Applied Numerical Mathematics, 218, 275–297, (2025).
[37] E. Tohidi, A. H. Bhrawy, and K. A. Erfani, Collocation method based on Bernoulli operational matrix for numerical solution of generalized pantograph equation, Appl. Math. Model., 37, 4283–4294, (2013).
[38] F. Toutounian and E. Tohidi, A new Bernoulli matrix method for solving second order linear partial differential equations with the convergence analysis, Appl. Math. Comput., 223, 298–310, (2013).
[39] S. Bazm, Bernoulli polynomials for the numerical solution of some classes of linear and nonlinear integral equations, J. Comput. Appl. Math., 275, 44–60, (2015).
[40] S. Nemati, P. M. Lima, and D. F. M. Torres, Numerical Solution of Variable-Order Fractional Differential Equations Using Bernoulli Polynomials, Fractal Fract., 5, 219, (2019).
[41] E. Keshavarz, Y. Ordokhani, and M. A. Razzaghi, Numerical solution for fractional optimal control problems via Bernoulli polynomials, J. Vib. Control, 22, 3889–3903, (2016).
[42] A. H. Bhrawy and M. A. Zaky, A method based on the Jacobi tau approximation for solving multi-term time-space fractional partial differential equations, J. Comput. Phys., 281, 876–895, (2015).
[43] Y. M. Chen, L. Q. Liu, B. F. Li, and Y. N. Sun, Numerical solution for the variable-order linear cable equation with Bernstein polynomials, Appl. Math. Comput., 238, 329–341, (2014).
[44] F. Costabile, F. Dellaccio, and M. I. Gualtieri, A new approach to Bernoulli polynomials, Rend. Mat. Ser. VII, 26, 1–12, (2006).
[45] G. Arfken, Mathematical Methods for Physicists, 3rd ed., Academic Press, 1985.
[46] E. H. Doha, A. H. Bhrawy, and S. S. Ezz-Eldien, A new Jacobi operational matrix: an application for solving fractional differential equations, Appl. Math. Model., 36, 4931–4943, (2012).
[47] E. H. Doha, A. H. Bhrawy, and S. S. Ezz-Eldien, A Chebyshev spectral method based on operational matrix for initial and boundary value problems of fractional order, Comput. Math. Appl., 62, 2364–2373, (2011).
[48] S. Mashayekhi, Y. Ordokhani, and M. Razzaghi, Hybrid functions approach for nonlinear constrained optimal control problems, Commun. Nonlinear. Sci. Numer. Simulat., 17, 1831–1843, (2012).
[49] S. A. Yousefi and M. Behroozifar, Operational matrices of Bernstein polynomials and their applications, Inter. Systems Sci., 32, 709–716, (2010).
[50] W. Labecca, O. Guimaraes, and J. R. C. Piqueira, Dirac’s formalism combined with complex Fourier operational matrices to solve initial and boundary value problems, Commun. Nonlinear Sci. Numer. Simul., 19(8), 2614–2623, (2014).
[51] M. Razzaghi and S. Yousefi, Legendre wavelets method for the nonlinear Volterra-Fredholm integral equations, Math. Comput. Simul., 70, 1–8, (2005).
[52] H. Danfu and S. Xufeng, Numerical solution of integro-differential equations by using CAS wavelet operational matrix of integration, Appl. Math. Comput., 194, 460–466, (2007).
[53] S. H. Behiry, Solution of nonlinear Fredholm integro-differential equations using a hybrid of block pulse functions and normalized Bernstein polynomials, J. Comput. Appl. Math., 260, 258–265, (2014).
[54] A. Saadatmandi and M. Dehghan, A new operational matrix for solving fractional-order differential equations, Comput. Math. Appl., 59, 1326–1336, (2010).
[55] A. Saadatmandi, Bernstein operational matrix of fractional derivatives and its applications, Appl. Math. Model., 38, 1365–1372, (2014).
[56] M. H. Atabakzadeh, M. H. Akrami, and G. H. Erjaee, Chebyshev operational matrix method for solving multi-order fractional ordinary differential equations, Appl. Math. Model., 37, 8903–8911, (2013).
[57] A. H. Bhrawy and A. S. Alofi, The operational matrix of fractional integration for shifted Chebyshev polynomials, Appl. Math. Lett., 26, 25–31, (2013). 
Volume 10, Issue 1
December 2025
Pages 1-16
  • Receive Date: 08 September 2025
  • Revise Date: 15 October 2025
  • Accept Date: 20 October 2025
  • Publish Date: 17 December 2025