Taylor polynomial solution of difference equation with constant coefficients via time scales calculus

In this study, we present a practical matrix method to find an approximate solution of higher order linear difference equation with constant coefficients under the initial-boundary conditions in terms of Taylor polynomials. To obtain this goal, we first present time scale extension of previous polynomial approach, then restrict the formula to the Integers with h step. This method converts the difference equation to a matrix equation, which may be considered as a system of linear algebraic equations.