In this paper, a numerical method which produces an approximate polynomial solution is presented for solving the high-order linear singular differential-difference equations. With the aid of Bessel polynomials and collocation ...