Numerical approximation of center, stable and unstable manifolds of multiscale/stochastic systems