“…Equation (1.1) includes many important kinds of equations (see for example [4,6,12,13,15,18,25,29]) and this method can be used to obtain numerical solutions of high order differential equations (for {k 1,v } k−1 v=0 = {k 2,v } k−1 v=0 = {M v } k−1 v=0 = 0), high order integro-differential equations (for {k 2,v } k−1 v=0 = {M v } k−1 v=0 = 0), high order delay differential equation (for {k 1,v } k−1 v=0 = {k 2,v } k−1 v=0 = 0). There are many existing numerical methods for solving Volterra integrodifferential equations, such as Legendre spectral collocation method [26], Runge-Kutta method [7], spectral method [27,28], Polynomial collocation method [8][9][10]23], Tau method [14], operational matrices [2], Homotopy perturbation method [16,24], Haar wavelet method [21], Taylor polynomial [20,22]. The aim of the present paper is to apply a direct collocation method based on the use of Taylor polynomials so that we generalize the Taylor collocation method for delay VIEs [3], first and second order delay IDEs [4,5], and high order VIDEs [19].…”