  We present the theoretical basis and sample Monte Carlo data for the YFS
exponentiated calculation of $e^+e^- \to Z Z \to f_1\bar f_1 + \bar f_2 f_2$ at
and beyond LEP2 energies, where the left-handed part of $f_i$ is a component of
an $SU_{2L}$ doublet, $i=1,2$. The calculation is performed for both SM
couplings and for anomalous $ZZV$ triple gauge boson couplings in the
conventions of Hagiwara et al.. Our formulas, which are gauge invariant, are
illustrated in a proto-typical YFS Monte Carlo event generator YFSZZ.
