Suppose a vector-valued diffusion process X satisfies the stochastic differential equation $$dX_t = b(X_t)dt + \sigma(X_t) dW_t,$$ in which $W$ is a Brownian motion and $b,\sigma$ are such that strong existence and uniqueness of solution to the SDE hold. Assume also that $X$ has a stationary distribution. Sufficient condition for the existence of a stationary distribution can be found in, e.g. Khasminskii, Stochastic Stability of Differential Equations. In a recent book by Mao and Yuan, Stochastic differential equations with Markovian switching, a Euler Mariamman method is proposed to numerically approximate the stationary distribution of $X$. But there is no mentioning on the efficiency of the algorithm, say, the convergence rate. What will be the convergence rate? And are there any other methods in the literature available? Many thanks for your comments and suggestions.