Liouville's Theorem: The volume of any closed surface in phase space remains constant as that surface moves through phase space. This is a consequence of trajectories never crossing.

⇒∇⋅v=∇⋅z˙=0

Apply divergence theorem to Hamilton's equations:

z˙=(∂H∂p,−∂H∂q)∇⋅z˙=∂∂q∂H∂p+∂∂p(−∂H∂q)=0