Numerical Computation of the Stable and Unstable Manifolds of