  We categorify the RK family of numerical integration methods (explicit and
implicit). Namely we prove that if a pair of ODEs are related by an affine map
then the corresponding discrete time dynamical systems are also related by the
map. We show that in practice this works well when the pairs of related ODEs
come from the coupled cell networks formalism and, more generally, from
fibrations of networks of manifolds.
