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) = \log \left(\frac{y}{x^5}\right)$, $g(x,y,z) = \sin (y)$, and $h(x,y,z) = \sinh \left(x^5+z\right)$
from sympy import *
from sympy.abc import x, y, z

f = log((y/(x**5)))
g = sin(y)
h = sinh(x**5+z)
X = Matrix([f, g, h]) 
print(X.jacobian((x, y, z)))
