this paper presents discrete galerkin method for obtaining the numerical solution of higher even-order integro-differential equations with variable coefficients. we use the generalized jacobi polynomials with indexes corresponding to the number of homogeneous initial conditions as natural basis functions for the approximate solution. numerical results are presented to demonstrate the effectiven...