Solution
= Solution
Integrating the first component in $x$ gives $\phi=x^3y+h(y,z)$. Matching the second requires $h_y=-2yz^2$, so $h=-y^2z^2+k(z)$. Matching the third gives $k'=0$. Thus the field is <conservative>, with <scalar potential>
$$
\boxed{\phi=x^3y-y^2z^2+C.}
$$
Direct differentiation verifies $\boldsymbol F_1=\nabla\phi$ on all of $\mathbb R^3$.