Solution (source code)

= Solution

$\nabla\times F=(3y^2,3x^2,3x^2+3y^2)$, whose divergence is $0$.