A result of BEYN on the replication of the phase portrait of an ordina
ry differential equation in the neighborhood of a hyperbolic steady st
ate by a one-step numerical method is generalized under appropriate as
sumptions to random differential equations with hyperbolicity defined
in terms of the multiplicative ergodic theorem of Oseledets. The proof
is complicated by the need to compare cocycles rather than underlying
vector fields and the use of random norms to obtain uniform estimates
. As an example we consider small random perturbations of a linear ord
inary differential equation with a hyperbolic null solution.