“…Some methods for determining the numerical solution of integro-differential equations include: Bernstein Method [14], Adomian decompositions method [2,3], Finite difference-Simpson method [17], Collocation method by [4,5,6,7,21,22], Hybrid linear multistep method [8,9], Chebyshev-Galerkin method [10], Bernoulli matrix method [11], Differential transform method [12], Lagrange Interpolation [13], Differential Transformation [15], Block pulse functions operational matrices [19] Chebyshev polynomials [16], Optimal Auxiliary Function Method (OAFM) [18] and Spectral Homotopy Analysis Method [20]. We consider first order Volterra integro-differential equation of the form 𝑦 ′ (𝑥) = 𝑔(𝑥) + ∫ 𝑘(𝑥, 𝑡)𝑦(𝑡)𝑑𝑡 𝑥 0 (1) with the initial condition 𝑦(0) = 𝑞…”