\phi(x, y, z) =
\int (xy - z)dz = xyz - \frac{1}{2}z^2 + g(x, y)