We describe an effective computational procedure of obtaining estimate for the global error in the long-term numerical integration with one-step embedded explicit Runge-Kutta methods applied to an ODE system. The paper relies on the existence of an asymptotic expansion for the global error. An example is given using a particular ODE system modeling the behavior of a heterogeneous catalytic reaction with complex dynamics.
