The jacobian command computes the Jacobin matrix M where A is a vector of functions and X is a vector of variables. The matrix has elements M[i][j] equal to the derivative of A[i] with respect to X[j].
A = [x^2*y,sin(y)]; jacobian(A,[x,y]); [[2*x*y, x^2], [0, cos(y)]]