Problem:
Find the jacobian of the vector valued function $\begin{pmatrix} f(x,y,z)\\ g(x,y,z) \\ h(x,y,z)\end{pmatrix}$ where $f(x,y,z) = \tanh ^{-1}\left(\frac{z^4}{x}\right)$, $g(x,y,z) = \sqrt[3]{x y}$, and $h(x,y,z) = y+z^4$
from sympy import *
from sympy.abc import x, y, z

f = atanh((z**4)/x)
g = cbrt(x*y)
h = y+z**4
X = Matrix([f, g, h]) 
print(X.jacobian((x, y, z)))
