The interaction between Gaussian beam and indefinite anisotropic metamaterial (AMM) slabs are simulated with auxiliary differential equation method in finite-difference time-domain method (FDTD). The directions of wave vector and energy flow densities in three types of AMM slabs are investigated. The field intensity distributions and phases of TE-polarized wave for the interaction between Gaussian beam and AMM slabs are given. Numerical results show that the sign of x component of relative permeability tensor of AMM determines the refraction is regular or anomalous, whereas the sign of z component of relative permeability tensor of AMM determines the energy flow is positively refracted or negatively refracted. Positive or negative Goos-Hanchen shift associated with Total Cutoff media are also shown.