The variable-order fractional differential equations appear in modeling diverse physical problems. The main issue we address in this paper concerns an accurate numerical solution of a class of variable-order differential equations. The given problem is transformed into a system of algebraic equations using the so-called operational matrix of variable-order differentiation and the shifted LegendeGauss-Radau collocation approach. Accordingly, the effort performed in calculations can be reduced.
Numerical simulation for a specific problem is presented to demonstrate the computational efficiency
and accuracy of the proposed algorithm.