Solution (source code)

= Solution

Differentiate the three components with respect to their matching <Cartesian coordinates>:
$$
\boxed{\nabla\cdot\mathbf F=y^2z^2+4y^3+x^3y^2}.
$$
The $3xz^2$ term contributes nothing to the $y$ <derivative> because it is independent of $y$.