In this paper, a numerical matrix method, which is based on collocation points, is presented for the approximate solution of a system of high-order linear differential-difference equations with variable coefficients under ...