In this work the discretization of the Hénon-Heiles system obtained by applying the Monaco and Normand-Cyrot method is investigated. In order to obtain dynamically valid models, several approaches covering from the choice of terms in the difference equation originated from the discretization process to the increase of the discretization order are analyzed. As a conclusion it is shown that discretized models that preserve both the symmetry and the stability of their continuous counterpart can be obtained, even for large discretization steps.