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

f = (1/((z**2)**(3/2)))
g = tan(z**2)
h = z**2-y
X = Matrix([f, g, h]) 
print(X.jacobian((x, y, z)))
