<abstract><p>In this paper, we introduce a mixed numerical technique for solving fractional differential equations (FDEs) by combining Chebyshev collocation methods and a piecewise quadratic quadrature rule. For getting solutions at each integration step, the fractional integration is calculated in two intervals-all previous time intervals and the current time integration step. The solution at the current integration step is calculated by using Chebyshev interpolating polynomials. To remove a singularity which belongs originally to the FDEs, Lagrangian interpolating technique is considered since the Chebyshev interpolating polynomial can be rewritten as a Lagrangian interpolating form. Moreover, for calculating the fractional integral on the whole previous time intervals, a piecewise quadratic quadrature technique is applied to get higher accuracy. Several numerical experiments demonstrate the efficiency of the proposed method and show numerically convergence orders for both linear and nonlinear cases.</p></abstract>