This study is aimed to develop a new matrix method, which is used an alternative numerical method to the other method for the high-order linear Fredholm integro-differential-difference equation with variable coefficients. This matrix method is based on orthogonal Jacobi polynomials and using collocation points. The improved Jacobi polynomial solution is obtained by summing up the basic Jacobi polynomial solution and the error estimation function. By comparing the results, it is shown that the improved Jacobi polynomial solution gives better results than the direct Jacobi polynomial solution, and also, than some other known methods. The advantage of this method is that Jacobi polynomials comprise all of the Legendre, Chebyshev, and Gegenbauer polynomials and, therefore, is the comprehensive polynomial solution technique