Numerical study of ergodicity for the overdamped Generalized Langevin Equation with fractional noise