Abstract
We discuss various numerical approaches for studying the chaotic dynamics of multidimensional Hamiltonian systems, focusing our analysis on the chaotic evolution of initially localized energy excitations in the disordered Klein-Gordon oscillator chain in one spatial dimension.