  We present an algorithm to analyze numerically the bounce solution of
first-order phase transitions. Our approach is well suited to treat phase
transitions with several fields. The algorithm consists of two parts. In the
first part the bounce solution without damping is determined, in which case
energy is conserved. In the second part the continuation to the physically
relevant case with damping is performed. The presented approach is numerically
stable and easily implemented.
