Discrete or digital systems that are only known at discrete moments are very important in various fields of physics and engineering such as robot and digital optimal control. Now everything around us depends on digital systems for its operation and function. Maple programming contains important packages such as the linear algebra package. Using this algebra package and some commands in Maple, we will build a simple procedure to compute the exact solution of linear discrete systems with time-varying coefficients and containing forward operators and forward difference operators. First, we will find the closed form of a discrete linear system of first order, and then generalized to higher orders. Then we build the algorithm and turn it into a maple procedure. Finally, to illustrate our procedure, we will give some examples with different types of operators. By entering the system parameters into each example, we will obtain the exact solution graph.