Propagation of a two-dimensional spatio-temporal electromagnetic beam wave is analysed. In parabolic (paraxial) approximation the exact analytical results for a spatio-temporal Gaussian impulse can be obtained. For solution of the full wave equation the numerical simulation has to be used. The various facets of this simulation are discussed here.