The multiterm fractional variable-order differential equation has a massive application in physics and engineering problems. Therefore, a numerical method is presented to solve a class of variable order fractional differential equations (FDEs) based on an operational matrix of shifted Chebyshev polynomials of the fourth kind. Utilizing the constructed operational matrix, the fundamental problem is reduced to an algebraic system of equations which can be solved numerically. The error estimate of the proposed method is studied. Finally, the accuracy, applicability, and validity of the suggested method are illustrated through several examples.