  This paper presents both a numerical method for general relativity and an
application of that method. The method involves the use of harmonic coordinates
in a 3+1 code to evolve the Einstein equations with scalar field matter. In
such coordinates, the terms in Einstein's equations with the highest number of
derivatives take a form similar to that of the wave equation. The application
is an exploration of the generic approach to the singularity for this type of
matter. The preliminary results indicate that the dynamics as one approaches
the singularity is locally the dynamics of the Kasner spacetimes.
