Abstract：In order to solve a class of variable fractional order nonlinear differential-integral equations, the numerical solution is proposed. Function approximation based on Shifted Jacobi polynomials, combining Caputo-type variable order fractional derivate definition, which are used to get the operational matrixes of Shifted Jacobi polynomials, are the main characteristic behind this method. With the operational matrix, the original equation is translated into the products of several dependent matrixes, which can be regarded as a system of nonlinear equations after dispersing the variable. By solving the nonlinear system of algebraic equations, the coefficients of Shifted Jacobi polynomials are got, then the numerical solutions of the original equation are acquired. Finally, some numerical examples illustrate the accuracy and effectiveness of the method.