Solution (source code)

= Solution

Integrating the first component and checking the remaining <partial derivatives> gives
$$
\boxed{\phi_3=x^2yz-\cos x+C,\qquad\mathbf F_3=\nabla\phi_3.}
$$
This explicit <scalar potential> proves that the field is <conservative>.