Ornstein-Uhlenbeck is a path-dependent SDE. The negative \(\mu\) in its standard stochastic calculus is representing mean-reversed force or restoring force. The disadvantages of Euler-Maruyama method have problems of weak convergence and explosion with time evolution for path-dependent SDE. The above five figures show that when negative \(\mu\) is growing larger and larger, the solution of approximation with Euler-Maruyama method will explode with time evolution. Therefore, Euler-Maruyama method is not the best method to approximate non-stationary path-dependent SDE.