We consider the dynamics of N bosons which are initially in a Bose-Einstein condensate and interact via pair interactions in the mean-field scaling regime. We derive a sequence of N-body functions which approximate the true many-body dynamics in norm to arbitrary precision in powers of 1/N. Joint work with Natasa Pavlovic, Peter Pickl, and Avy Soffer.