This group considered the stability of matrix Lie group valued stochastic differential equations, dXt = AXtdt + BXtdt. Random dynamical systems such as this arise in many applications (e.g., oceanic turbulence, helicopter blade motion, light in random channels, wireless networks) in which stability is of practical and theoretical concern. The stability of the zero solution, Xt = 0, is determined by the top Lyapunov exponent. But in practice, analytic calculations of the Lyapunov exponent are often impossible, so time discrete approximations and simulations are necessary. Oceledet’s famous multiplicative ergodic theorem ensures that the Lyapunov exponent of Xt is almost surely constant. But the Lyapunov exponent of a time discrete approximation is itself a random variable. The mean of this random variable has been studied. The group provided estimates on the variance, distribution and rate of convergence in certain numerical approximation methods. In addition, they expanded on and provided simulations for recent results regarding the top Lyapunov exponent of certain Lie group valued SDEs.