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}{x^2}$, $g(x,y,z) = (x+y)^3$, and $h(x,y,z) = e^{\frac{x+y}{z}}$
import math

from sympy import *
from sympy.abc import x, y, z

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