Como \( \mathbb{R}^3 \) es simplemente conexo entonces \( F\cdot \,d r \) será exacto si y solo si es cerrado, es decir, si y solo si \( d(F\cdot \,d r)=0 \), entonces
\[
\begin{align*}
d(F\cdot d r)&=d((z+x^2)\,d x+2yz^2\,d y+g\,d z)\\
&=-dx\,dz-4yz\,dy\,dz+\partial _xg\,dx\,dz+\partial_yg\,dy\,dz\\
&=0\\
\iff &\partial _xg=1\,\land\, \partial _yg=4yz
\end{align*}
\]
Entonces tenemos una ecuación diferencial parcial. Al integrar tenemos que
\[
g(x,y,z)=g(x,0,z)+2y^2z\,\land\, g(x,y,z)=g(0,y,z)+x
\]
De ahí vemos que una posible solución al problema es \( g(x,y,z)=x+2y^2z \).∎