An exact analytical solution is obtained to the reflection and refraction of a Gaussian beam by a biaxial anisotropic slab. The electromagnetic waves in different regions are expanded in terms of rectangular vector wave functions or represented as superposition of eigen plane waves, and the unknown expansion coefficients or angular spectrum amplitudes are determined by virtue of the boundary conditions. For a localized beam model, the normalized field intensity distributions are calculated, and the propagation characteristics are discussed briefly.