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) = \tan \left(x+z^3\right)$, $g(x,y,z) = \tanh \left(y z^3\right)$, and $h(x,y,z) = \cos ^{-1}\left(\frac{z^3}{y}\right)$
from sympy import *
from sympy.vector import *

C = CoordSys3D('C')
x, y, z = C.x, C.y, C.z
f = tan(x+z**3)
g = tanh(y*z**3)
h = acos((z**3)/y)
field = f * C.i + g * C.j + h * C.k
print(divergence(field))
