Problem:
Find the divergence 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) = \cos \left(y^2 z^2\right)$, $g(x,y,z) = \sqrt[3]{y^2-z^2}$, and $h(x,y,z) = \left(x y^2+z^2\right)^{3/2}$
from sympy import *
from sympy.vector import *

C = CoordSys3D('C')
x, y, z = C.x, C.y, C.z
f = cos(y**2*z**2)
g = cbrt(y**2-z**2)
h = (x*y**2+z**2)**(3/2)
field = f * C.i + g * C.j + h * C.k
print(divergence(field))
