Para la aproximación de sistemas rígidos de EDOs que surgen de la cinética química, los integradores implícitos emergen como buenos candidatos. Este artículo propone un enfoque variacional para este tipo de sistemas. Además de presentar la técnica, presentamos sus propiedades más básicas y probamos su rendimiento numérico a través de algunos experimentos. La principal ventaja con respecto a otros métodos implícitos es que nuestro enfoque tiene una convergencia global. Los otros enfoques necesitan garantizar la convergencia del esquema iterativo utilizado para aproximar las ecuaciones no lineales asociadas que aparecen por la implicitud. Tenga en cuenta que estos métodos iterativos, para estas ecuaciones no lineales, tienen cuencas de atracción limitadas.