“…Many research works have been allocated to providing efficient numerical procedures for solving constant order fractional models. For example, a general formulation based on the Hamiltonian function associated with optimal control of fractional problems without delay (Agrawal 2004), a collocation method based on the Bessel functions (Tohidi and saberi 2015), two-dimensional Müntz-Legendre hybrid functions (Sabermahani et al 2020), Müntz-Legendre polynomials (Kheyrinataj and Nazemi 2020a), a hybrid of orthonormal Taylor polynomials (Marzban and Malakoutikhah 2019, a hybrid of the conventional Legendre polynomials (Marzban 2021a), combining fractional-order Legendre functions with the block-pulse functions (Marzban 2021b), fractional Chebyshev functions (Kheyrinataj and Nazemi 2020b), Genocchi polynimials (Chang et al 2018), a neural network scheme (Yavari and Nazemi 2019), Bernstein polynomials (Nemati 2018), two-dimensional Müntz-Legendre wavelets (Sabermahani 2020), Ritz’s method (Jahanshahi and Torres 2017), a hybrid method with the use of Hermite cubic spline multi-wavelets (Mohammadzadeh and Lakestani 2018), Bernoulli wavelets (Rahimkhani et al 2017), Euler–Lagrange equation (Rakhshan and Effati 2020), second Chebyshev wavelets technique (Baghani 2021), Legendre wavelet approach (Yuttanan et al 2021), Hartley series (Dadkhah and Mamehrashi 2021). The main idea and fundamental concepts of variable-order fractional operators have been extended and studied in the excellent works (Samko and Ross 1993) and (Lorenzo and Hartley 2002).…”