We formulate a numerical procedure to calculate Hawking radiation during non-equilibrium black hole formation. The procedure is applied to a static string in thermal AdS and it is shown that for an arbitrary initial state, the final state is an equilibrated heavy quark string. The fluctuations in the quark string are transmitted from the horizon to the boundary leading to Brownian motion in the boundary theory.