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(\frac{y}{x^4}\right)^{3/2}}$, $g(x,y,z) = \sqrt[3]{y+z^2}$, and $h(x,y,z) = \frac{1}{z^2}$
from sympy import *
from sympy.vector import *

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

