Divergence:
div(F) = ∂(3x²y² + 2z)/∂x + ∂(3x² + y + yz)/∂y - ∂(3x²z² + 2yx + xyz)/∂z
div(F) = 6xy² + 1 + z - 6x²z - xy
Curl:
curl(F) = [∂(3x² + y + yz)/∂z + ∂(3x²z² + 2yx + xyz)/∂y] i
… … … … … - [∂(3x²y² + 2z)/∂z + ∂(3x²z² + 2yx + xyz)/∂x] j
… … … … … + [∂(3x²y² + 2z)/∂y - ∂(3x² + y + yz)/∂x] k
curl(F) = (-xz - 2x - y) i + (6xz² + yz + 2y + 2) j + (6x - 6x²y) k