Se describe un método especial para la estimación de errores sistemáticos y accidentales acumulados en un proceso de integración numérica de un sistema de ecuaciones diferenciales ordinarias. El método se ha aplicado con éxito, en particular en el caso de la integración numérica de las ecuaciones de variación de los parámetros de una órbita. Se utilizará también en un estudio comparativo de los métodos más usuales en el cálculo de órbitas.