In this paper, we present numerical method for solving integro-differential equations of fractional order based on a hybrid of block-pulse functions and Taylor polynomials. Fractional derivative is described in the Caputo sense. Some numerical examples are presented to demonstrate the theoretical results.
