  We have recently constructed a numerical code that evolves a spherically
symmetric spacetime using a hyperbolic formulation of Einstein's equations. For
the case of a Schwarzschild black hole, this code works well at early times,
but quickly becomes inaccurate on a time scale of 10-100 M, where M is the mass
of the hole. We present an analytic method that facilitates the detection of
instabilities. Using this method, we identify a term in the evolution equations
that leads to a rapidly-growing mode in the solution. After eliminating this
term from the evolution equations by means of algebraic constraints, we can
achieve free evolution for times exceeding 10000M. We discuss the implications
for three-dimensional simulations.
