  The Hermite pseudospectral method is applied to solve the Navier-Stokes
equations on a two-dimensional infinite domain. The feature of Hermite function
allows us to adopt larger time steps than other spectral methods, but also
leads to some extra computation when the stream-function is calculated from the
vorticity field. The scaling factor is employed to increase the resolution
within the region of our main interest, and the aliasing error is fully removed
by the 2/3-$rule$. Several traditional numerical experiments are performed with
high accuracy, and some related future work on physical applications of this
program is also discussed.
