Mathematica can be used to calculate the Jacobian matrix in exact form for any given system of equations. There is currently no built-in definition, but a simple construction using the partial differentiation operator D and the generalised outer product Outer is possible [15]:
JacobianMatrix[f_List,x_List] := Outer[D,f,x]