  The existence and uniqueness of the stationary distribution of the numerical
solution generated by the stochastic theta method is studied. When the
parameter theta takes different values, the requirements on the drift and
diffusion coefficients are different. The convergence of the numerical
stationary distribution to the true counterpart is investigated. Several
numerical experiments are presented to demonstrate the theoretical results.
