syms x y
jacobian([x^2+2*y-8,x+y^2-y-3],[x,y])