A new method for computing the unsteady motion of a water surface, including the overturning of water waves as they break, has been developed. It is based on a Cauchy theorem boundary integral for the evaluation of multiple time derivatives of the surface motion. The numerical implementation of the method is efficient, accurate and stable.