from sympy import *
from sympy.abc import *
x,y=symbols('x,y', real=True)
f1=-y
f2=x - 3*y*(1-x**2)
X=Matrix([f1,f2])
Y=Matrix([x,y])
print("Variables=",Y)
print("function=",X)
print("Jacobian=",X.jacobian(Y))
