Problem:
Find the curl of the vector field $f(x,y,z)\uvec{i} + g(x,y,z)\uvec{j} + h(x,y,z)\uvec{k}$ where $f(x,y,z) = \frac{1}{\left(x^3+z\right)^2}$, $g(x,y,z) = y z$, and $h(x,y,z) = \sqrt[3]{x^3-z}$
from sympy import *
from sympy.vector import *

C = CoordSys3D('C')
x, y, z = C.x, C.y, C.z
f = (1/((x**3+z)**2))
g = y*z
h = y*z
print(curl(f * C.i + g * C.j + h * C.k))

