  We follow the time evolution of nonabelian gauge bosons from far-from-equilibrium initial conditions to thermal equilibrium by numerically solving an effective kinetic equation that becomes accurate in the weak coupling limit. We consider initial conditions that are either highly overoccupied or underoccupied. We find that overoccupied systems thermalize through a turbulent cascade reaching equilibrium in multiples of a thermalization time $t\approx 72./ (1-0.12\log \lambda)/\lambda^2 T$, whereas underoccupied systems undergo a "bottom-up" thermalization in a time $t\approx (34. +21. \ln(Q/T))/ (1-0.037\log \lambda)(Q/T)^{1/2}/\lambda^2 T$, where $Q$ is the characteristic momentum scale of the initial condition. We apply this result to model initial stages of heavy-ion collisions and find rapid thermalization roughly in a time $Qt \lesssim 10$ or $t\lesssim 1$ fm/c. 
Approach to equilibrium in weakly coupled nonabelian plasmas