Jacobian(L:LIST):MAT where L is a list of polynomials.
Use R ::= Q[x,y]; L := [x-y,x^2-y,x^3-y^2]; Jacobian(L); Mat([ [1, -1], [2x, -1], [3x^2, -2y] ]) -------------------------------