An outstanding problem in numerical analysis is to justify long-time simulations of molecular systems. Typically, these systems are chaotic and numerically computed trajectories are accurate for only short periods of time. Nevertheless, qualitative features of computed trajectories are often quite accurate. I will discuss some recent work on this problem. The main innovation is the introduction of a weak ergodicity that suffices for the purposes of numerical simulation but is easier to work with than the standard mathematical definition.