This paper deals with the numerical continuation of invariant manifolds, regardless of the restricted dynamics. Typically, invariant manifolds make up the skeleton of the dynamics of phase space. Examples include limit sets, co-dimension 1 manifolds separating basins of attraction (separatrices), stable/unstable/center manifolds, nested hierarchies of attracting manifolds in dissipative systems and manifolds in phase plus parameter space on which bifurcations occur. These manifolds are for the most part invisible to current numerical methods. The approach is based on the general principle of normal hyperbolicity, where the graph transform leads to the numerical algorithms. This gives a highly multiple purpose method. Examples of computation...