Over the last decades Set Oriented Numerical Methods have been developed for the analysis of global behavior of dynamical systems. Instead of computing single long term trajectory, these method rely on the information obtained from several short term trajectories and allow to compute invariant sets and invariant manifolds. This talk concerns the basic technique and implementation of Set Oriented numerics with some application to astrodynamics.