In this paper, we solve the exact output tracking problem For square nonlinear discrete-time systems through stable inversion. This stable inversion problem involves finding a bounded solution of a dynamical system driven by a bounded signal. Such a solution is given by a difference representation formula, whose existence, under the proper assumptions, can be proven using a Picard process.