In this thesis, we discuss the numerical approximation of random periodic solutions (r.p.s.) of stochastic differential equations (SDEs) with multiplicative noise. We prove the existence of the random periodic solution as the limit of the pull-back flow when the starting time tends to ∞ along the multiple integrals of the period. As the random periodic solution is not explicitly constructible, it is useful to study the numerical approximation. We discretise the SDE using the Euler-Maruyama scheme and modified Milstein scheme. Subsequently we obtain the existence of the random periodic solution as the limit of the pull-back of the discretised SDE. We prove that the latter is an approximated random periodic solution with an error to the exact one at the rate of √Δ t in the mean-square sense in Euler-Maruyama method and Δ t in the modified Milstein method. We obtain the weak convergence result in infinite horizon for the approximation of the average periodic measure.